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This  work  extends  single  inpu£ -output  linear  time  invariant  minimum-phase 
’'’quantitative  feedback  synthesis^'to  two  new  complex  plant  structures  with  in- 
ternal sensing  points.  One  is  the  triangular  structure.  The  second  consists 
of  parallel  branches,  each  with  cascaded  sections.  Due  to  uncertainty,  the  plan 
parameters  are  elements  of  given  sets.  The  system  response  must  satisfy  specif i 
time  or  frequency  domain  tolerances.  The  basic  problem  is  how  to  divide  the 
feedback  burden  among  the  available  loops  so  as  to  minimize  the  net  rms  effect 
at  the  plant  input,  of  the  various  sensor  noise  sources. 


Abstract 


This  work  extends  single  input-output  linear  time  invariant 
minimum-phase  "quantitative  feedback  synthesis"  to  two  new  complex 
plant  structures  with  internal  sensing  points.  One  is  the  triangular 
structure.  The  second  consists  of  parallel  branches,  each  with 
cascaded  sections.  Due  to  uncertainty,  the  plant  parameters  are 
elements  of  given  sets.  The  system  response  must  satisfy  specified 
time  or  frequency  domain  tolerances.  The  basic  problem  is  how  to 
divide  the  feedback  burden  among  the  available  loops  so  as  to  minimize 
the  net  rms  effect  at  the  plant  input,  of  the  various  sensor  noise 
sources . 

Frequency-response  formulations  are  presented  which  provide  a 
deep  understanding  of  the  trade-off  among  the  feedback  loops.  One 
vital  feature  is  "free  uncertainty",  wherein  a loop  optimized  to  cope 
with  uncertainty  , may  in  fact  for  some  frequency  ranges , handle 
uncertainty  U >>  . A second  is  "bandwidth  propagation",  wherein 

the  loops  take  turns  in  dominating  the  design  over  the  frequency  range. 
Together,  they  locate  the  frequency  regions  in  which  the  respective 
loops  dominate,  and  the  key  trade-off  parameters  among  them.  1 Design 
Perspective"  then  enables  the  designer  to  very  rapidly  find  a close 
approximation  to  the  precise  design  based  on  any  choice  of  these 
parameters.  Numerous  design  examples  with  very  large  uncertainty, 
illustrate  the  design  procedures  and  the  advantages  of  multiple-loop 
design. 
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CHAPTER  1 


STATEMENT  OF  PROBLEM  AND  PRELIMINARY  BACKGROUND 
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1 . 1 Introduction 

This  work  deals  with  the  problem  of  making  a system  perform 
satisfactorily  despite  uncertainty,  in  the  following  context: 

Equipment  has  been  assembled  by  specialists  in  the  area  of 
concern,  in  order  to  achieve  certain  objectives,  for  example,  an 
airframe  and  engine  for  accomplishment  of  certain  aeronautical 
objectives,  a chemical  plant  for  production  of  certain  chemicals, 
etc.  This  assembly  of  equipment  is  denoted  as  the  plant.  The  plant 
has  the  ability,  the  muscle  so  to  speak,  to  achieve  the  objectives. 
However,  it  does  not  have  the  accuracy  needed.  This  is  manifested 
by  uncertainty  in  the  parameters  of  the  mathematical  relations  des- 
cribing the  plant.  For  example  in  Figpre  1.1,  suppose  the  relation 


Plant 


Figure  1-1.  Plant 


between  plant  output  y and  input  x is  given  by  a linear  time- 
invariant  (Iti)  ordinary  differential  equation 

+ A y(n  ♦ ...  + A y *■  K(x*m^+B  x*1"-1^  + ...  + B x] 
l n 1 m 


- • -- 


i 


(1.1-1) 


i.e.,  the  , B ^ are  constants.  There  are  q physical  para- 
meters k,,...,k  and  K , A.  , B.  are  functions  of  these  para- 
1 q i j 

meters.  The  values  of  these  parameters  are  not  known  precisely,  but 

1 2 

it  is  known  that  they  lie  within  certain  bounds  k^  € lk^,k^]  . 

Hence,  each  possible  parameter  vector  k , each  combination  of 
values  i - l,...,q  , gives  a different  plant  transfer  function, 
generating  a set  = (P(s) } of  possible  plant  transfer  functions. 

Such  a formulation  of  the  uncertainty  problem  may  appear  naive 
because  one  might  argue  that  often  the  parameter  values  change  with 
time-giving  uncertain  linear  time-varying  relations,  because  the 
rate  of  variation  is  uncertain.  Also,  the  tti  description  is 
usually  an  approximation  of  a nonlinear  relation.  We  are  really 
assuming  £ti  relations  with  the  above  uncertainty  form,  in  order 
to  be  able  to  rigorously  use  Laplace  transforms  and  frequency  res- 
ponse methods.  However,  it  has  been  rigorously  proven  [HI]  that 
a large  class  of  practical  uncertain  linear  tine-varying  plant  problems 
is  reducible  to  the  above  tti  uncertainty  form,  and  even  uncertain 
nonlinear  time-varying  plants  to  a certain  extent  [H2],  can  be  so 
reduced.  Hence  the  above  modelling  is  highly  justifiable,  because  the 
techniques  developed  for  the  tti  synthesis  problem,  can  be  rigorously 
used  for  these  classes  of  linear  and  nonlinear  time-varying  synthesis 
problems.  Finally,  one  must  begin  somewhere  with  the  development  of  a 
scientific  synthesis  theory  for  uncertain  systems,  and  "scientific" 
implies  quantitative".  The  tti  case  is  obviously  where  to  start. 

The  objective  is  to  achieve  certain  apriori  specified  per- 
formance objectives  V P € . Since  the  overall  system  is  to  be 


.1 

i 


tti  , it  can  be  characterized  by  its  response  to  any  input,  and  the 


step  response  is  very  popular  because  it  combines  within  it  both  the 


fastest  hind  of  input  (an  abrupt  change)  and  the  slowest  (no  change) 


Time  domain  specifications  are  reasonable  in  many  cases,  as  in 


Figure  1.2a,  the  step  response  is  to  be  inside  the  bounds 


with  additional  bounc..  of  similar  nature  on  the 


upper-^ 

bound/ 

_ b,  / 


lower  bound 


Figure  1-2 (a)  . Time  domain  step  response  specification 


first  and  perhaps  higher  derivatives.  Our  design  technique  is  in 


the  frequency-domain,  so  we  must  translate  such  t-domain  bounds  into 


equivalent'1  u-domain  bounds  on  the  system  frequency  response 


T(jw)  . If  the  system  is  minimum-phase  (H3) 


we  restrict  ourselves  here  to  such  systems.  This  translation  is,  as 


of  this  date,  an  engineering  art  rather  than  a science.  Advice  on 


how  to  translate  is  scattered  in  the  literature  [H3,  Kl,  SI].  Very 


good  results  have  been  obtained  with  only  moderate  effort.  We  shall 


0 


Figure  1-2 (b) . Frequency  domain  specification. 


assume  in  this  work  that  the  translation  has  already  been  done.  It  is 
worth  noting  that  it  has  been  shown  [H2]  that  for  minimum-phase 
systems,  time-domain  specifications  on  the  step  response  and  on  its 
derivatives  of  the  following  nature 

bj^  (t)  j cU)  (t)  $ bjx>  (t)  , i - 0,1, ...  ,n  , t € 10, »)  (1.1-2) 

can  always  be  satisfied  by  means  of  to-domain  bounds  of  the  following 
nature 

B2(w)  $ | T ( j o> ) | < B^u)  . (1.1-3) 

In  our  work,  the  bounds  on  system  performance  will  have  this  form. 
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1.1.1.  Previous  work. 

The  quantitative  aspect  of  our  work  cannot  be  over-emphasized. 

The  sensitivity  reduction  capability  of  feedback  is  very  well-known. 

Hundreds  of  books  and  thousands  of  papers  have  been  written  on  the 

subject,  but  the  number  of  these  which  are  quantitative  in  nature  is 

extremely  small,  i.e.,  with  uncertainty  bounds  and  performance 

bounds  explicitly  included  in  the  problem  statement.  It  is  as  if 

the  mere  use  of  a feedback  configuration  around  the  uncertain  plant, 

suffices  to  scare  it  into  docile  behavior.  In  the  vast  majority  of 

the  techniques  the  uncertainty  is  completely  ignored,  and  there  are 

none  or  extremely  crude  performance  specifications.  One  presumably 

emerges  with  the  same  design  whether  the  parameter  uncertainty  is 

x%  or  1000x%,  and  irrespective  of  whether  the  bounds  , B2  in 

Figure  1-2 (b)  are  narrow  or  wide  apart.  There  is  no  concern  with 

the  'cost  of  feedback'  - which,  aside  from  the  sensors,  lies  in  the 

bandwidth  of  the  loop  transfer  function,  and  little  concern  with  the 

extremely  important  matter  of  sensor  noise  [see  Sec.  1.2].  These 

points  have  been  emphasized  in  [H4,  H12]. 

* 

Our  work  follows  closely  in  the  tradition  of  'quantitative 
synthesis'  recently  established  [H3-6,  Wl],  To  appreciate  the 
present  work,  it  is  important  to  be  aware  of  the  highlights  of  this 
previous  work.  Quantitative  synthesis  was  first  developed  for  a 
plant  with  only  one  variable,  the  plant  output  c(t)  in 
Figure  1-3,  available  for  feedback  [H3] . The  system  command  input 
r(t)  was  also  assumed  accessible,  so  the  processing  of  these  two 
signals  provides  two  independent  compensation  functions  to  the 
designer.  An  infinitude  of  canonical  two-degr ee-of-freedom 


structure  may  be  used  [H7] . The  design  procedure  developed  in  [H3] 
used  Figure  1-3  (a),  but  suppose  the  sensor  transfer  function  is 
H(s)  , then  one  can  use  Figure  1-3 (b)‘,  letting  G^H  (of 
Figure  1-3 (b) ) = G (of  Figure  1-3 (a)),  in  order  to  have  the  same 
loop  transmission  function  L(s)  = GP  = GjPH  , and  F^G^  = FG  in 
order  to  have  the  same  system  transfer  function 


J 

0 


found  that  an  independent  sensor  measuring  c(t)  (e.g.,  a tacho- 
meter in  a position  servo)  gives  less  noise  than  the  differentiation 
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of  a position  .sensor,  so  both  sensors  are  used,  as  in  Figure  1-4, 
with  the  two  sensor  transfer  functions  , and  say  the 

structure  in  Figure  1-4  is  used.  This  is  a two-loop  structure 


Figure  1-4.  2-loop,  2-D.O.F.  structure. 


physically,  but  in  terms  of  fundamental  feedback  design  it  is  a two- 
degree-of- freedom  system,  so  the  quantitative  design  theory  of 
Figure  l-3(a)  is  used,  giving  G and  F . It  is  required  that  the 
loop  transmission  around  F , be  the  s^ame  in  both  cases,  i.e., 


I 

L - PG  (Fig.  1-3 (a))  - 


and 


1 


so 


T • F TO?  ‘ "b  i.Tio I 


0 - vvw  • w, 

I 

M^,  11^  are  known,  so  o,Je  must  decide  how  to  split  G ■ G (Hj+G^Hj) 
between  G^  and  G^  . [this  is<  done  by  considering  the  effect  of 
Honsor  noise  , N2  it  the  plant  input, 

vvKw 


(1.1-5) 


(1.1-6) 


(1.1-7) 


XN(lw)  " UPlC^+G^n 

\ 


<)  \ 


(1.1-8) 


mmmm 
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given  that 

G - jr(ju>)  - Ga(H1  + H2Gb)  (1.1-9) 

is  fixed  by  the  quantitative  design  technique  of  (H31 . 

The  objective  is  to  minimize  / IXjJ  dw  , subject  to  the 

o 

above  constraint.  This  is  a straightforward  optimization  problem 
which  can  be  solved  outside  the  realm  of  quantitative  feedback 
synthesis.  The  latter  only  provides  the  design  with  the  feedback 
loop  transmission  ( L ) needed  around  the  plant,  and  the  prefilter 
(F)  needed  to  process  the  command  input  r(t)  . The  state-of-the- 
art  in  sensors  and  in  filter  synthesis  determines  how  L and  F 
are  to  be  realized.  In  fact,  in  the  above  context  one  might  con- 
sider use  of  an  accelerometer  in  a 3-loop  feedback  structure.  But 
from  our  point  of  view  the  structure  remains  that  of  a two-degree- 
of-freedom  system  and  we  shall  continue  to  associate  the  latter 
with  a single-loop  system. 


1 . 2 Review  of  Two-Degree-of-Freedom  Quantitative  Design 
Theory 

GP 

Figure  l-3(a)  is  used  with  T ■=  F -y—  . It  is  assumed  that 
the  compensation  network,  whose  power  level  can  be  very  low  (as 
the  plant  contains  the  power  elements),  can  be  constructed  with 
negligible  uncertainty  in  their  transfer  functions.  Hence,  due  to 
the  uncertainty  in  P , 

GP  I 

Un  T * Hn  - A In  — # L - GP  (1.2-1) 

and 

L(ju>) 

1+1.  ( j w) 


A In  l T ( Jw)  t - A tn 


(1.2-2) 
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Given  that  the  specifications  require  that  Af.n|T(ju>)|  ( 61<lb  at  “i  in 

Figure  1-2 (b),  what  are  the  resulting  constraints  on  L(ju^)  ? It 

is  convenient  to  pick  a "nominal"  plant  Pq(s)  , and  derive  the 

bounds  on  the  resulting  "nominal"  loop  function  L *■  P G . These 

o o 

bounds  can  be  found  by  means  of  a digital  computer,  but  it  is  very 

useful  for  insight  to  see  it  done  on  the  Nichols'  chart 

(logarithmic  complex  plane  with  abscissa  in  degrees,  ordinate  in 
decibels  - 20  log1Q) • The  procedure  is  illustrated  for  the  case 

P(S)  " a(s>a)  S k 6 11,101  » a 6 11,101  ' (1-2-3) 

and  say  k-1  , a - 10  are  chosen  as  nominal,  giving  Pq  - 10/($  (s+10)}  . 
At  hi  ■ 2 rps  , P(j2)  lies  within  the  boundaries  given  by  ABCD  in 
Figure  1-5.  Since  In  L « In  G + (n  P , the  pattern  outlined  by 


|U 


AUCD  may  be  translated,  but  not  rotated,  on  the  Nichols'  chart 


the  amount  of  translation  being  given  by  the  value  of  tn  G(j2)  . 
For  example,  if  a trial  design  of  L(j2)  corresponds  to  the  tem- 
plate P(j2)  at  A'B'C'D'  in  Figure  1-5,  then 

l°u2'U  • lu’2,ldb  - lp,sj>ldb 


- (-2.0)  - (-13.0)  = 11.0  db 


(1.2-4) 


Arg  G(j2)  = Arg  L(j2)  - Arg  P(j2) 

= (-60°)  - (-153.4°)  = 93.4°  . 


(1.2-5) 


1.2.1.  Bounds  on  L(jai)  in  the  Nichols'  chart. 

The  templates  of  P(ju>)  are  manipulated  to  find  the  position 
of  L(jw)  which  results  in  the  specifications  of  Figure  1-2 (b)  on 
fn  | T ( jou)  | being  satisfied.  Taking  the  u>  = 2 template,  one 
tries,  for  example,  positioning  it,  as  shown  in  Figure  1-5,  at 
A'B'C'D’  . Contours  of  constant  fnlL/(l+L)|  are  available  on  the 
Nichols'  chart.  Using  these  contours,  it  is  seen  that  the  maximum 
change  in  fn|  L/(l+L)  | which  from  Fig.  1.2-2  is  the  maximum  change  in 
tn  | T | is,  in  this  case,  very  closely  (-0.49)  - (-5.7)  = 5.2  db  , 
the  maximum  being  at  point  C'  , the  minimum  at  point  A'  . 

Suppose  that  the  specifications  tolerate  a change  of  6.5  db  at 
w = 2 , so  the  above  trial  position  of  |LQ(j2) | is  in  this  case 
more  than  satisfactory.  The  template  is  lowered  on  the  Nichols' 
chart  to  A"B"C"D",  where  the  extreme  value  of  fn|L/(l+L)|  are 
at  C"  (-0.7  db)  , A" (-7.2  db)  . Thus,  if  Arg  LQ(j2)  = -60°  , 
then  -4.2  db  is  the  smallest  magnitude  of  Lc(j2)  which  satis- 
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fics  the  6.5  db  specification  for  A £n  |t|  . Any  larger  mag- 
nitude is  satisfactory  but  represents  over-design  at  that  frequency. 
The  manipulation  of  the  u ■ 2 template  is  repeated  along  a new 

vertical  line,  and  a corresponding  new  minimum  of  |l  ( j 2 ) | found. 
Sufficient  points  are  obtained  in  this  manner  to  permit  drawing  a 
continuous  curve  of  the  bound  on  Lq(j2)  , as  shown  in  Figure  1-5. 
The  above  is  repeated  at  other  frequencies,  resulting  in  a family 
of  boundaries  of  permissible  Lo(ju>)  . 

1.2.2.  Nature  of  the  bounds  on  L ( jio) . 

■ — i -i-  ■ i-—  — — - O ' *• 

A typical  set  of 'bounds  is  shown  in  Figure  1-6.  The  bounds 


-360°  -240°  -120°  0“ 

DEGREE 

Figure  1-6.  Typical  bounds  on  Lq  in  Nichols'  chart. 
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tond  to  mo vo  down  in  tho  Nichols'  chart  (become  less  onerous), 
obviously  because  as  u increases,  greater  change  in  | T < ju>)  | 
is  permitted,  as  in  Figure  1-2 (b).  It  is  in  fact  essential  that 
at  large  enough  u>  , the  uncertainty  in  |T(ju)|  (i.e.,  the  bounds 

on  | T ( ju»)  | ) be  greater  than  the  uncertainty  in  P(jw)  , because 
the  net  sensitivity  reduction  is  always  zero  in  any  practical 
system  as  was  long  ago  [Bl]  shown  by  Bode, 

/ *n  | ( jw)  |dw  - - / <n  |l+L(jw)|  d<*i  - 0 (1.2-6) 

o o 


T DT /T 

where  sp  “ alv?  is  tfle  sensitivity  function. 

U A 

In  the  above  example  as  w -*■  » , P , so 

A In  |p|  -»  A fn(Ka)  ■ 40  db  . Note  in  Figure  1-2 (b)  that  the 
permitted  A tn  |T(jw)|  40  db  as  w > 50  . Such  large  toler- 

ances on  | T ( jw) | at  large  m are  tolerable  because  |T(jw)|  is 
negligible  at  large  to  , e.g.,  if  |p(jw)|  can  change  at  most 
by  40  db  at  large  u>  but  | T ( jtu)  | changes  by  52  db  , who  cares 
if  this  52  db  change  is  from  |t| 

400  x 10  6 . In  return,  one  can  concentrate  the  sensitivity 
reduction  over  the  bandwidth  of  T(jw)  . Thus,  although  |p(jw)| 
in  this  region  varies  by  say  40  db  , |T(jw)|  may  be  controlled 
to  vary  by  only  4 db  , or  0.04  db  if  desired. 


, • 10  & to  | T | 

min  1 'max 


Universal  high-frequency  (UHF)  boundary. 

As  noted,  in  the  high-frequency  range  A tn  | T ( ju»>  | must 
realistically  be  allowed  to  be  x*  A tn  |r(jw) | , and  this  is 
reflected  in  the  bounds  on  LQ(iw)  tending  to  a very  narrow  pencil 
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in  the  arithmetic  complex  plane  (if  -♦  — — - as  s ■*  ®)  as  in 

° se 

Figure  1-7 (a)  and  as  in  1-7 (b)  in  the  Nichols'  chart.  In 


Figure  1-7 (b) , the  boundary  is  drawn  for  the  case  A fn  L - A In  k» 

20  db  , A in  | T ( ju>)  | - A in  |l/(1+L)|  $ 23  db  at  • However,  the 

resulting  peak  value  of  |l/(1+L)|  is  23  db  ■ 14.1  arithmetic  at 

k » k , indicating  a highly  under-damped  pole  pair  at  the 

corresponding  frequency  with  damping  ratio  £ * 0.034  , when 

k - * This  tremendous  peaking  does  not  appear  in  the  system 

response  to  the  command  inputs  R , because  it  is  filtered  out  by  the 

pre-filter  F in  Figure  1-3 (a).  But  the  system  response  to  a disturb- 

C -1 

ance  D in  Figure  1-3,  is  given  by  Td  ■ ^ ■ (1+L)  . Disturbance 

attenuation  generates  its  own  requirements  on  L , which  may  lead 


L L (OEGREE) 

Figure  1-7 (b) . Bounds  on  L(jw)  on  Nichol's  chart. 

to  more  stringent  bounds  on  L than  those  due  to  T(ju>)  . The 
final  contours  used  in  the  design  [H3]  must  be  the  most  stringent 
composite  of  the  two.  However,  even  if  D is  very  small,  it  is 
usually  certain  that  a peak  | | of  14.1  is  ihtolerable.  It  is 
reasonable  to  add  a requirement  |~j-|  f y a constant,  for  all 
a>  and  over  the  whole  range  of  P parameter  values.  The  resulting 
constraining  contours  denoted  by  B^  are  shown  in  Figure  1-7 (b) 


.1 

i 

I 


6 
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for  tho  cast;  A In  It  • 20  db  , and  fur  y *=  2.3,  3.5,  5 db  (all 
those  contours  arc  symmetrical  with  respect  to  the  vertical 
Arg  L “ -180“  on  the  Nichols'  chart).  If  y = 5 db  is  used,  then 
B(w  ) indicates  the  composite  contour  shown  in  Figure  1—7 (b) . For 
u > , | A T(ju>)  | increases  while  y remains  the  same,  so  that 

sooner  or  later  there  is  reached  a frequency  3 B(w)  «■  ( y ) , 

V . This  boundary  B^  is  called  the  "universal  high  frequency" 

(UHF)  boundary. 

1.2.4.  The  optimum  L(ju>)  . 

It  has  been  showy  [H6,  H8)  that  a realistic  definition  of 

optimum  in  the  tti  system  is  the  minimization  of  k , defined  by 
» 0 

lim  L(s)  ■ k s , where  e is  the  excess  of  poles  over  zeros 
s x* 

1 assigned  to  L(s)  . 

It  has  been  proven  [H6,  H8]  that  the  optimum  L lies  on  its 
boundary  B.^  at  each  or  and  that  such  an  optimum  exists  and  is 
unique.  Most  important  for  the  present  purpose,  is  that  in  signi- 
ficant plant  ignorance  problems  the  ideal  optimal  L has  the 

properties  shown  in  Figure  1-8,  i.e.,  over  a significant  range  it 

0 

follows  B^  along  UV  up  to  the  point  J at  which  it  abrqptly 
jumps  to  infinity  along  WW'W”  and  returns  on  the  vertical  line 
Y '£  , whose  phase  is  (-90*)  e . Such  an  ideal  L(jw)  is,  of  course, 
impractical.  A practical  suboptimum  L is  shown  in  Figure  1-8 
(see  Section  1.4.3,  for  more  discussion). 


jj 


u 


L L (degree) 
-100* 

i I — Upt_ 


Practicol  subopti  mum  L 


figure  1-8.  Bounds  on  L and  Optimum  L on  Nichols'  Chart 


Some  results  of  a numerical  design  example  are  shown  in 


Figure  1-9.  They  were  derived  for  the  following  problem 


1.2.5.  Numerical  example.  (Figure  l-3(a)) 


Plant: 


Plant  ignorance: 


3 S+B  * D * *3  * 

Performance  Specification:  Shown  in  Figure  1-2 (b) , were  originally 

derived  from  time  domain  bounds  [SI]. 


Disturbance  response: 


This  example  is  used  later  in  a demonstration  of  the  advantages  of 


multiple- loop  design.  The  L 


transfer  functions  are  given  in 


Appendix  1 of  this  chapter 


In  significant  plant  ignorance  problems,  there  is  a strong 


tendency  for  the  design  to  be  such  that  N 


so  highly  an*>lified  as  to  saturate  the  plant  input  at  X . The 


noise  response  function  is  (see  Figure  l-3(a)) 


range  where  |L  <<  1 


N represents  the  square  root  of  the  noise  power  spectrum. 

The  noise  response  of  the  numerical  design  example  of  the  last 


section  is  show  in  Figure  1-10.  Notice  that  the  noise  component 


in  Figure  1-2 (a),  is  most  important  in  the  high-frequency 


range  where  the  useful  command  and  disturbance  components,  due  to  R 


and  are  relatively  small,  rather  than  in  Mie  low  frequency  range 

where  the  latter  are  relatively  large.  This  is  further  enhanced 


by  the  fact  that  arithmetic  scales,  shown  in  Figure  1-11,  must  be 


• )dw  , ♦ * N : noise  power  spectrum  (1.3-2) 


Figure  l-9a.  Single  loop  LgQ  and  outer  loop 
of  the  numerical  example  in  Section  1.2.5. 
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Hence,  it  is  desirable  to  decrease  | L | vs  w , as  fast  as 
possible  in  the  high  frequency  range.  Even  a saving  which  is  small 
in  the  logarithmic  scale  near  A in  Figure  1-10,  can  be  signi- 
ficant in  rms  sensor  noise  effect.  Reduction  in  this  ’cost  of  feed- 
back' is  the  primary  motivation  for  turning  to  multiple-loop  design. 


1.3.1.  Reduction  in  Cost  of  Feedback  by  means  of  linear 

time-varying  compensation  and  nonlinear  compensation. 
To  reduce  the  hf  sensor  noise  effect,  one  way  is  by  linear 
time-varying  compensation  if  the  problem  has  time-varying  features 
[Rl].  Another  is  by  nonlinear  compensation.  Actually  the  so-called 
“adaptive1  system  is  in  the  category  of  nonlinear  compensation. 

They  may  or  may  not  be  better  than  fti  compensation  in  reducing 


u 

u 


the  'cost  of  feedback'.  It  is  noteworthy  and  scandalous  that  in 
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tho  vast  literature  on  adaptive  systems,  there  is  hardly  ever  any 
quantitative  comparison  between  the  adaptive  design  promoted  and  a 
proper  iti  design  accomplishing  the  same  design  objectives.  One 
could  excuse  this  not  being  done  in  a general  manner,  because 
there  is  hardly  any  'adaptive'  method  permitting  quantitative 
design  in  the  sense  here  defined.  However,  it  could  at  least  be 
done  experimentally.  Occasionally  one  sees  a comparison,  with  an 
'ordinary'  or  so-called  'classical'  design.  But  the  comparison 
is  usually  greatly  biased,  because  generally  some  very  naive  fti 
design  is  used,  and  thpre  is  no  statement  cf  specifications  - 
even  made  up  after  the  fact.  There  is  not  recalled  a single  com- 
parison, on  the  part  of  the  proponents  of  adaptive  systems,  with 
the  tti  quantiative  design  technique  [H3]  discussed  here.  Some 
nonlinear  compensation  techniques  for  which  a quantitative  design 
theory  exists  to  a greater  or  lesser  extent  have  appeared  in  the 
literature  [H5,  H9,  H10,  K2]  for  which  such  comparisons  are 
possible.  It  is  noteworthy  that  these  were  expressly  motivated  by 
the  desire  to  reduce  the  'cost  of  feedback',  so  that  such  comparisons 
were  a natural  by-product. 

1.3.2.  Multiple- loop  feedback. 

Another  method  of  'cost  of  feedback'  reduction,  in  the  con- 
text of  iti  design,  is  by  means  of  multiple-loop  feedback, 
restricted  to  those  cases  where  an  additional  plant  variable  (besides 
the  plant  output)  are  available  for  feedback  purposes.  Such  a 
multiple  loop  design  technique  was  first  developed  (H4J  for  the  cas- 
caded structure  of  Figure  1-12. 


-23- 


In  thia  thesis,  quantitative  synthesis  is  extended  as  follows: 

(a)  In  Chapter  2,  to  the  2-branch  parallel  structure  in  Figure  13a, 
in  which  one  branch  has  n cascaded  sections. 

(b)  In  Chapter  3,  to  the  natural  generalization  of  Figure  13a  for 
n-2  , into  the  structure  of  Figure  13b. 

(c)  In  Chapter  5,  to  the  fairly  general  structure  of  Figure  14  of  m 
parallel  branches  with  the  i-th  branch  consisting  of  a cascade 
of  n^  sections. 

1.3.3.  Plant  Modification. 

It  is  noted  that  in  Figures  1-12  to  1-14,  each  feedback  loop  is 
returned  to  the  plant  input  X . No  feedback  is  allowed  to  any  internal 
plant  variables  e.g.  from  C to  C2  or  more  generally  from  any 
to  Cj  , j > i and  j - 2,3,...  n-1  . The  reason  is  that  any  such 
internal  feedback  constitutes  "plant  modification".  The  plant  has  been 
assembled  by  its  specialists  to  deliver  some  maximum  output  C and  the 
permissible  levels  of  C 2 - C/Pj  , C3  - C/PjPj , . . . ^ - C/P1«..Pi_1 

are  thereby  determined.  Suppose  there  is  feedback  from  c vis 

r 

to  Cj  , as  shown  in  the  insert  in  Figure  1-12.  Now  X2  - — (1  + 
with  signal  level  possibly  much  greater  than  the  previous  C/P^  , which 
the  plant  may  perhaps  not  be  able  to  handle. 

We  thus  assume  that  the  "feedback  specialist"  is  called  in  to 
design  the  feedback  network  around  the  plant,  after  the  latter  has 
been  built  by  the  "plant  specialist".  This  is  the  situation  very 
often  in  practice.  If  the  feedback  specialist  does  his  job  properly, 
i.e.  achieves  the  system  response  function  T(s)  within  its  tolerances 
V P € ^ , then  the  signal  levels  inside  the  plant  will  be  within  the 
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values  allowed  by  the  plant  specialist,  so  long  as  the  command  input 
functions  r(t)  are  in  the  set  for  which  the  system  was  designed. 

Recently,  quantitative  design  has  been  extended  for  the  first 
time  [Wl] , specifically  to  the  single  branch  cascade  plant,  with  plant 
modification  allowed.  The  degree  of  modification,  in  a rms  sense,  of 
the  internal  plant  variables  was  added  as  one  more  design  specification 
in  addition  to  those  listed  here.  It  was  shown  that  the  loop  band- 
widths  can  thereby  be  significantly  decreased,  beyond  that  possible  in 
"no  plant  modification"  designs.  This  indicates  that  in  significant 
plant  uncertainty  problems,  it  is  definitely  advantageous  to  have  the 
"feedback  specialist"  participate  with  the  "plant  specialist"  in  the 
design  of  the  plant.  In  this  thesis,  as  previously  noted,  plant 
modification  is  not  allowed  because  this  is  very  often  the  situation 
in  practice.  Also,  it  is  necessary  to  pose  such  designs  in  order  to 
compare  the  results  with  plant  modification  designs. 
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The  multiple-loop  no  plant  modification  synthesis  theory 
developed  here,  is  based  to  a large  extent  on  the  design  theory  for 
the  cascaded  multiple-loop  design  theory  for  Figure  1-12.  Hence,  it 
is  very  important  to  first  thoroughly  understand  the  essentials  of  the 
above  cascaded  design  philosophy,  which  is  therefore  next  presented. 

1 .4  Review  of  Cascaded-Loop  Design  Theory 

1,4.1.  Cascaded  2-loop  design. 

The  basic  idea  is  to  use  the  inner  loop  Ln  « GjP., 

Figure  l-15(a),  to  minimise  the  effect  of  sensor  noise  at  the 


Figure  1-15  (a).  Cascaded  2 -loop  system  — no  P.M. 


plant  input  x . { is  the  square  root  of  the  power  spectrum  of 

the  sensor  noise  source.)  This  effect  is 


G1G2 


N1  1+G2r2+GlG2Pir2 


GlG2/(l+G2P2) 


P1P2G1G2 

1+G2P2 


1 + 


(1.4-1) 


J 

i 


(1.4-2) 


? 
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whore  P.. 

2o 


(:2P2 

HG2r2 


V(iW 

1 + L. 


* P1CIP2. 


(1.4-3) 


(1.4-5) 


1 

- p-  — in  the  hf  range  where  |l  | « l . (1.4-6) 

1 r 2 1 

x 

Hence  to  reduce  jj—  in  the  hf  ranye,  on<  must  try  to  reduce  |l  I . 

1 1 

But  L.  must  cope  with  the  uncertainty  i P p “ P,  (P-G_/l+P  G ) ) . 

Assuming  the  worst  case  of  uncorrelated  t\ certainties  in  P^  and 

P,  , the  best  that  can  be  done  by  the  in'.er  loop ' P,G,  is  to 

wipe  out  the  uncertainty  of  P,^  , s.  that  need  only  cope 

with  the  uncertainty  of  Fj  . Physic lly,  this  makes  sense  - for 

obviously  the  inner  loop  cannot  take  are  of  the  uncertainty  of 

P^  . There  is  then  left  the  single-1  >p  system  of  Figure  l-15(b) 

and  can  be  designed  to  handle  t » ignorance  of  P^  only.  The 

resulting  Lj  is  therefore  more  eco  incal  in  bandwidth  than  its 


\ 


Figure  l-lS(b).  Equiv  . ont  single*  loop  structure 
of  ascaded  2- lo>\i  system. 
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counterpart  in  the  single  loop  system,  for  the  realization  of  the 
same  C/R  specifications.  For  example,  compare  the  appropriate 
TN^  in  Figure  1-16  (a)  (logarithmic  scale)  and  in  Figure  1-16  (b) 
(arithmetic  scale)  of  the  numerical  example  in  Section  1.2-5. 


This  looks  very  good,  but  the  obvious  question  is:  What  of 
the  effect  of  sensor  noise  ^ not  present  in  the  single  loop 
design?  It  would  first  appear  that  the  inner  loop  G2P2  would 
have  to  be  enormous  in  magnitude  and  bandwidth,  leading  to 
tremendous  effect  of  N2  . But  this  is  not  so.  And  the  basic 
reason  is  available  if  one  studies  the  mechanics  of  sensitivity 
reduction  by  frequency  response  methods.  The  reason  is  that 
the  outer  loop,  even  if  optimally  designed  to  handle  a 
certain  definite  amount  of  uncertainty,  is  nevertheless  able 


] 

.1 

I 

I 
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Figure  1-16  (b) . Comparison  of  noise  response 

| TN  | on  arithmetical  scale  plot. 

to  handle  a "much  larger"  amount  of  uncertainty.  This  is  nicely 
seen  in  the  Nichols'  chart  in  Figure  1-17. 

Thus,  in  Figure  1-17,  suppose  the  uncertainty  in  P^Cjw)  is 
given  by  the  template  shown  of  F^  which  is  not  a point  (it  would 
be  a point  if  there  was  no  uncertainty)  but  a region.  Suppose  the 
specifications  require  the  closed  loop  response  uncertainty  to  be 


//'/''Template  of  P 

Y/'jZO////// 


The  range  of  P,  »gnoronc 
can  be  handled 


Figure  1-17.  The  template  P.  and  bound  on  L 


not  more  than  3 db 


and  the  optimum  design  locates  L 


in  the  position  shown,  i.e.  It1 


0.2  db  . Note  now  that  the  actual  uncertainty  in  P 

^ 1 

the  entire  shaded  region  lying  in  between  the  loci  of 
I T • | « .2  db  and  | T * | « -2.8  db.  « 


could  be 


Thus  P 


can  in  practice  have  significant  uncertaint 


even 
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1.4.2.  Nature  of  "free"  uncertainty. 

In  Figure  1-18,  suppose  arg  Lo  * -130°  , the  template  of  p is 

a vertical  line  20  db  in  length,  | AT | allowed  is  8 db,  and  | P | is 

chosen  |p|  . It  is  easily  seen  that  the  minimum  satisfactory 

|L  | “ -7.7  db  . But  it  is  also  easily  seen  that  this  minimum  L is 

satisfactory  for  any  uncertainty  (for  the  vertical  line  template) 

> 20  db  , even  infinite.  In  fact,  the  template  could  even  expand  to 

the  right  significantly,  as  shown  by  the  dashed  lines  in  Figure  1-18. 

This  kind  of  situation  is  true  for  all  -180°  < Arg  L < -90°  . It  is 

o 

not  true  for  Arg  Lq  € (0,-90°l  , but  for  significant  uncertainty  it 

is  almost  true.  Thus,  in  Figure  1-18,  suppose  Arg  Lq  « -50°  with 

the  same  plant  template  and  |at|  allowed  is  3 db,  giving 

|L  | * 4.8  db  . If  5.3  db  is  used  for  II  ( instead  of  4.8  db, 

o min  1 o 

then  this  Lq  is  satisfactory  (for  a vertical  line  template)  for  any, 
even  infinite  uncertainty. 

In  view  of  the  above  phenomenon  there  is  usually  little  advantage 
in  having  the  outer  loop  ignore  the  uncertainty  in  of 

Figure  1-15  (a)  in  to  ranges  in  which  is  in  the  above  regions  in 

the  Nichols  chart.  The  exception  is  when  the  uncertainty  in  is 

small  in  the  low  and  intermediate  (relative  to  the  overall  system 
bandwidth)  to  regions,  but  the  uncertainty  in  P2  is  large.  However, 
in  control  systems  the  power  levels,  and  with  them  the  plant  non- 
linearities  (equivalent  to  uncertainties,  see  (Hill)  and  uncertainties, 
usually  increase  as  one  proceeds  from  plant  input  to  output. 

The  situation  changes  radically  after  has  "turned  the 

corner",  say  at  point  A in  Figure  1-19.  At  point  B (with 
|PQ|  " lPljnin  ^ and  a vert*cal  template,  it  can  handle  uncertainty 
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of  BB'  only  and  not  any  bit  more.  Hence,  in  this  range  it  is  very 

advantageous  to  have  the  outer  loop  handle  only  the  uncertainty 

in  P^  in  Figure  1-15 (a) . Suppose  in  this  relatively  high-frequency 
)c  K 

region  P = , P = ~~  , k € [4,1250]  , 

1 ei  2 e2  1 

s s 

k2  € (10, 000j 312,900]  . In  a single-loop  design  the  "universal  high- 
frequency"  (UHF)  boundary  of  Section  1.2.3  is  H AVJ  KU  in  Figure  1-19, 
with  practical  optimum  Lg  shown.  The  length  of  B'BJ  is 
20  log1Q  9778.125  » 80  db.  In  the  2-loop  design,  need  handle  P^ 
uncertainty  of  50  db  only,  giving  the  UHF  boundary  HABKU  with 
B'B  = 50  db  in  length,  with  practical  optimum  outer  loop  in 

Figure  1-19.  In  practice,  it  is  best  to  compromise  a little  in  this 
range  and  "overdesign"  the  outer  loop  by  amount  db,  because  small 
overdesign  of  considerably  eases  the  burden  on  the  inner  loop 

L20  in  this  range  (see  Figure  1-22  for  the  trade-offs) . 

1.4.3.  Advantage  of  2-loop  over  single-loop  design. 

The  improvement  achieved  by  the  above  is  much  better  seen  in 

Bode  plots  of  L and  . First,  it  is  important  to  see  on  the 

Bode  plot,  the  nature  of  the  "practical  optimum"  Lgo  of  Figure  1-8. 

This  is  shown  in  Figure  1-20.  Lgo  is  moving  down  HAV  in  Figure  1-19 

(at  maximum  premissible  phase  lag  180-0  , in  order  to  decrease  as 

m 

fast  as  possible) . Its  average  slope  is  approximately 
-(180  - 0 ) 

jgo x 40  db  per  decade  [H4]  . If  the  shaping  of  Lgo  is  done 

by  hand  then  such  an  average  slope  can  be  achieved  by  means  of 

alternating  negative  real  poles  and  zeros,  i.e.  over  one  interval 

L (s)  has  an  excess  of  2 poles  over  zeros,  followed  by  one  with  an 
so 
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excess  of  only  one.  Near  the  corner  at  V in  Figure  1-19,  one 
arranges  for  an  excess  of  two,  then  introduces  a zero  to  provide  some 
extra  phase  lead,  followed  by  some  underdamped  complex  pole  pairs  as 
soon  as  possible  (see  [H7,  p.  384]).  The  greater  eL  , the  excess  of 
poles  over  zeros  of  L , the  further  off  these  complex-pole  pairs. 

b 

If  it  is  done  by  an  automatic  computer  program,  then  the  pole-zero 
pattern  which  emerges  depends,  of  course,  on  the  nature  of  the  program. 
However,  the  result  in  the  to  domain  obtained  by  several  programs 
tried,  is  very  similar  to  that  obtained  by  hand,  but  the  latter  is 
generally  more  economical  in  the  number  of  poles  and  zeros.  In  any 
case , in  this  work  we  use  eL  = 5 and  the  "universal  high-frequency 
cut-off  characteristic"  (UHFCOC)  shown  in  Figure  l-24a. 

The  following  is  very  important  in  revealing  the  basic  simplicity 
of  multiple- loop  design.  The  essential  difference  between  and 

L^0  is  that  UHFCOC  for  L^p  is  introduced  (fr-a^)  db  higher  up  than 
for  L , qivinq  BB'B"  instead  of  JJ'J"  in  Figure  1-20.  Hence,  if  a 

SO 

single-loop  design  giving  Lsq  has  been  found  to  solve  the  problem, 

it  is  then  very  easy  to  find  the  outer  loop  L^p  for  a two  (or 

higher) -loop  structure,  without  a detailed  design  of  L^p  . One  simply 

translates  the  UHFCOC  portion  of  Lgo  vertically  by  (4>-a^)  db  (and  to 

the  left  to  merge  with  L as  in  Figure  1-20) , where  $ is  the  high- 

so 

frequency  (hf)  uncertainty  of  P2P3>..Pn  , on  the  assumption  that  in 

this  to  region  P„P-...P  , is  closely  equal  to  its  high-frequency 

2 3 n 

eI+e2+‘ ’ "en 

asymptote  k_k  ...k  /s  — an  assumption  made  throughout  this 


work. 


Improvement  in  sensor  noise  effect 


It  was  shown  (Equation  1.3-1)  that  in  the  single-loop  design 

(X/N. ) = -L  /P  P in  the  hf  region,  whereas  (Equation  1.4-6)  in  the 

1 S s i z 

2-loop  structure  (X/N^^  = _L1/p1P2  • Ifc  easy  to  compare  the  two 

noise  responses  in  Figure  1-20.  Sketch  PioP20^u^  * clearly 

(X/N^)  < (X/N1>2  in  the  region  BB',  but  (X/N^)2  <<  (X/N1)g  in  much 

of  the  region  [oj  ,°°]  . In  any  case,  this  simple,  easy  construction 
K. 

permits  one  to  see  at  once  the  advantage  of  the  2-loop  design,  with 

respect  to  effect  of  sensor  noise  . Of  course,  a new  sensor  with 

2 

its  noise  power  spectrum  N2  has  been  introduced  and  its  effect  at  X 
must  be  considered.  For  this  purpose,  the  nature  of  the  inner  loop 
L_  is  next  reviewed. 


^AB' 


u>*  rM 
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\r  teix 


Figure  1-20.  Improvement  in  sensor  noise  effect 


Design  of  first  inner  loo 


The  outer  loop  has  been  designed  to  handle  all  uncertainty 

in  P^  and  P2  , up  to  the  corner  frequency  denoted  by  in 

Figure  1-19,  at  which  turns  the  corner  of  its  UHF  boundary. 

Hence,  for  u>  < the  demands  on  L^q  are  very  minor  (see 

Figure  l-9b  for  examples  of  the  bounds  B2(uj)  on  L2Q  * G2P20  in 

Figure  l-15a) . But  B (oi)  are  significant  for  w > id  , , where  L,  _ 

2 xl  10 

handles  the  uncertainty  only  in  P^  . To  see  this,  note  that  in 
Figure  1-21,  it  is  required  that  the  actual  (due  to  P and  P 

h 

uncertainties)  - ■■  ■ does  not  penetrate  into  the  locus  y.  , where 

S 1 So 

L1  = G1P1  1+lT  ' L2  = G2p2  ' With  the  nominal  L10  = G1P10  I+L  ' 

2 20 

L20  = G2P20  ‘ Due  to  P^  uncertainty  only  (with  L2  fixed  at 
L20  ^ ' L1  exten<^s  up  to  the  dashed  part,  e.g.  at  up  to  point  E 

s 

in  Figure  1-21.  It  is  therefore  necessary  that  t— — does  not  cause 


Figure  1-21.  Explanation  of  finding  an  acceptable  region  for 
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further  extension  into  , e.g.  EF^  is  not  allowed,  but  EF2  is 
allowed.  In  this  way  an  acceptable  region  for  ■— £—  at  each  w , is 

L2 

found.  What  are  the  bounds  on  Lori  such  that  — ; — stays  in  its 
acceptable  region  for  all  P2  ? This  is  almost  precisely  a repeat  of 
the  single-loop  problem  and  is  solved  in  the  same  manner.  It  is  also 
necessary  to  introduce  y2  (analogous  to  y^  for  the  single-loop 
problem)  and  with  it  a UHF  boundary  for  L2q  . The  reason  is  that  in 
Figure  1-1 5a, 


1 

(1+Ll) (1+L2) 


(1.4-7) 


Hence,  for  the  same  reason  as  in  Section  1.2.3,  it  is  necessary  to 

L„ 


assign  bounds  to 
Bh2  ' 


1+L„ 


$ y2  , resulting  in  the  L2Q 


UHF  boundary. 


Typical  bounds  B2(w)  on  L20  , which  result  from  the  above, 
are  shown  in  Figure  l-9b.  The  reason  for  their  nature  is  discussed 
in  [ H4  ] . The  resulting  "practical  optimum"  L2Q  (in  the  sense 
previously  discussed  in  1.4-1  for  L^g  ) is  also  shown  in  Figure  l-9b. 
It  is  useful  to  note  that  at  some  frequency  w = u>m2  (note  its 
location  with  respect  to  Bhl  in  Figure  l-9a),  |L2Q(jiu)|  has  a peak 
value  M2  at  which  Arg  L20^tllm2^  ~ “90°  • The  value  of  M2  is  a 
function  of  0^  , the  phase-margin  of  (related  to  y^  of 

Figure  l-9a)  and  of  the  overdesign  of  . Empirically 

obtained  graphs  of  M2  vs.  these  two  variables  are  shown  in  Figure 
1-22,  and  they  will  be  very  useful  later.  Note  also  the  existence  of 
the  "corner"  frequency  wx2  in  Figure  l-9b  (analogous  to  u>xI  of 
Figure  l-9a) , at  which  L2Q  turns  the  corner  at  the  bottom  of  Bh2  . 
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Figure  1-22.  Peak  value  of  inner  loop  as  function  of  0 


Note  that  for  id  > idx2  , the  shape  of  the  practical  optimum 
1^0  (j“>)  is  very  similar  to  that  of  ( jcu)  for  <d  > <dx1  . In  these 

respective  ranges  the  ( joj)  assume  their  respective  UHFCOC  (recall 
Section  1.4.3),  which  are  similar  functions  of  the  phase  margin  0 

mi 

and  of  eLi  the  excess  of  poles  over  zeros  of  . 


Sensor  noise  effect  in  first  inner  loop. 


In  Figure  l-15a,  it  is  readily  found  that 


X L2^P2  • L2 

" (l+I^)  (1+I<2 ) " 


(1.4-8) 


in  the  hf  region  where  |l*  |,|l*2|  <<  1 • Hence,  on  a Bode  plot 
containing  |l*20l  ' one  onlV  needs  to  sketch  |p20|  ' to  see  I X/N2 1 . 
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On  the  same  Bode  plot  let  L , L, . , L_rt  , PinPon  , P__.  be  sketched. 

SO  1U  ZU  1U  ZU 

One  then  sees  at  a glance  the  advantage,  if  any,  of  a 2-loop  design 
over  the  single-loop.  If  the  major  effect  is  in  the  hf  region,  the 
effect  of  in  a single-loop  design  is  the  difference  between  Lsq 

and  P'l0P20  • In  a two-loop  design  it  is  the  difference  between 

and  p.gP20  for  N1  ' and  the  dlfference  between  L2Q  and  P2Q  for 

N-  (see  Figure  1-23). 


1.4.5.  Design  of  second  and  higher  inner  loops. 

/ 

Suppose  n - 3 in  Figure  1-12,  so  a second  inner  loop  = G3P3 
is  feasible.  If  is  not  used,  then  from  1.4.4,  the  first  inner 

loop  handles  the  hf  (u>  > wxl)  uncertainty  of  P2P'3  • L3 

used  then  the  relations  between  L3  and  L2  are  very  similar  to 
those  between  L2  and  , respectively:  L2  is  allowed  to  handle 

the  uncertainty  in  P2P3  UP  to  wx2  L20  ^ ’ This  a,X2  is 

obtained  in  the  same  manner  as  u)^  of  L^0  , i.e.  B^2  , the  UHF 
boundary  for  L2Q  is  obtained  by  assigning  to  L2Q  only  the  hf  un- 
certainty of  P2  (not  of  P2P3  ) and  some  extra  margin  a2  (analogous 
to  of  L10  ).  By  definition,  I<20  turns  the  comer  of  Bh2  at 

ux2  (point  in  Figure  l-9b) . 

The  bounds  on  are  very  minor  for  w < u>x2  (just  as  those 

°n  L20  were  ^or  u < uxi  )•  They  become  significant  for  u>  > u>x2 
(see  Figure  l-9c)  and  the  resulting  |l3q|  reaches  a maximum  at 
Ar9  1<30  ~ -90°  at  u ■ , with  L20  located  relative  to  Bh2 

exactly  as  Lio^(i>m2^  Was  re^atlve  to  (compare  points  u>ro2  , 

wm3  in  Fi9ures  l-9a,  9b).  The  value  of  = lL30^wm3^  is  the  same 

function  of  0^  and  o2  as  M2  - I L20 ( j“m2)  I is  of  % and  al  ' 
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so  the  graphs  in  Figure  1-22  may  be  used.  There  is  needed  a Bh3  , 
the  UHF  boundary  for  L^q  , exactly  like  B^2  for  L^q  • The  shape 
of  Ljq  is  therefore  very  similar  to  that  of  L20  , except  that  it 
unfolds  in  a higher  frequency  range.  This  is  called  the  bandwidth 
propagation  effect.  The  effect  of  sensor  noise  N3  is  found  to  be 

_X_  L3/P3 . L3  . 

n3  * (1+1^)  (1+l2)  (i+l3)  = p3 

in  hf  range  where  each  ]l^|  <<  1 . Thus  one  needs  only  sketch 
|l30|  , |P30!  on  a Bode  plot  to  see  the  order  of  magnitude  of 
| X/M3 | (see  Figure  1-23). 

1.4.6.  Bandwidth  propagation  effect  in  cascade  design. 

It  was  shown  in  Section  1.4.3  how  the  UHFCOC  of  may  be 

obtained  from  that  of  L . The  peak  point  of  (l_.J  at  hi  . lies 

so  20 1 m2 

near  the  end  of  the  almost  horizontal  part  of  the  UHFCOC  of  L^0  . 

It  is  then  necessary  for  | l2q  ( jn>)  | to  decrease  until  it  reaches  the 
u>x2  point  of  Bh2  at  which  its  UHFCOC  takes  place.  Near  the  end  of 
the  almost  horizontal  part  of  the  latter,  there  is  the  value  at 

which  t.30(jw)  has  its  peak.  L3Q(ja >)  must  then  decrease  till  it 

reaches  the  u>x3  point  of  B^3  , and  then  its  UHFCOC  takes  place,  etc. 
At  each  step  one  is  pushed  to  a higher  frequency  range.  But  it  is 
interesting  to  see  in  Figure  1-23  that  the  UHFCOC  of  I.  and  of  the 
last  inner  loop  LnQ  = in  Figure  1-23,  are  not  far  apart  in 

frequency,  with  the  latter  somewhat  in  a higher  range,  because  there 
is  inevitably  some  extra  frequency  range  consumed  at  each  inner  loop 


stage. 
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1.5  A Simple  Fast  Technique  for  Cascade  Design  Perspective 
The  cascade  synthesis  procedure  has  a highly  useful  property. 

One  can  achieve  excellent  design  perspective  by  means  of  simple 
straightforward  sketches.  Only  an  initial  single-loop  design  need  be 
done  in  some  detail.  These  preliminary  perspective  sketches  enable 
the  designer  to  decide  which,  if  any,  and  how  the  available  plant 
internal  variables,  should  be  used  for  feedback  purposes.  He  can  then 
perform  his  detailed  multiple-loop  design  with  considerable  confidence 
that  iteration  will  not  be  necessary.  Most  of  the  design  parameters 
needed  for  'design  perspective'  have  been  defined: 

0 . , the  phase  margin  for  L (Figures  l-9a,b,c)  related  to  y.  , 

mi  io 

the  maximum  permitted  value  of  |Lj/(1+L^)(  . 
ok  , the  overdesign  of  (Figures  l-9a,b,c). 

u>  . , the  frequency  at  which  L . turns  the  "corner"  at  the  bottom  of 
xi  io 

Bhi  extended  downward  by  the  additional  cr  value. 

* the  frequency  at  which  |L^(ja))|  has  its  maximum  value,  and 

Mi  = lLio(ju,mi)  I ' t*ie  co^resP°n^in9  peak  value  of  lL^Ql  for 
i — 2, . . . ,n  . 

UHFCOC  (L^q) , the  universal  high-frequency  cut  off  characteristic 
effective  for  a)  > for  the  outer  loop.  That  shown  in  Figure  l-24a 

is  used  here,  obtained  for  e,  » 5 , 0 , ■ 50°  . Corresponding 

ml 

characteristics  for  the  inner  loops  are  shown  in  Figure  l-24b,  labelled 

L.  (HF)  . There,  they  are  given  for  different  0 values,  but 
io  m 

eL  * 5 for  all. 

1 A 

In  addition,  we  need:  A.  = Min  (B.  .)  - a.  , where  Min  (B.  ,) 

i ni  i ni 
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I 


is  the  point  of  minimum  magnitude  of  Bhi  (=  -54.5  db  in  Figure  l-9a 
for  i - 1 ) . 

Also,  we  need  the  value  of  u>  . and  the  shaoe  of  L.  for 

mi  10 

u)  6 • The  latter  was  obtained  by  studying  numerous  designs. 

It  is  given  in  Figure  l-24a,  labelled  L^o(IF)  , since  it  applies  in 
the  intermediate  frequency  (IF)  region.  The  arrows  in  Figures  l-24a,b 
are  used  to  locate  , as  explained  below. 

1.5.1.  Procedure  for  "design  perspective". 

1.  Make  a single-loop  design  Lgo  to  handle  the  entire  problem, 

e.g.  L in  Figure  1-23  for  the  example  of  Section  1.2.5.  It  is 
so 

assumed  that  the  hf  portion  of  Lgo  is  reasonably  close  to  that  shown 
i.e.  UHFCOC(L1Q)  in  Figure  l-24a. 

2.  Let  4>-  be  the  hf  uncertainty  of  P_P_...P  , and  allow  a. 

£ £ 3 n 1 

an  overdesign  margin  for  L1Q  . obtain  L1Q  by  shifting  the  UHFCOC 
portion  upwards  by  ^2"’al>  ^ aS  shown  -*-n  Figures  1-20,  1-23. 

3.  The  value  of  w is  at  the  arrow  on  UHFCOC (L,„)  in  Figure 

mz  10 

l-24a.  The  value  of  M2  = lL2o^a,m2^  is  available  from  the  graphs 
in  Figure  1-22.  Locate  M2  at  w = a)m2  an<*  this  be  point  Q in 
i.Jre  1-23.  Point  Q is  | L2Q ( ) | . Draw  a horizontal  line  in 
Figure  1-23  at  A2  magnitude  (=  -25  db  for  above  example) . Trans- 
parencies of  Figures  l-24a,b  are  assumed  available.  Place  L.  (IF) 

10 

of  Figure  l-24a  on  Figure  1-23  so  that  the  two  Q points  coincide.  Find 
where  L^q(IF)  of  Figure  l-24a  intersects  the  Aj  line  of  Figure  1-23 
and  label  this  point  of  intersection  as  point  C in  Figure  1-23. 

4.  Pick  the  L^o (HF)  curve  in  Figure  l-24b,  according  to  the 


I 
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1 


0m2  value  being  used  for  L2Q  . Lay  (HF)  on  Figure  1-23,  such 

that  the  two  C points  (of  L.  (HF)  and  that  obtained  in  Step  3) 

10 

coincide.  L_.  consists  of  L.  (IF)  of  Figure  l-24a  in  the  inters 

10 

mediate  u>  range  and  of  L^q(HF)  of  Figure  l-24b  in  the  high- 
frequency  range.  Use  the  portion  of  L^q (HF)  to  the  left  of  C to 
obtain  a smooth  continuous  curve  for  L2Q  . 

5.  Steps  3,  4 are  repeated  in  order  to  determine  L3Q  . use 

the  arrow  on  L.  (HF)  of  Figure  l-24b,  to  locate  cj  , . Use  Figure 
io  mi 

1-22  to  obtain  M3  , giving  a new  point  Q in  Figure  1-23.  Then  lay 
Lio(IF)  of  Figure  l-24a  on  Figure  1-23,  so  that  the  Q's  coincide,  etc. 
A horizontal  line  of  value  £>3  is  drawn,  etc.  The  entire  process  is 
repeated  until  all  the  loops  are  exhausted. 

The  results  in  Figures  1-23,25  are  examples  of  the  excellent 
agreement  obtained  in  practice  between  these  quick  estimates  and  the 
actual  final  detailed  designs. 

6.  After  each  L.  is  obtained,  it  is  a good  idea  to  sketch 

io 

the  effective  P values  to  use  for  the  sensor  noise  effect.  Thus 

after  L is  obtained,  sketch  |p  I =1 P. rtP_rt. . .P  I in  Figure  1-23. 

If  there  is  little  sensor  noise  amplification  ( lLS0/^0l  not  large 

over  a large  u>  range) , there  may  be  no  point  in  using  more  feedback 

loops.  After  L^q  has  been  obtained,  it  is  easy  to  see  the  saving 

in  sensor  N,  noise  effect,  by  using  L_  . Sketch  |?__ I «| P_rt. . .P 

x i 20  20  no1 

to  see  the  hf  N2  sensor  noise  effect  ( «|  ^20^2o\  * * similarly 
lLi</^iol  ' I Piol  =l  Pio"  * *PmJ  ' ^ives  the  hf  IT  sensor  noise  effect. 

The  designer  has  to  decide  which  sensor  points  to  use  and  the 
corresponding  trade-off  values.  It  is  emphasized  that  this  design 


J 


i 

I 

I 
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Figuro  l-25c.  Demonstration  of  design  perspective. 


perspective  approach  does  not  give  him  this  information  immediately. 
Design  perspective  is,  however,  a very  useful,  fast,  simple,  tool  for 
making  his  choice.  It  makes  it  easy  for  him  to  try  different 
combinations.  For  example,  he  may  find  it  questionable  to  use  a sensor 
at  C2  perhaps  because  N2  is  very  large.  He  can  then  consider  using 
a larger  value  of  , thereby  throwing  more  burden  on  L^0  and  less 
on  L2q  . He  need  not  execute  a detailed  design.  The  above  design 
perspective  gives  him  the  answer  in  a few  minutes.  He  may  decide  to 
forego  the  sensor  at  C2  . If  so,  in  our  terminology  the  new  - the 
old  P^P2  ' *■•••  t*'e  °utcr  1°°P  must  now  handle  V <d  the  uncertainty 
in  the  former  P^^  . 

Our  objective  is  to  develop  similar  design  perspectives  for  the 
more  complex  structures  of  Figures  13a,  13b,  14. 
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Appendix  1 

Transfer  Functions  of  the  Numerical  Example  in  1.2.5. 

Loop  transmissions  for  satisfying  the  appropriate  boundaries  were 
obtained  sometimes  by  means  of  a computer  program,  and  at  other  times 
by  hand.  The  former  are  characterized  by  much  higher  order  expressions. 
For  the  numerical  example  of  Section  1.2.5  (Figures  l-9a,b,c,  1-23), 
they  are: 


5(1  + 1,05s)  (1  + 1.087s) (1  + 0,05764s) (1  + ,03845s) 

10  " s(l  ♦ 1.954s) (1  + 1.89s) (1  + .164s) (1  + .0265s) (1  + .026s) 


(1  + .0074s)  (1  + .0028s)  [1  (2^0)2] 


<i..0O04s)2u*.  005,)  u*355§*  l^)2l  11*35555*  <2i5oo> 2)  2 


L20  " 


0.5 


1 + 


s 


13800 


[1  0.4s  . > 2i  [i  + °-,6s  + ( — - — 

1 1920  l1920'  1 1 16000  l16000 


2 2 
) ] 


So  ■ 


900 


0.645 


19200 


1 + 


45000 


(1  + -2.:AS  x ( 

1 11250  1 


11250 
1 


•>2]  (7  8 


45000  '45000 


>2) 


(1  + 


2s 


+ 


180000  '180000 


r>2] 


J 


L 
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C HAP TER  2 

THE  CASCADE  STRUCTURE  IN  PARALLEL  WITH  ONE  BRANCH 

2.1  Introduction 

Quantitative  feedback  design  theory  was  defined  in  Section  1.1. 
Prior  to  this  work  it  was  confined  to  the  cascaded  plant  structure  of 
Figure  1-12.  This  chapter  extends  quantitative  synthesis  to  the  plant 
structure  (heavy  lines)  of  Figure  2-1,  consisting  of  two  parallel 
branches,  one  of  which  has  n internal  sensing  points  C, ,C,,...C  . 

It  is  assumed  the  system  input  r and  output  c can  also  be  sensed, 
giving  n+2  degrees  of  freedom,  in  the  form  of  n+2  independent  data 
processing  of  these  measurements.  Figure  2-1  shows  one  canonic 
structure  which  exploits  this  freedom,  but  there  is  an  infinitude  of 
others  — see  Section  2.4. 

It  is  assumed  each  plant  section  has  independently  uncertain 
parameters  giving  for  each,  a set  {P^  £ . , and  the  objective,  as  in 


-I 


Figure  2-1.  System  structure  (Heavier  lines  indicate  plant  with  sensor 
noise  source's)  . 
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Section  1.1,  is  to  guarantee  that  these  objectives  are  satisfied 
V €^)  i . Just  as  in  the  cascade  case,  no  analytic  theory  is  avail- 
able but  an  engineering  design  theory  is  developed  for  a large  class 
of  practical  plants.  The  enphasis  is  on  the  minimization  of  net  sensor 
noise  effect  at  the  plant  input  X . A transparent  design  theory  is 
developed  which  is  very  similar  in  its  results,  to  that  of  the  cascade 
system,  although  its  development  is  more  difficult.  There  is  the  same 
'bandwidth  propagation'  phenomenon  as  in  Section  1.4.6.  Also,  "design 
perspective"  (Section  2.5)  is  possible  by  relatively  simple,  fast 
sketches,  based  on  only  a preliminary  single-loop  design  for  the 
problem.  This  enables  the  designer  to  decide  which  loops  to  use  and 
how,  and  gives  him  a very  good  idea  of  the  final  design  — without  a 
detailed  design.  Such  rapid  perspective  design  (Section  2.5)  is 
compared  with  detailed  design  for  a specific  3-loop  numerical  example 
with  large  uncertainty.  The  agreement  is  excellent  (Figure  2-5) . The 
treatment  is  devoted  entirely  to  minimum-phase  systems. 


2.2  Design  of  the  Outer  Loop 
It  is  recalled  that  the  design  objectives  are: 
T(jw)  ■ C(ju)/R(ju>)  satisfy  specified  bounds  of  the 

B2(u)  * |T(Ju»)  | < B1(w) 


and 


-L, 

l+L 


( Y a constant,  V ui 


with  » D2  ' Y apriori  given  and  L is  the  outer 
(It  is  more  convenient  in  this  chapter  to  use  L for 
for  the  first  inner,  etc.,  unlike  the  notation  in 


that 

form 

(2.2-1) 

(2.2-2) 

loop  transmission, 
the  outer  loop. 
Chapter  1 . ) Just 
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aa  In  the  cascade  case,  (2.2-1,  2)  are  achievable  [H3]  with  a single 


loop  G.  ■ 0 , 1 - l,...,n  , G-G_j<0  in  Figure  2-1.  But  the  result- 


ing Lg  • GgP  nay  then  require  very  large  bandwidth,  causing  great 


anplification  of  the  sensor  noise  N , as  in  a later  example  — Figure 


2-6a,b.  The  sinplistic  approach,  later  justified,  is  to  design  the 


outer  loop  L from  C to  cope  only  with  the  uncertainty  in  , Pc 


which  can  give  an  L much  more  economical  than  in  a single-loop  design 


in  which  L must  cope  with  all  P . The  first  inner  loop  from 

may  be  designed  to  cope  only  with  Pj^  , with  possible  great  saving 
compared  to  an  L^  which  copes  with  piP2"‘Pn  * Similarly  the  second 


inner  loop  need  cope  only  with  P,  , etc.  The  result  has  considerable 


transparency  and  insight,  giving  the  designer  a good  perspective  of  the 


optimum  division  of  the  feedback  burden  among  the  loops,  even  before 


execution  of  the  design  details  (Section  2.5).  Sinplifications 


initially  made  in  order  to  concentrate  on  the  essentials,  are  covered 


in  Section  2.6.  In  Figure  2-1,  let 


"a  S 'nVr"Vl  ' P ' Vb  * Pc  ■ 


2>  - (1  + P Gn  + P P .G_  . + ...  + P ...P.G.)  + PG 
n n n n-l  n-i  n 11 


- ^ + PG  - S&^l+L)  , L - gp  - PgG 

nf.i  M fGP/J). 


(2.2-3a-f) 


T(S)  - SisL  _ F®  _ _ F _L_ 

R(a)  £ l+(GP/j£)1)  1+L  ’ 


x G 

N £ 


_s y* 

(1+L)  1+L 


In  (2.2-3f)  the  sensor  noise  effect  is  examined  at  the  plant  in- 


put X where  it  tends  to  be  large  (Figure  2-6) , causing  plant 


saturation.  In  the  high  frequency  range  (denoted  as  hf ) , (2.2-3f)  + L/P 


^a...  a, i.K-vt: 
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where  |l(  jut)  | <<  1 but  |l/p|  can  b « very  large  (Figure  2-5)  — the 

9 

hf  range  is  the  major  trouble  source.  For  example,  in  Figure  2-6b,  the 
lowest  <i>  range  with  large  and  sharp  peaking  of  | X/N  | is  300  rps  at 
which,  from  Figure  2-5,  the  nominal  |l|  ~ -48  db  but  |p|  ~ -83  db  . 
Hence,  the  major  effort  in  sensor  noise  effect  minimization  will  be 
made  in  hf.  This  highly  justified  fact  is  very  helpful  in  simplifying 
the  problem.  Since  P is  contrained  in  (2.2-3b) , such  minimization 
requires  |l|  minimization.  But  from  (2.2-3d),  L must  cope  with 
pe  “ P/jOj  uncertainty.  Therefore,  for  maximum  economy  of  L , choose 
the  G^  , i ■ l,...,n  in  of  (2.2-3c)  to  minimize  the  uncertainty 

in  Pe  - P/^1  . 

Consider  accordingly  the  uncertainty  in 

P1P2...PnPb  ♦ Pc  e 

Pe  “ r+  Pncn  Pn---P2P1G1  in  hf<  Where  each  Pi  * ki/s  1 * ei 

the  excess  of  poles  over  zeros  of  ?i  . since  P^.-.P^  parallels 

Pc  in  Fi9ure  2-1,  it  is  assumed  in  the  meantime  (see  Section  2.6)  that 

(*1  * *2  + *n]  *b  " ea  + eb  “ ec  * Henc*»  »t  hf 


P - 


V, 


< k. +k 
a d c 


s [1+k  g +**+k  g .1 
n n a l 


k 

e_ 

•c 

s 


where 


(2.2-4a,b) 


_ A n 

G » g s 

n n 


n-1 


Vi8 


e +e  , 
n n-1 


G1  “ V * 


The  notation  in  (2.2-4b)  is  used  to  simplify  the  expressions.  The 
resulting  Gi  do  not  necessarily  have  zeros  at  the  origin  - see 
Section  2.3. 

The  range  of  ki  is  taken  as  la^bj  # b.^  > a.^  > 0 . In  the 
logarithmic  complex  plant  (Nichols  chart) , Pe  is  not  a point  but  a 
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set  (P  } because  of  the  uncertainty.  For  any  fixed  k_,***,k 
® In 

values  the  set  {p  } , due  to  (k.  ) , (k  ) in  (2.2-4a) , is  a vertical 

e dc 

r k Mb  -I 

line  whose  length,  Lgth  {P  ) - i ~ — is  a function  of  k , 

e l k a.  ta  l „ a 

1 a b c J db 

(see  Figure  2-2a)  and  is  maximum  at  k - a if  b /a  * b.  /a.  (at  b if 

a a c c b d a 

b^Jb^/s^  ).  The  former  is  assumed  because  P^  is  in  parallel  with 

P^P^  — see  Section  2.6.  Hence,  due  to  j*1_1_  the  k^  uncertainty 

sets  Lgth(P  } * (a  b.  +b  )/(a  a.  ♦ a ) , with  equality  iff  3 g.  such 
e a d c a b c i 

that  sets  {kgdi^.k^k^  } of  (2.2-4a)  5 as  the  ki 

independently  range  over  la^b^  . For  example,  gj^  - 1/2  fc^/b^,  has 
the  effect  shown  in  Figure  2-2b,  so  that  the  range  of  k is  AE  > AB  . 


It  is  readily  seen  that  such  gA  exist,  e.g.,  g2 


9n  ” 0 ' 


bfa/bc  ( « ab/,ac  ' cowpstible  with  the  previous  bc/a^  * b^/a^  . 

(In  the  case  bc/ac  ( b^/a^  tlie  ana^-°9ous»  compatible  condition  is 
b^/b^,  * g^  * a^/a^,  Tbus  at  hf  the  best  the  inner  loops  can  do  for 
the  outer  loop  L , leads  to  it  coping  with  a gain  uncertainty  set. 


{8^  + ^}  -]/'{PM  , of  Lgth]/* 


a b +b 
a b c 

a a +a 
a d c 


(2.2-5a,b) 


For  example,  if  n«2  , all  ai - 1 , bi»40,  10,  60,  200  for 

i»l,2,b,c  then  in  a single-loop  design  L must  handle  hf  { p > of 

2) 

length  l <Vb  + bc,/( Vb  + Vdb  " 81  • 7 wh«reas  (2. 2- 5b)  gives 
42,3  db,  a saving  of  39.4  db. 

The  hf  region  is  most  important  for  sensor  noise,  and  the  hf  form 
of  P^  in  (2.2-4)  greatly  simplifies  the  problem  there.  In  Chapter  1, 
Section  1.4.4,  it  was  shown  that  the  saving  in  the  outer  loop  is  very 
minor  for  u>  < u>x  (the  frequency  at  which  L turns  the  corner  - 
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Figure  l-9a).  Therefore,  just  as  in  the  cascade  case,  let  L handle 

the  entire  plant  uncertainty  problem  up  to  . The  value  of  w is 

obtained  by  having  L handle  the  smaller  hf  uncertainty  of  (2.2-5b). 

The  region  of  concern  for  the  inner  loops  is  thus  at  hf,  where  each 
e 

Pj  ■ kj/s  . Therefore,  to  simplify  the  presentation  and  notation  we 
e-« 

write  Pj  ■ kj/s  , j ■ a,b,c  as  if  (2. 2-4, 5)  apply  for  all  u . 
Outer-loop  design  is  now  a single-loop  problem  with  the  equivalent 
plant  P#  of  (2.2-4a)  denoted  by 


. P „P  +P 

o A aO  b c 

e ™ 1+P  nG  +•••+? 

nO  n nO 


i k. +k 
a d c 


*P20P10G1 

A k<> 
A e 


(2.2-6a,b) 


(1+a  g +***+a  g )s  c s 
n n a l 


o o 

and  L ■ P G . 

e 


The  super-oh  on  Pg  , k^  , L indicates  P , kg  , L with  P^ » P2 ’ * * * ' Pn 

at  their  nominal  values.  The  problem  is  to  find  G or  equivalently  a 

o o . 

nominal  (the  sub-oh  indicates  nominal  values)  Lq  * peoG  30  that 

A_  (u>) 

(2.2-2)  and  At  |T(ju) I S t . , . of  (2.2-1),  are  satisfied.  The 
n n Aj  (to) 

optimum  design  for  this  single-loop  problem  [H3]  was  reviewed  in 
Section  1.2,  so  we  proceed  directly  with  an  example  which  is  very  help- 
ful in  explaining  the  multiple- loop  design  theory. 
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2.2.1.  Design  example. 

In  Figure  2-1  let  n-2  , P , - k ,/s  , j«l,2,b  j P -k  /a3  i 

11  c c 

* 20  , »2  - 50  , ^ « 1 , ac  - 1000  , ^ - 800  , b2  - 500  , ^ - 60  , 

bc- 200,000  . Time-domain  bounds  on  the  acceptable  step  response  are 

shown  in  Figure  2-3a,  and  their  translation  into  bounds  on  |T(jw) | 

are  shown  in  Figure  2- 3b.  Such  translation  is  somewhat  of  an 

engineering  art,  but  in  practice  good  results  have  been  obtained  with 

moderate  effort  [H3,  Kl,  si].  Also,  y®2.3  db  in  (2.2-2),  corresponding 

to  "-23.5%  overshoot  for  the  second  order  model.  The  nominal  plant 

values  are  taken  as  a^^  (with  no  loss  in  generality).  Note:  The 

specifications  must  be  consistent  with  physical  reality,  i.e.,  it  is 

crucial  (H8]  that  3 u)q  , such  that  for  u>  > Wq  » the  largest  variation 
A2  (w) 

Pe  < *n  A~(io)  ' *n  or<*er  that  no  sensitivity  reduction  be  needed  at 

large  enough  u>  , permitting  L(ju>)  -*  0 as  to  -*  <»  . From  Figure  2-3, 
o>0  »»  55  for  a multiple- loop  design,  280  for  a single-loop  design, 
revealing  the  advantage  of  the  former  with  respect  to  the  outer  loop. 

Single- loop  design  theory  was  reviewed  in  Section  1.2.  The 
bounds  and  the  resulting  loop  transmission  on  L°(ja))  are  shown  in 
Figure  2-4  and  its  Bode  plot  in  Figure  2-5. 

It  is  easy  to  derive  the  single-loop  Lq  ( j<*> ) which  would  be  necessary 
for  the  same  problem  and  this  is  included  in  Figure  2-5,  to  emphasise 
the  saving  in  bandwidth  and  sensor  noise  effect. 

Returning  to  the  multiple  loop  design,  note  that  the  nominal 
L of  (2. 2- 3d),  denoted  by  LQ  is  identical  to  L°  and  the  nominal 
peo  of  (2'2_3d)  to  p°0  * justifying  G - LQ/Pe0  - L°/P°Q  . Here, 
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Lq(s) 


(1.6) IQ10  (s+1 .96)  (s+15.8)  <8+39. 5)  (s2+4 ,6s+10.7)  (s2+7.8s+22.2) 

2 2 2 
s (s+2 .16) (s  +4.7s+7.86) (s  +5.6s+25.1) (s  +25.2s+371)\Ms) 


<»(s)  * (s2+324s+153,800)  (s2+324s+158,700)  . 


(2.2-7) 


The  trade-off  between  complexity  of  the  functions  and  the  economy  in 
bandwidth,  is  again  emphasized.  Simpler  Lq(s)  may  conceivably  be 
as  good  or  ever,  better  than  the  above  in  satisfying  the  B(w)  , as 
no  expertise  in  this  art  is  being  claimed  here.  But  for  a given  skill 
in  the  art,  the  larger  the  order  of  LQ(s)  , the  more  economical  it  is 
in  bandwidth,  up  to  the  point  where  the  trade-off  becomes  very  small. 
The  designer  must  use  his  own  engineering  judgement  where  to  stop. 

Use  of  an  inner  loop  permitted  a maximum  reduction  of  39.4  db 


b2^5lbb  + bc 

(20  log  r—  ■ 20  log  93.7)  here,  giving  a saving  in  bandwidth 

*2*1^  + bc 


of  ~ 40/29  decades,  because  along  B,  in  Figure  2-4, 

U 

[Arg  Lq ( jk>)  ] 1*130°  , with  corresponding  average  d|L0(jw)  |/dw  = 

- (40)  » -29  db/decade  [H7] . Thus,  a>xs  (single-loop)  i 1900  , 

while  u ■ 80  , a difference  of  1.38  decades.  The  reduction  in  sensor 
x 

noise  effect  at  X is  enormous  (Figures  2-6a,b),  because  the  rms  noise 
value  is  obtained  by  integrating  arithmetic  values  on  an  arithmetic 
u scale. 


For  later  use,  the  following  o>  intervals  in  Figures  2-4,  5 


are  emphasized: 


- f0, u)  ) « [0,90)  , I_  - [u  ,w  ) » (90,330) 

/i  X d x z 


*C  " w [330, «)  . 


(2.2-8) 
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Figure  2- 6a.  Sensor  noise  effects  at  X — in  decibels 
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Figure  2- 6b.  Sensor 

noise 

effects  at 

X — arithmetic 

scale. 

The  design  of  the  first  inner  loop  is  decisively  influenced  by 

these  intervals  of  the  outer  loop  LQ  , and  it  is  very  important  to 
understand  their  role. 

2. 3 Design  of  Inner  Loops 
2.3.1.  First  inner  loop  . 

In  Section  2 the  inner  loops  were  apparently  sacrificed,  in  order 

to  obtain  the  most  economical  outer  loop  and  thereby  minimize  the 

e 

effect  of  sensor  noise  N at  X . =0  , G = s b /b  were 

2 3 r.  1 be 

found  satisfactory  for  this  purpose.  The  obvious  criticism  is  that 
this  G^  , besides  being  impractical,  would  tremendously  amplify  h£ 
noise  effect  at  X (Figure  2-1)  and  likely  more  than  cancel  the 
benefit  gained  for  L . The  answer  is  that  while  these  G^  are 
satisfactory,  there  are  other  much  smaller  acceptable  values.  This  is 
due  to  the  mechanics  of  sensitivity  reduction  such  that  L^ju^) 
optimally  designed  to  handle  an  uncertainty  set  can  in  practice 

handle  a set  much  larger  than  yf^)  (e.g.,  Figure  1-17) . 

Tliis  phenomenon,  of  crucial  importance  in  multiple-loop  design  using 
frequency  response,  was  emphasized  in  Section  1.4.1  in  connection  with 
cascade  design.  So,  the  next  step  is  to  find  the  bounds  B^w)  on 
the  first  inner  nominal  loop  Llf)(ju>)  such  that  the  economical  LQ 
of  Section  2.2  is  satisfactory.  The  bounds  L1Q  are,  in  fact,  very 
modest. 

For  this  purpose  (2.2-3c,f)  are  extended  as  follows.  Let 

$ - « (l+L)  - [(1+P  G +***+P  •••P.G.)  +P  G ] (1+L) 
l nn  n 2 2 a 1 

- < <1+L>  ■ ^2(1+L1)  (1+L)  ' 


(2. 3-la, b) 
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P G,  y G1 

a 1 , X I 


Ll/P. 


N " j£  " ££(1+1^)  (1+L)  (1+Ll)  (1+L) 


. (2.3-lc) 


L./P  in  the  crucial  hf  . Hence  to  minimize 
X A 


I X/N. 


at  fixed  L , 


minimize  IlJ  . But  Lj  must  cope  with  the  uncertainty  in 


P ,P  ,»»»,Pn  ignored  by  L 


However,  if  G2  can  cope  with 


P ,»»*,P  then  L need  only  cope  with  P,  . L is  designed  accord- 
2 n 1 xi 

ingly  and  denoted  by  L°  to  indicate  its  neglect  of  Pj,***  ,PR 
uncertainty.  So  now,  P°  of  (2.2-6)  is  replaced  by 


.1  A 


P,  P,rt‘**P  +p 
1 20  nO  b c 


1+P  _G  +• • *+P  • 

nO  n nO 


’P20P1G1 


V.oV’c 
»20  IWlV 


10 


(2. 3-2a,b,c) 


with  v*  -ire*)  - { } 


instead  of  - (a  of  (2.2-5a)  . In  (2.3-2)  L^g  is  the 

nominal  L°Q  = I*1Q  (cf  L°  = LQ  ) and  Q 20  ”©2  at  nominal  Pio  ' 
for  i * 1, . . . ,n  . 


Lq  was  designed  to  handle  P°  with  its  ‘\f°  , but  now  it  must 

handle  P*  with  its  V1  3 v*  . What  are  the  bounds  B (u)  on 
© 1 

1*10  (j“)  80  that  the  original  LQ  remains  satisfactory?  This 

question  may  be  answered  by  simply  trying  values  and  checking  if 

(2. 2-1, 2)  are  satisfied.  It  is  found  that  the  B^ (w)  are  decisively 
influenced  by  the  intervals  I , I , I of  Lrt  in  (2.2-8).  The 
results  are  stated  here  and  their  explanation  in  Appendix  1,  at  the 
end  of  this  chapter.  Note  the  similarity  of  these  bounds  to  those  on 
the  first  inner  loop  in  cascade  design  — Chapter  1,  Figure  1-9. 


i 


Nature  of  bounds  on  L^Q  . 

(1)  For  a)  € I = [0,90]  , { oi)  are  upper  bounds,  i.e., 

| L^Q(ja))|  must  be  < some  value,  which  is  a function  of 

ArgL1Q(jo))  — Figure  2-7a. 

(2)  For  id  6 I_  = [90,330]  , B,  to ) are  lower  ones  precluding 

B 1 ~ 

L1qs  0 (Figure  2-7b) . 

(3)  For  (o  € Ic  = [330,®]  , B^  to ) in  Figure  2-7b  are  closed 

curves  in  the  Nichols  Chart  which  tend  to  a vertical  line 

/bl  \ 

B of  length  — J at  ArgL  = - ir  . 

In  “ 


(?) 


Just  as  in  the  design  of  LQ  , so  the  optimum  would  lie  on 

B^to)  at  all  (o  but  is  in  practive  approximated  by  a rational 
function  — Figures  2-4,  2-7.  One  may  define  intervals  of  similar 

t0  TA  ' XB  ' IC  °f  L0  , i.e.,  in  Figures  2-4,7:  Ilft  = '* 

X1B  " l“lJ'"lx)  ? Ilc"tol*'*)  ’ Here  J1B  * ° because  Bi«  has  zero 
width.  In  practice,  one  would  likely  (in  addition  to  (2.2-2),  assign 

bounds  in  Figure  2-1,  on 


D.  (1+L) (1+L  ) • • • 1+L  ) 
i In 


(2.3-3) 


leading  to  9 . values  and  finite  width  B.  and  larger  I.  . Such 


finite  B are  easily  added  in  Figures  2-7,8,  but  are  omitted  here 
lH 

for  simplicity.  L. n(jw)  is  sketched  in  Figure  2-7b  and  is 


Lio- 


°-14  (1+0OO)(1  + 2fe> 


(1+S5o  ,t1  + ifes4'(33o)2  + (nfe> 


— > 2 f 

inn'  J 


The  constraints  on  the  inner  loops,  due  to  2.3-3,  are  considered  in 
detail  in  Chapter  3,  Section  4. 
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2.3.2.  Second  inner  loop  L^  . 

The  above  discussion  is  repeated  for  L 2 , but  now  P2  uncert- 
ainty is  included,  with  (2.3-1)  extended  to 

" ^2(1+L1)  (1+L)  = l(1+PnGn+“*+Pn**'P3G3)+Pn“’P2G21  (1+L1)  (1+L) 


6 (4^1+P  • • *P_G_)  (1+L  ) (1+L)  ^•^(l+L. 


(2. 3-4a,b,c) 


3 ‘n * 2^2'  v -1  TiJi ' ^ (1+L^)  (1+L)  , 


A VP2G2 


-X  G2 


VPn“*P2 


"2_  ~S\  ■ (1+L2) (1+1^) (1+L) 

" L2^P2P3  tbe  crucial  hf.  To  minimize  the  latter  it  is  best  to 
let  L2  handle  ?2  uncertainty  only,  leading  to  (cf  (2  2a)) 


A, A P P +p 
1 2 aO  b c 


,2  A 

■>  s - — ■ ■ 

e 1+P  _G  +***+P  „p„  G_ +P  • ••&  P G +P  •••p  ppg 
nO  n nO  30  3 nO  30  2 2 nO  30  2 1 1 


(M  P oPk+P  ) 
1 2 aO  b c 


& 30 11+X2L20+X1X2L10 (1+L20)  1 ' 2 P 


2_ 

20 


1 P 


with 


10 

(2 . 3-5a,b,c) 


Y2  =m  p2}  «{ — 1 


instead  of  the  smaller  uncertainty  set  I/1  of  (2.3-2c) 


The  next  step  is  to  find  B2 (io)  , the  bounds  on  L2Q  ■ L2Q  , so 

that  Lq  , L1q  designed  for  , remain  satisfactory  for  ir2  . The 

resulting  B„(u>)  are  similar  to  B„  (w)  : upper  bounds  in  I,.  = 10, u,  ) , 
* 1 1A  lx 

lower  ones  in  X and  closed  curves  merging  into  a B2h  etc.  of 

/b2\ 

length  I — I — see  Fiqure  2-8.  The  explanation  is  given  in  Appendix 

XVdb 

2.  A rational  L20(s)  *s  found  which  satisfies  the  B2  (id)  , as  shown 
in  Figures  2-8,  2-5  with 
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Figure  2-8.  Bounds  B2  (to) 


on  L2Q(jw)  - upper  in  Ilft  . 


2.3.3.  Third  and  higher  inner  loops. 

One  can  continue  indefinitely  in  this  manner.  The  resulting  L2Q 
(Figures  2-5,8)  has  three  intervals  I2ft  , I2B  , I2C  which  decisively 
influence  the  bounds  on  a L^0  designed  to  handle  uncertainty, 

etc.  The  general  forms  for  the  ^ etc.  are  for  i = l,*,,,n  . 


A -A 


W 


P •••P.G. 
n ii 


i+l 


(1+1^)  (1+Li_1)*“  (1+L) 


(2.3-6a-d) 


'S' 


•X,P  -P.+P  ) 
i aO  b c 
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Note  that  F in  Figure  2-1  is  available  from  Equation  (2.2-3e)  as  soon 

as  LQ  is  known,  by  associating  a nominal  TQ(s)  with  the  nominal 

L.(s)  . But  G is  not  known  until  L , L are  known.  Thus 

vi  n n-1  1 

from  (2 . 3- 6b) 


Ln0<^+1,0 


nO 


“ ^O^nO  ' Gn-1  “ Ln-l,0^0/Pn0Pn-l,0 


with  J&) 


1 and  <0nO  = 1+LnO  ' etC‘ 


2.3.4.  Design  completion. 

The  numerical  design  is  completed  for  the  above  example. 

The  previous  steps  gave  rational  functions  for 

L0  • *.0  ■ ■ L10  ■ P.0°/^20  ■ L20  ■ P2oV>30  "ith 

^30  = 1 because  n =2  . One  starts  from  G 2 = L2o/,P20  ' avai*ahle 

because  L^q  and  P^  = 50/s  are  known.  Then  <#20  Of  (2.3-la)  = 

1 + P2oG2  *s  obtained»  whose  numerator  is 

[s2+1260s+(5.14)106] [s2+7860s+ (1. 90) 10®] [s2+9930s+ (2 . 58) 10®]  . 

The  denominator  is  known  from  that  of  • 

The  next  step  is  to  evaluate  G^  = , L10»*^0  known  and 

PaQ  = 1000/s2  . This  gives  &1Q=  1 + p2oG2 + P20P10G1  ' whose  numer- 

ator  is 

(s+605) [s2+132s+(1.32)105] (s2+1260s+ (5. 10) 106] [s2+1200s+(6.05)  106] 

[s2+1300s+(6.50)106] [s2+7860s+ (1 . 90) 10®] [s2+9930s+ (2 . 58) 10®]  . 

Next  find  G = LR  ^q/Pq  • with  P0  = 2000/s3  . F in  Figure  2-1  is 
obtained  precisely  as  in  single-loop  design  ( [H3] ) . Here  it  was 
chosen  [with  F(0)  • 1 ], 
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. (1.59)105(s+5.05) 

F (s)  - — 2 ~2 • 

(s  +289S+103000) (s  +3.67s+7.80) 

The  system  was  simulated  on  the  digital  computer  with  the  results 
shown  in  Figures  2-9a,b.  Runs  with  very  similar  outputs  are  shown  by  a 
single  curve  labelled  S . There  is  a slight  excursion  beyond  the 
upper  bound  in  Figure  2-9a  for  some  parameter  combinations,  and  for 
small  t of  the  lower  bound.  But  this  is  due  to  inexact  translation 
of  time-domain  bounds  into  equivalent  w-domain  bounds,  because  the 
(D-domain  specifications  are  all  satisfied  by  the  design.  The  slight 
violations  in  Figure  2-9a  are  not  atypical  in  the  experience  of  the 
authors.  As  previously  emphasized,  such  translation  is  only  an  art  as 
yet.  With  a little  more  effort  it  is  possible  to  precisely  satisfy 
the  time-domain  bounds  — this  has  been  proven  ( [Hill ) . But 
generally  such  slight  violations  are  acceptable.  The  maximum  over- 
shoot in  the  disturbance  response  in  Figure  2-9b,  is  25%  which  is  very 
well  correlated  with  y of  2.3  db  specified,  which  gives  23.5%  over- 
shoot for  a second  order  system.  The  extent  to  which  the  design 
results  fill  up  the  bounds,  indicates  a design  very  close  to  optimum. 

2.4  Generality  of  Structure 

In  the  system  considered,  input  R in  Figure  2-1  and  n + 1 plant 
outputs  are  available  for  processing,  permitting  an  infinitude  of 
n + 2 degree-of-freedom  structures  . The  n + 2 fundamental 
system  functions  are  the  system  transfer  function  T(s)  ■ C/R  and  the 
n + 1 loops  L,L  ,***,L  . In  any  acceptable  structure,  L is  gotten 

by  cutting  the  outer  loop  just  after  the  C sensor,  giving  in  Figure 

2-10,  L - pQHHiH2^Dl  * *0  i " 1 + P2H2  + P1P2H1H2  ’ KeeP+n9  the  first  cut 

0 


and  with  another  cut  after  the  sensor,  gives  - P2PlHlH2  *1- 


« i+p2h2  . 

T ( s) 

is  always  of  the 

form 

T ■ tyL/  ( 1+L)  , * 

independent  of 

Pi  ' 

- 1/H 

in  Figure 

2-10. 

The  design  technique 

provides  T and 

the 

nominal 

Li0  from 

which 

the  compensations 

(of  Figure  2-1)  or  (of  Figure  2-10)  or  those  of  any  other  structure 

are  derived.  The  excess  of  poles  over  zeros  assigned  to  T(s)  , eT  , 

must  be  compatible  with  the  structure.  In  Figure  2-1, 

eT  = eF + eL  = ep + ep + eQ  each  a positive  integer  but  in  Figure  2-10, 

e ■ Ee.  “ e +e  , i * Q,H  ,H  ,P,H  . 

L i T H 1 2 

2.5  Practical  Design  Procedure  and  Trade-off  — Design  Perspective 

Sections  2.2,3  described  a design  procedure  based  on  the  best 

(most  economical)  L , subject  to  the  best  L subject  to  the 

n n-i 

best  L^  , in  turn  to  the  best  L . First  preference  is  given  to  L , 
then  L^  , etc.  This  section  shows  how  Sections  2.2,3  provide  the 
perspective  for  making  reasonable  trade-offs  between  the  loops  early 
in  the  game,  without  a detailed  design.  Design  perspective  enables 
the  designer  to  very  quickly  consider  design  alternatives,  making  it 
comparatively  easy  to  make  his  decision  which  loops  to  use  and  the 
ou  trade-off  value.  The  procedure  is  very  similar  to  that  for  the 
cascade  design  described  in  Chapter  1,  Section  1.5.  In  fact,  Figures 
l-24a,b  are  used  here  in  the  same  manner.  The  display  in  Figure  2-5 
is  used. 

Design  Perspective  Procedure. 

1.  The  first  step  is  to  design  the  single-loop  LgQ  to  hadle 
the  entire  uncertainty  of  P in  (2.2-3b).  Let  aQ  be  the  overdesign 
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margin  of  Lq  , then  the  outpr  loop  Lg  of  a multiple-loop  design,  is 


obtained  by  shifting  the  UHF  characteristic  of 
the  net  saving  due  to  the  multiple-loop  design, 

(aQ  ■ 0 was  used  in  our  example)  . 4>g  » 20  log 


Ls0  uPward  by  *s  ' 


here  39.4  db  - aQ 

XlxaaVbd 
aabb  + bc  J'a° 


db  - 


see  Figure  2.2-5. 


2.  Find  (at  which  L^g  has  its  peak  value)  by  placing  a 

transparency  of  Louter  (of  Figure  l-24a)  over  LQ  of  Figure  2-5  and 
noting  the  arrow. 


3.  Locate  (the  peak  position  of  L^g  ) at  , with  value 

of  (“  -10  db)  obtained  from  Appendix  3 of  this  chapter.  Point 

is  lLio^wml^  I * Draw  a horizontal  line  in  Figure  2-5  at  A^ 
magnitude  with  A^  ■ min  (B^)  ~ (“  “36  or  -56  db  for  the  example) , 

where  min  (B^)  is  the  point  of  minimum  value  of  . Lay  L^g  (IF) 

of  Figure  l-24a  over  Figure  2-5,  so  that  Q coincides  with  . Find 
the  intersection  point  (or  CJ  ) of  this  LiQ(IF)  with  Ax  (or 

A|  ) horizontal  line  in  Figure  2-5. 

4.  Pick  the  LiQ(HF)  curve  in  Figure  l-24b,  according  to  the 
value  being  used  for  L^g  . Lay  L^g(HF)  on  Figure  2-5,  such 

that  C coincides  with  C±  (or  CJ  ).  L1Q  consists  of  LiQ (IF)  of 
Figure  l-24a  in  the  intermediate  u range  and  of  LiQ(HF)  of  Figure 
l-24b  in  the  high  frequency  range.  Use  the  pattern  of  LiQ(HF)  to  the 
left  of  C to  obtain  a smooth  continuous  curve  for  L^g  . 

5.  Steps  3,4  are  repeated  in  order  to  determine  L^g  . Use  the 
arrow  on  L^g(HF)  of  Figure  l-24b,  to  locate  u>m2  . Use  (2.A3-6)  to 
obtain  M2  , giving  a new  point  Q2  in  Figure  2-5.  Then  lay  LiQ(IF) 


of  Figure  l-24a  on  Figure  2-5,  so  that  Q2 = Q • A horizontal  line  of 

value  » Min  (Bh2>  - a2  (at  -22  db  in  the  example)  is  drawn,  etc., 

where  Min  (B^)  tAe  P°int  of  minimum  value  of  B^2  . The  entire 

process  is  repeated  until  all  the  loops  are  exhausted.  The  dashed  lines 

in  Figure  2-5  are  the  results  of  using  this  perspective  technique.  They 
are  in  excellent  agreement  with  the  actual  detailed  design  (solid  lines 
in  Figure  2-5)  obtained  by  using  the  computer  to  obtain  the  inner  loop 
bounds,  finding  a rational  LiQ(s)  to  satisfy  these  bounds,  etc. 

6.  After  each  is  obtained,  it  is  a good  idea  to  sketch  the 

effective  P values  to  use  for  the  sensor  noise  effect.  Thus  after 
LsQ  is  obtained,  sketch  I pq I = I paoPbO  + Pc0 I in  Fi9ure  2“5-  If  there  is 
little  sensor  noise  amplification  (|Ls0/po!  not  lar9e  over  a large  w 
range) , there  may  be  no  point  in  using  more  feedback  loops.  After  LQ 
has  been  obtained,  it  is  easy  to  see  the  saving  in  sensor  N noisr 
effect,  by  using  . Sketch  |pag|  to  see  N1  sensor  noise 

effect  (at  |L10/pa0|).  Similarly  |L2c/P2ol'  9ives  the  hf  N2  sensor 
noise  effect. 

The  designer  has  to  decide  which  sensor  points  to  use  and  the 
corresponding  or  trade-off  values.  Design  perspective  enables  him  to 
try  out  various  designs  very  quickly  and  easily,  and  thus  arrive  at  a 
suitable  trade-off,  after  which  he  can  proceed  with  a careful  detailed 
design. 


2.5.1.  Bandwidth  propagation  and  similarity  with  the  cascade 
plant  structure. 

Let  the  bandwidth  BW(L^)  be  arbitrarily  defined  as  that  at  which 
Li0  achieves  its  final  asymptotic  slope:  u^0  for  Lq  , w2l  for 
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L10  ' ub2  for  L20  in  2-5.  BW(L^)  increases  with  i . This 

phenomenon  occurs  in  precisely  the  same  manner  in  the  cascade-system 

(Chapter  1,  Figure  1-23).  The  relations  between  the  LiQ  , the  role 

of  * the  sensor  noise  effects  and  trade-offs,  etc.  are  very 

similar  in  the  two  structures.  However,  the  values  of  L..I 

10  max 

are  different  and  the  derivation  is  more  difficult  here.  Here,  at  each 
new  stage,  one  must  use  a moro  complex  form  of  . In  the 

cascade  system,  the  step  from  i to  i + 1 is  identical  to  that  from 
i-1  to  i . But  the  final  results  are  remarkably  similar. 

In  Figure  2-5,  BW(L  _)  • ui  is  comparable  with  BW(Lon)  at  ui  , 
nu  7.2  ou  zs 

a little  larger  due  to  the  extra  few  db  of  gain  margin  needed  per 
section.  Thus,  the  final  cut-off  frequency  for  a single-loop  design 
is  comparable  to  that  for  a multiple-loop  design,  but  they  are 
associated  with  different  loops  so  there  can  be  a qreat  improvement  in 
sensor  noise  effect.  Thus,  in  Figure  2-5,  (Z„-Z  )db  ■ -22  + 89  • 67  db  , 
while  |p_J  ..  - IpJ  „ - 127  db,  an  improvement  if 

20  (Id  0 C1D 

| N0/n I < 127-  67  ■ 60  db.  In  practice  it  is  reasonable  to  assume  that 
tho  plant  power  levels  and  with  them  the  sensor  noise  levels  increase 
in  proceeding  from  input  to  output.  Tho  design  procedure  is  highly 
transparent  |>ermittinq  a good  estimate  of  the  optimum  division  between 
the  feedback  loops,  without  a detailed  design. 

2 . V 2 . High- frequency  unrot ta inty . 
learly.  multiple  loop  design  ran  be  highly  superior  to  sinqle- 
. ^ p ant  uncertainty.  The  linearized  plant  model  is 

n m .....  t nun  I meat  plant  about  an  operating 


mm  Iu»  iff  1l f 
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operating  points,  e.g.,  in  flight  control  ([01])  where  values  > 1000 
have  been  reported. 


It  has  been  proven  that  in  a large  class  of  linear  and  nonlinear 

time-varying  uncertain  plants  the  latter  can  be  represented  for 

synthesis  purposes  by  an  equivalent  linear  time-invariant  uncertain 

plant  set  P^fs]  ([Hi,  Hll]).  The  set  equivalence  is  exact  with 

respect  to  a prescribed  acceptable  plant  output  set.  Linear  time 

invariant  design  applied  to  the  Pe^[s]  problem  is  guaranteed  to 

work  for  the  original  nonlinear  problem.  A nonlinear  plant  with  no 

uncertainty  can  thus  qenerate  larqe  hf  uncertainty  in  P f s ] , e.g., 

eq 

3 3 

consider  y = k x , x the  input  and  y the  output.  Suppose  a fairly 

3 -t  3 

linear  response  is  desired  for  y = A (1-e  ) , A £ [0.5,5]  . To  find 


P [s]  , evaluate  77^7- 
eq  X(s) 


6kA 


eq  (s+2) (s+3) 


in  this  case.  Since 


A € [0.5,5]  , the  hf  gain  of  P varies  by  a factor  of  100,  due 

eq 

2 • 1/3 

to  A . For  a simple  dynamic  example,  consider  y + By  sgn  y = kx  , 

.3 


giving  P 


6kA 


eq  (s+3) [BAs+6A3+2BA] 
factor  of  100. 


2 2 

-»  6kA  /Bs  at  hf  , with  uncertainty 


2.6  Justification  of  Assumptions 

General  Plants.  This  section  is  devoted  to  the  justification 

ei 

of  simplifying  assumptions  in  2.2,3.  One  was  use  of  Pi  =»  ki/s  for 

all  a)  , not  just  in  hf  where  it  is  applicable.  Recall  in  Section  2.2, 

the  first  step  was  to  find  the  smallest  {P°}  of  (2.2-6a),  by 

minimizing  over  G, , • • * ,G  and  the  values  of  P, n,‘",P  _ . Suppose 
in  10  no 

P^»k^/(s+q^)  with  k^  , uncertain.  This  minimization  problem  is 
extremely  difficult  at  medium  a>  . Fortunately,  it  makes  little 
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difference  if  it  is  not  done  at  all.  The  reason  is  that  which  made  the 


inner  loops  unnecessary  for  w < wx  * in  cascade  design  (Section  1.4.4), 


and  in  this  chapter  unnecessary  in  IB  , L_-.  in  I..  , etc.. 


i.e.,  under  certain  conditions  there  is  little  difference  in 


I L_  = GP_ I . needed,  whether  {P}  = set  S,  or  set  S„  <<  S, 

1 u u 1 min  1 2 1 


Figure  2-4,  suppose  that  instead  of  AB  (A  at  X-  ) , the  much  larger 


uncertainty  set  is  ABEFG  with  E,F  extending  even  to  °°  . LQ  at 


X2  results  in  almost  the  same  AHn|T|  for  both  (23.85  db  instead  of 


23.5  db) . 


It  is  therefore  concluded  that  for  oj  < , { P } of  (2.2-3b)  be 


used  for  LQ  design,  just  as  in  Lg  design.  P°  is  used  only  for 


oj  > and  it  is  assumed  that  in  this  hf  range  P . is  well  approximated 
s . 3 

by  k ,/e  ■*  . This  has  been  verified  for  several  numerical  examples. 


e.g.,  for  n=l  with  P = k /(s+q  ) , P,  = k./s  , 


P = k 


;/s(s+qc)  , ka  € [1,400]  , kb  € [1,60]  , k,  € [1,200]  , 


q^  € [0,5.2]  , all  independently  uncertain.  The  maximum  difference  in 


the  two  B(id)  is  only  three  db  even  though  the  difference  between 


( P } ard  fp°}  is  ^<40  db.  If  this  conclusion  is  incorrect  for  an 
e 


unusual  case,  then  it  is  also  likely  that  the  obligation  on  and 


I will  be  greater  too.  By  using  { P } in  medium  w , one  is  certain 


that  the  obligations  on  L in  I will  be  negligible,  as  in 


Section  2.2.  The  simple  and  transparent  forecasting  of  Section  2.5  may 


then  be  used.  If  these  indicate  less  than  desired  saving  in  sensor 


noise  effect,  then  one  can  return  to  check  if  greater  saving  is 


possible  with  P in  I 
e A 
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Another  assumption  in  Section  2.2  was  bc/ac  > 1^/^  • If  the 
opposite  is  true  then  minimum  Lgth  {P°}  is  at  k = b of  value 

Q a a 

r ^ + bc  1 

(Wbc,/(baab+ac)  because  d[kaaa--acJ  / dka  * ° as  shown  in 
Figure  2-2c.  There  exist  a set  of  g^  , e.g.,  g2  = ...  = gn  = 0 , 

% bb 

— $ g.  $ — , which  achieve  this  and  the  procedure  is  precisely  the 

cl  1 D 

C C 

same  as  before.  A third  assumption  is  that  ea+eb  = ec  giving  (2.2-4) 

with  {P  } in  hf  a vertical  line  in  the  Nichols  chart.  Suppose,  how- 
e 

ever,  e +e,  = e +t|»  , i p an  integer.  Then  in  (2.2-5b),  after  proper 
d o c 

\h 

factoring  out,  b is  replaced  by  b (ju)  , a is  replaced  by 
c c c 


ac  ( j uj)  ^ . The  template  vT' 


is  no  longer  a constant,  but  a function 


of  a)  . However,  if  ip  is  even  it  is  still  a vertical  line,  so  the 
simplicity  and  transparency  of  the  design  procedure  is  retained.  The 
design  of  the  outer  loop  is  only  slightly  affected.  One  proceeds 
precisely  as  before  and  obtains  the  bounds  B(u>)  on  L°  , etc.  It  is 
conceivable  that  P(jw)  =0  at  some  discrete  finite  to  values  at  some 
plant  parameter  combinations,  if  * 2 for  example.  Of  course,  this 
means  that  over  some  intervals  in  the  hf  region,  the  two  parallel  paths 
are  basically  working  against  each  other  — which  is  probably  unlikely 
in  a well  designed  plant.  But  if  this  is  so,  the  specifications  on 
T(jw)  and  C/D(jto)  must  allow  for  it.  Of  course,  the  modification  of 
b , a by  (ju)1^  persists  in  the  inner  loop  design  considerations 
(2.3-2c,  2.3-5c,  etc.).  However,  for  <p  even,  the  design  transparency 
is  not  affected.  The  arguments  given  in  Appendices  1,  2 are  applicable 
as  before,  inasmuch  as  they  apply  at  each  a>  and  are  based  on  the  two 


terms  in  the  numerators  of ' 


in  (2.3-2c,  2.3-5c,  etc.)  being  real. 


Separate  treatment  is  required  for  iji  odd  as  it  is  a much  lengthier  problem. 


! I! 
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Another  assumption  was  that  the  disturbance  attenuation  was  a 
minor  problem,  dealt  with  by  (2.2-2,  2.3-3).  The  procedure  is  basically 
the  same  if  it  is  a major  problem,  for  then  CD  «■  C/D  must  satisfy 
|Cp(jw) | s y(w)  no  longer  a constant,  over  {P}  . This  can  be  trans- 
lated into  bounds  on  L0(ju>)  . The  more  stringent  of  BD(w) 

and  of  B (w)  due  to  (2.2-1),  is  used  but  thereafter  the  design 
procedure  is  the  same. 

2.6.1.  Unstable  and  nonminimum-phase  plants. 

Minimum-phase  plants  were  assumed  in  Sections  2.2,  2.3  for 

simplicity.  But  clearly  the  design  procedure  applies  so  long  as  the 

Li0  exist  which  satisfy  the  B^fw)  . Consider  LQ  first.  It  must 

handle  {P°}  giving  (Section  2.2)  a single-loop  problem.  The  latter  is 
6 

solvable  if  {P°}  contains  open- loop  poles  whose  range  of  uncertainty 
includes  part  of  the  right-half  as  well  as  the  left-half  plane  [H3,  H6] . 
If  however,  {p°}  includes  nonminimum-phase  elements  then  LQ  exists 
only  if  the  performance  specifications  are  compatible  with  the  now 
limited  bandwidth  of  . The  same  conclusions  apply  to  the  inner 
loops.  Again,  right  half-plane  P^  poles  pose  no  problem,  but  such 
zeros  impose  limitations  on  . 

2.7  Summary 

For  a class  of  feedback  systems  with  large  uncertainty,  a 
multiple-loop  design  can  result  in  sensor  noise  sensitivity  much 
smaller  than  in  a sinqle-loop  design  satisfying  the  same  specifications. 
The  designer  can  divide  up  the  feedback  burden  among  the  loops  in  a 
sensibly  optimum  manner,  wherein  the  uncertainties  of  the  plant 


J 


sections,  their  levels  and  associated  sensor  noise  sources  play 
important  roles.  An  important  feature  of  the  design  techniques  is 
its  transparency.  In  return  for  learning  the  mechanics  of  sensitivity 
reduction  in  the  language  of  frequency  response,  there  is  gained 
excellent  insight  into  the  trade-offs  between  loops  and  the  overall 
cost  of  design  in  terms  of  bandwidth  and  noise  sensitivity  — even 
without  performing  the  detailed  design.  It  is  interesting  that 
although  the  derivation  in  this  chapter  is  much  more  difficult  than 
for  the  cascade  system,  the  results  are  remarkably  similar.  This  is 
revealed  by  the  great  similarity  of  the  Design  Perspective  procedure 
of  Section  2.5,  to  that  for  the  cascade  structure,  of  Chapter  1, 


Section  1.5 


Section  2.3  presented  without  explanation  the  bounds  B1  (u>)  in 
terms  of  the  intervals  Ift  , I , Ic  of  LQ(ja))  . The  explanation  is 
available  by  considering  the  uncertainty  of  variation  set  (2.3-2c) . 


IT1  ■ { ^7  }^V°  ■ 


+ k } of  (2 . 2-5c) . 
c 


is  the  line  AB  in  Figures  2-Ala, b,c,  whereas  \f 1 is  the  larger  set 

ABC  D , a function  of  L, „ and  X , . The  point  A is  always  the 

j j 10  lmax 

nominal  LQ  , X1=l  , ^ , kc  = % irrespective  of  the  value  of  I<10 , 

because  that  is  the  objective  of  the  B^w)  . Attention  is  focused  on 
the  range  -it  < ArgL1Q  < 0 . The  following  properties  of  ^*are  important 

(PI)  In  Figure  2-Ala,  as  |l^q|  is  increased  at  fixed  Arg  , 

boundaries  BCi  , AEL  shift  downward  - compare  BCjCj  at  0 db  with 

BC  C'  at  -20  db  and  BC,C'  at  -40  db;  and  similarly  the  AD . D . 

2 2 11  ii 

(P2)  For  fixed  L1Q  , the  effect  of  increase  in  X^^  is  extension  of 

the  BCi  , ADi  , i.e.,  widening  of  the  regions  by  decreasing  amounts, 

to  a maximum  of 

i „ -1  m sin  6 i 

A9  = 0 - tan  — 

1 1 + m cos  0 1 

as  Xlx*«  , where  L^g  * m /0  . This  can  be  seen  as  follows  from 
Figure  2-Ald. 


A0  = tan 


X,  m sin  0 
lx 

1 + X,  m cos  0 
lx 


-1  m sin  0 

i ^ r 

1 + m cos  0 


-+  0 - tan 


-1  m sin  0 
1 + m cos  0 


This  effect  of  large  X^  is  important  in  explaining  the  nature  of 


Figure  2-Ala.  Explanation  of  nature  of  (w)  in  — family  of 

at  fixed  ArgL^  ■ -90°  , for  various  |liq| 
and  X,  . 


(P3)  For  given  X and  | L | , T/*  at  /Lir>  = ® the 


image  (about  AB  ) of  l/1  at  /Lin  = -0  • 


The  upper  bounds  of  B (u)  in  I are  explained  by  property 

(PI)  in  Figure  2-Ala . A family  of  2/1  at  fixed  L^g  = -90°  tr*ed 

at  a)  * 40  € I,  , i.e.,  point  A of  is  set  at  L (j40)  = -32  db 

A 0 

/-130°  (from  Figure  2-4)  at  which  Figure  2-1  requires 


& $ 34.3  db  . 

T . 

1 min 

At  w = 40  (2. 2-1, 2)  are  precisely  satisfied.  It  is  seen  in  Figure  2-Ala 
that  at  ArgL^g  = -90°  , |L10I  < -20  db  is  OK  while  lL10l  5 0 db  is  not 
because  A | T | = | — 34—2 . 3 1 = 36.3  db  and  larger  |L^g|  gives  larger 
a|t|  . The  upper  bound  here  is  between  0 db  and  -20  db.  From  a study 
of  the  shape  of  constant  |l/1+l|  loci  on  the  Nichols  chart,  it  is 
seen  that  this  result  applies  for  all  w £ Ift  at  which  ArgLg  £ -90°  . 

In  Figure  2-4,  there  is  a small  interval  in  which  ArgLg > -90°  and  in 
general  there  may  be  a low  frequency  region  where  ArgLQ > -90°  . 

However,  the  final  result  is  basically  the  same,  because  of  the  very 
small  sensitivity  of  the  loci  of  constant  |l/1+l|  on  the  Nichols 
chart  at  | l|  » 1 . 

It  is  worth  noting  that  in  Ift  , if  L^Q  did  not  exist  at  all , then  * , 
i-l,...,n  would  only  be  a much  longer  vertical  line  with  lowest  point 
at  A . From  Figure  2- Ala,  both  (2.2-1)  and  (2.2-2)  would  still  be  satis- 
fied. Thus  for  ui  € I , LQ  designed  for  Pg  , Pg  uncertainty  only, 
automatically  handles  p^***pn  uncertainty  as  well.  However,  L^  is 
needed  in  I , precluding  L. n ■ 0 in  I and  giving  there  upper 
bounds  as  in  Figure  2-Ala.  Similarly  note  that  in  Ift  , B^w)  are 


hardly  affected  by  large  Increase  of  X 


There 


fore,  L could  handle  the  entire  uncertainty  of  P 


place  of  X if  G ■ G 


Figure  2-Alb.  Explanation  of  nature  B^(w)  in  IB 

at  fixed  ArgL  ■ 0 . 


family  of 


Property  (PI)  also  explains  in  Figure  2-Alb  the  lower  bounds  in 
1B  . At  |liq|  -n^  , penetrates  into  M<2.3  db  , violating 
(2.2-2).  Thus  in  Figure  2-Alb,  at  ArgL1Q  ■ 0 , I L10 1 mln  “ m2  • In 
range,  (2.2-2)  easily  dominates  so  there  is  no  danger  of  violating 
(2.2-1)  (cf  Figure  2-Ala)  except  possibly  at  very  large  |l  | , which 


-86- 


would  not  be  used  anyhow.  Here  too,  X^  could  be  increased  to  • with- 
out affecting  B^(u>)  — recall  (P2) , the  effect  of  large  X^  on  i/1 
in  Figure  2-Ala  and  the  critical  factors  in  Figure  2-Alb.  Thus  there 

is  no  need  for  L_,...,L  in  as  well.  (Pi)  also  explains  in 

d£  n b 


Figure  2-Alc  the  upper  and  lower  bounds  in  I 


Thus  at  /L^q  = 0 , 


) L.  n | must  be  either  < m_  or  > m . From  (P2)  the  width  of  \P 


'10 

is  < | ArgL 


10' 


2 5 

Hence,  Figure  2-Alc  shows  that  as  u increases  in 


Ic  , the  value  of  -Argl»10  for  which  all  lL10l 


are  acceptable , in- 


creases steadily,  explaining  why  the  B^tu)  closed  curves  shrink  to 

/bl  \ 

B in  Figure  2-lb.  B length  is  l — I because  at  ArgL  = -tt 
1H  1H  \ S]L  10 

(say  at  1,000  here)  1 + X^L^  = 1-X  |i.^q|  with 

al  1 

lLlol  < 1/Xlmax  “ b^  ' othervise  V'1  extends  in  length  to  “ and 
being  360°  wide,  must  intersect  with  the  forbidden  |t~t|  <Y  m 2.3  db 

i.T*L 


regions  located  at  ArgL  * ± nir  , n = 1 , 3 , . . 


This  is  also  seen  from 


(2.3-3),  for  let  L.  *X,L.  _ , and  P * P _ for  a i . Then  at 
i l io  a au 

ArgLlri>*-n  , |l  | < l/\, is  essential,  otherwise  (c^/D 


iO  ' 1 iO 1 

infinite  at  X, 

imax 


imax 


i'  i> 


is 


Increase  of  X^  affects  the  bounds  at  m^  , requiring 
l^iO  I < , 0 < 1 , but  not  the  lower  boundary  at  m,.  . In  Figure 

2- lb  it  is  seen  that  L1Q  lies  on  the  upper  part  of  B^ (u)  for  most 
of  Ij^  , so  L^q  designed  to  handle  P^  only,  can  also  cope  with 


V 


if  V 


■ L ■ 0 . 
n 


I 


Chapter  2,  Appendix  2 — Bounds  on  Second  and  Higher  Inner  Loops 

The  function  of  L^q  is  to  guarantee  that  L^0  is  satisfactory 

despite  its  design  on  the  basis  of  P^  of  (2.3-2a) . It  was  seen  in 

Appendix  1 that  for  cd  € Ift  , I and  part  of  Ic  , L10  suffices, 

i.e.,  may  be  zero.  This  is  so  only  for  u)  < u>  , at  which 

20 

/L,  ■ -it  . In  this  case  0 , * 0 , so  to  . is  denoted  by  m, 

/ 10  ml  xl  lir 

It  was  notod  also  that  ^ < — ' 30  L2  is  needed  for 

lmax 

m - iu,  _ . Hence,  L_ft  = 0 is  impossible  in  I and  it  is  not 

1 TT  lA 

surprising  that  the  B. (u)  there  are  upper  bounds  (recall  in 


Appendix  1 precisely  the  same  situation  for  L1Q(ja>)  in  1^  ) . At 

o)*u)  € I.  , 1 + X,  L, . = e>0  (.38  in  the  example),  so  the  denom- 

lir  IB  lmax  10 

inator  (2.3-5a)  of  P^  is  (1"X2+eA2+eX2L20)  ' and  f°r  ifc  * 0 


at  ArgL2Q  = 0 , U*20l>  — 


= 1.4  here.  So  there  is  a lower 


bound  on  iL^fja)^)  | which  is  a function  of  ArgL^  . 

2 

To  find  B_(to)  in  I,_  it  is  necessary  to  use  P of  (2.3-5a)  in 
2 ic  e 

1 2 

place  of  P of  (2.3-2a).  It  is  convenient,  however,  to  express  P in 
e e 

1*  1 

terms  of  P , defined  as  p with  X,X.  replacing  X,  , because 
e e 1 2 l 

IV  } is  easier  expressed  in  terms  of  Vi*?}  , while  'IT { P1  } is 
e e e 

easily  gotten  from  1/{P^}  shown  in  Figure  2-Ala  by  letting  X.  > X, 

^ e 1 lmax 

From  (2.3-5a)  and  replacing  X^  in  (2.3-2a)  by  X^X2  ' 


(WW 


(2.A2-1) 


in  Figure  2-A2,  as  follows.  Let  °Qa=^1X2L10  < ArgL^  < “n  in  IBl'ICl  ^ ' 

QV-1  , QDi-a  , |D.V|  - | aL2Q|  , ArgE.D.V  - ArgE.D.^  - ArgL20  , 
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DiCi*A2DiV“  A2aL20  ' 80  0V  * °Q  + QV  - XiA2L10  + 1 ' 

QC 


QC, 


QC.  + D.C. 
i 11 

QD.  + D.V 


a + aX  I, 

2 20 


1 

2 20 


a + aL 


20 


1 + L 


and 


20 


1 +X->L™ 

0C.  = OQ  + QC.  = XXL  + yTl  ~ ' ?ivin<*  (2.A2-1). 

20 

Figure  2-A2  was  sketched  for  u>  = 2000  , X^  = 40  , X2  = 10  , at  which 
(Figure  2- 7b)  /-230°  , LQ = -127  db  /-430°  , for  assumed 

ArgL_  * -117°  constant.  The  D.  describe  an  arc  of  a circle  as 

20  l 

1 1-20 1 vatied,  as  do  the  drawn  for  X2  = 10*=X2max  , i.e., 

* 10  D^V  . Clearly  for  |l2q|<<1  , 0V/0C-*1  and  for  |l2q|>>1  , 

1 0V/0C | < 1 , so  such  1 1*20 1 are  acceptable.  Obviously  9 < X2max  ' 

3 resulting  C^  circle  passes  through  0,  giving  infinite  OV/OC^  and 

the  resulting  L^{P^}  passes  through  M=2.3  db.  Thus,  3 upper 

and  lower  bounds  in  this  w range.  As  w increases,  lL10l  and  its 

angle  decrease,  so  the  arc  C^C_.  • • • does  not  extend  to  0 in  Figure 

2-A2  and  any  |l2^ | is  acceptable.  Hence,  the  B2 (o>)  tend  to  a line 

B at  -it  , from  0 to  (a./b.)  ..  . B_(u)  are  shown  in  Figure  2-8, 

Zn  Z Z ab  Z 

including  a sensibly  optimum  L2Q(jo))  with  its  intervals  I2ft  , , 

J2C  * 


For  the  third  inner  loop  (if  n>  2 ),  P is  needed  and  there 

e 

is  an  analogous  situation  with  respect  to  I . At  I , w (at 

ZPi  ZB 

which  ArgLjQ  = -it  ) is  very  large  (~6500)  and  as  before,  there  is 

a lower  bound  on  L^Q  at  . For  o>  >.  , |lq|  , |liq|«1  , so  in 

(2.3-6d),  Denom.  (P^) ” 1 + X,L  + X,X_L_  (1  + L,_)  , similar  to  Denom. 

e J oU  J £ £U  ju 

2 3 2* 

(P  ) if  i is  replaced  by  i - 1 . P /P  similar  to  (2.A2-1)  is 
e e e 


I 
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obtained  giving  a figure  similar  to  Figure  2-A2  and  analogous  B^(u)  . 

The  process  is  continued  to  B-(co),...,B  (u>)  . 

4 n 


Case  8,^0. 

mi-4 

The  above  is  all  based  on  0 . = 0 , i£l.  For  then  u>  =(u. 

mi  xi  iir 

But  if  0.^0,  there  are  demands  on  L._  for  w > <o  ... 

mi  10  x,i-l 

For  (i)  > u)xl  , the  constraint  on  I^g  is  that  due  to  ^ ^'*2x^' 

the  outer  loop  does  not  cross  the  Yg  locus.  Or,  alternatively,  the 

inner  loop  does  not  cross  the  y^  locus.  For  0ffl^ = 0 , the  y^  locus 

is  merely  the  single  point  -1.  From  Equation  (2.3-lb) 


P1P2G1 
1 + P2G2 


(case  n=2  ) 


P1G1L2 

G2(1+L2) 


which  is  of  the  same  form  as  in  the  cascade  problem.  Therefore,  the 
constraint  on  L2  here  is  similar  to  that  on  L2  in  the  cascade  case. 
The  formula  for  the  peak  value  of  |L2I  ' i-e"  | !*2 < j a>m2 ^ I used 

Design  Perspective  is  that  used  in  the  cascade  problem  (Figure  1-2.2). 
This  is  so  for  the  other  inner  loops  also. 
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Chapter  2,  Appendix  3 - Estimates  of  the  Peak  Values  of  the 

Inner  Loops 

2.A3-1  — Tha  first  lnnar  loop  L^0  . 

In  Figure  2-A3a,  the  frequency  u>ml  ia  approximately  at  point  A 


with 


£ 


(ui_, ) CM  -225°  , also  ArgL1Q(jwml)  « -90°  , and  the  bound  on 


0 'ml 


L (jui  ) corresponds  to  the  template  tangent  to  the  y locus,  say  at 
10  ml 

point  E.  BEC  is  the  locus  of  point  B at  , kc«bc  , due  to  the 
uncertainty  of  6 U**lxl  • Th®  difference  between  any  point  on  BEC 
and  point  B can  be  written,  from  (2.3-2c),  as 


(X . a b.  + b \ / / a b + b \ . 


(2.A3-1) 


Let  X^X*  at  Point  E,  the  vertical  difference  between  E and  B is 
|f(X*)|  - (aQ  + 0Q)  db  . Both  X*  and  L10  are  not  known,  but  we  have 
found  that  using  X*  ■ X partly  offsets  neglect  of  0Q  . With  these 
approximations,  Equation  (2.A3-1)  can  be  solved  for  lLio^uml^  ' giving 


(for  aQ  in  db) , 


/XlxaaVbc\2  1 
Vaabb+bc,/  1 + 


1+lLlcJ  „ 2 

T"  *a 


* r (a./20)  -.2 
- [ 10  ° ] (2.A3-2) 


x2  |L 

lx' 


10' 


Define  6 


. X,  a b.  + b 
A lx  a b c 

a b + b 
a b c 


and  6 • 6 /6  , then  Equation 

as  as 


(2. A3- 2)  becomes 


M1  " lL10(u’ml)| 


I l-(6as)2 

i ^asV2"1 


(2.A3-3) 


Due  to  the  shapes  of  the  inner  loops  (see  Figures  2-7b,c),  it  is  found 
that  better  results  are  achieved  by  adding  a few  db  (e.g.  1.5  db)  to 
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the  value  obtained  from  (2.A3-3).  In  the  design  example  of  Section 

2.2.1,  X ■ 40  , 6*1  , 6 *10,  gives  6 *0.1  and 

ix  au  s as 

|l*lo(<*»miH  “ 0.26  . Adding  1.5  db,  gives  |liq(u  . ) | at 

20  log  (0.26)  +1.5  at  -10  db  . The  above  calculation  was  used  to  find 

I Ll0(  j<*>mi)  I in  the  design  perspective  sketches  of  Figure  2-5. 

2.A3-2  — Second  and  higher  inner  loops. 

The  case  0 , ■ 0 is  considered  here  (6,^0  was  treated  in 
mi  ml 

Appendix  2).  The  peak  value  of  |l  | occurs  at  u)*u>  at  which  we 

m2 

have  found  approximately  “ -225°  . The  constraining  condition 

is  that  the  template  of  the  outer  loop  (due  to  X^  € [l,X^  ] , i*l,2  ) 
does  not  cross  the  locus  of  y . 


The  nominal  L^jw^)  is  at  point  A^  on  Figure  2-A3a,  which  is 
so  far  down  that  the  determining  condition  essentially  is  that  the 
denominator  of  lq(*i'*2>  **  ® • This  condition  can  be  written  (see 
2. 3- 5a) 


(1  + X2L20)  + (X1A2L10)‘ (1 +I'20)  * ° * (2.A3-4) 

At  w_wm2  » the  nominal  1 1 “ - (X^  + a^^  + 6^ db  (see  Figure  2-A3b) 

% 

where  db  is  the  overdesign  margin  of  , so  (we  drop  0^,  for 

the  same  reason  0Q  was  dropped  previously)  “ “(a^)db  /-225°  , 

,5  A ^0^ai^20^  , and  (2.A3-4)  becomes 


20  ' (/2  « - 1)  . jl  J 

0(1 

with  f replaced  by  ■ the  limiting  case,  giving 


'-■s'l 


I ^ Wlo  - 11 2 + 1 

/ 1/2  «al  - II2  ♦ 1 


(2.A3-5) 


Since 


giving 


"dT^  >°  ' We  UM  X2-X2x  ' 


/(/ 2 6 /X  -1)  +1 

M2  " lL2(“«2,l  " / 2 

2 2 m2  / (/2  6 - 1)  + 1 

al 


(2.A3-6) 


In  the  example  of  Section  2.2.1,  X^-10  ' °1  " 2 * (6^  *1.26), 

“ 0.18  db  , used  in  the  design  perspective  in  Figure  2-5. 

For  n-  3 , the  first  part  of  the  previous  is  the  same,  but  the 
denominator  of  L(X^,X2,X3>  is  now  (from  Equation  2.3-6d) 

1 + X3L30  + X3X2L20(1  + L30)  + X3X2X1L10 (1  + L20)  (1  + L30)  * 

However,  at  <<  ^L20^'^L30^  ' 80  the  above  is  closely 


1 + X3L30  + X3X2L20*X + L30)  ' 


which  is  the  same  as  (2.A3-4)  if  is  replaced  by  * for 

i * 1,2  . Hence,  in  (2.A3-6)  replace  i by  i+1  for  i*l,2  , giving 


wher.  6 . is  the  overdesign  margin  of  L.  in  arithmetic  value. 


07) 


CHAPTER  3 


THE  TRIANGULAR  MULTIPLE-LOOP  PLANT 

3.1  Introduction 

This  chapter  extends  quantitative  synthesis  to  the  plant  structure 
of  Figure  3-1,  denoted  as  the  "triangular  multiple-loop  plant",  again 
with  no  "plant  modification"  (Section  1.3.3).  This  plant  can  be 
considered  as  one  form  of  generalization  of  the  cascade-parallel  plant 
of  Figure  2-1,  in  which  the  cascade  branch  has  only  two  sections. 

Chapter  5 presents  another  generalization.  It  is  not  surprising  that 
quantitative  synthesis  for  Figure  3-1  contains  elements  of  cascade 
synthesis  (Chapter  1)  and  of  cascade-parallel  synthesis  (Chapter  2) . 

The  objectives  are,  as  before  to  satisfy  specifications  of  the  form 
(2. 2-1, 2),  with  minimum  net  sensor  noise  effects  at  the  plant  input. 

And,  as  before,  "design  perspective”  is  sought  which  permits  the 
designer,  with  relatively  little  effort,  to  see  the  important  trade- 
offs and  enables  him  to  decide  which  of  the  available  sensing  points 


The  following  notation  is  used: 


3.2  Design  of  the  Outer  Loo; 


In  Figure  3-1,  with  the  notation  of  3.1-1,  the  system  transfer 


function  is 


where 


The  effect  of  sensor  noise  N.  at  the  plant  input  X , is 


Hence,  to  minimize  this  effect  reduce  1 1,^ | as  fast  as  possible,  etc. 

As  in  previous  chapters,  the  conclusion  is  to  let  L.  handle  the  entire 


the  frequency  at  which  it  turns  the  corner  of 


uncertainty  up  to  u 


its  "universal  high-frequency  boundary 


handles  only  the  uncertainty  in  P 


with 


and  nominal 
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Note  the  similarity  of  this  notation  to  that  in  Chapter  2 (Equations 
2.2-6a,b) . 

P2  here  replaces  P#  ■ PiP2,,*I>n  ot  F^9ure  2-1.  rt  *s  a9*l-n 

assumed  that  in  w > uxl  , each  Pi , is  closely  equal  to  its  hf 

ei1  3 e2 

asymptote  j/s  J . The  value  of  k2  (P2  -*■  k2/s  ) to  be  used  as 

nominal  in  order  to  minimize  the  template  of  P°  , and  the  resulting 

requirement  on  , are  deduced  exactly  as  in  Chapter  2,  Section  2,  so 

the  discussion  is  not  repeated  here.  Thus,  the  design  of  the  outer  loop 

is  exactly  as  in  Chapter  2,  if  P there  is  replaced  by  P_  here. 
10  * * 


3.3  Design  of  First  Inner  Loop  L-, 

The  obvious  step  now  is  to  allow  for  the  uncertainty  in  P2b  ' P2c 
to  be  handled  by  the  loop  via  G2  . For  this  purpose,  let 

P2  * P30P2b  + P2c  “ X2(P30P2b0 +P2cO)  “ X2P20  (3.3-1) 


- 


(P2Plb  + Plc)Gl 


(X2P20Plb  + Plc)Gl 


(1 +PnOGn+ +P30G3)  + P2G2  Ao  + X2P20G2 


with 


•QsO*1  + *2L20* 

’ 1 + PnOGn  + + P30G3 


20 


P20G2 

&30 


(3.3-2a-c) 


(3.3-4) 


Equation  (3.3-2c)  should  be  compared  with  (2.3-2b).  Except  for  the 

o 

slight  difference  in  notation,  the  problems  are  identical.  was 

designed  to  handle  uncertainties  of  P^  , P^c  only  (in  Chapter  2 
Pfa  , Pc  only) , but  LXQ  must  handle  the  uncertainty  of 

pl  A X2P20Plb  * Plc 
1 


(3.3-5) 


which  ia  analogous  to  P*  of  (2.3-2) , with  L2Q  hare  analogous  to 
L1q  of  (2.3-2),  X*  hare  to  X^  there.  The  design  procedure  for  L^q 
(first  inner  loop  here)  is  therefore  identical  to  that  of  the  first 
inner  loop  in  Chapter  2 (denoted  there  by  L10  >'  and  requires  no 
further  discussion. 


3.4  The  Two  Constraints  on  Second  Inner  Loop 
At  this  point,  we  distinguish  between  two  cases.  The  first  is 
that  treated  in  Chapter  2,  where  the  only  specifications  are  those  in 
(2. 2-1, 2).  There  are  none  on  X/D1  of  (2.3-3),  and  therefore  no  phase 
margins  in  the  universal  hf  bounds  B^  . In  this  case,  the  latter  are 
vertical  lines  in  Figures  2-7b,8.  The  inner  loop  designs  in  Chapter  2 
were  all  done  on  this  basis.  In  the  second  case,  there  is  also  the 
specification  (2.3-3)  and  similarly  in  Figure  3-1,  on 


D1  (1-flj) (1+L2) ... (1+Ln) 


(3.4-1) 


with  L_,Lj,...L  defined  later. 


1+L, 


Due  to  the  bandwidth  propagation  effect,  the  peak  of  each 
I * occurs  in  a different  w range,  so  each  factor  can  be 


considered  separately.  A maximum  y^  (analogous  to  y of  2.2-2)  is 

Li 

, which  is  reflected  by  a phase-margin  6^ 


assigned  to  each 


1+L, 


analogous  to  0m  of  the  outer  loop  in  Chapter  2,  and  a finite  width 

. The  specifications  on  each  can  be  different,  because  larger  peak 
values  may  be  tolerable  at  higher  frequencies.  In  this  second  case, 
when  designing  , one  must  take  care  that  both  the  of  and 

of  L>2  are  not  violated.  The  latter  constraint  dominates. 


Similarly,  whan  designing  any  one  must  conside.  the  of 

for  i - 1,2,3,..., j-1  . The  last  constraint  (on  ) dominates. 

The  reason  for  the  above  domination  is  demonstrated  for  the 
cascade-parallel  problem  of  Chapter  2,  by  means  of  Figure  3-2,  where 
the  values  of  the  wx^  , relative  to  the  B^  are  shown.  The 

obligations  on  the  second  inner  loop  (of  Chapter  2)  L 2 due  to  the 
outer  loop  L begin  at  . But  | ( juxl ^ I *s  verV  small  (e.g.  in 

Chapter  2,  Figure  2-7b,  wxl  1000  ),  with  located  a good  ways 

down  the  hf  asymptote  of  LQ  . The  only  danger  to  L (due  to  L*2  ) is 
that  1+L  i 0 . The  resulting  bound  on  l2Q  in  the  Nichols  chart  is 
the  vertical  line  at  -180°  from  0 db  to  -(X2x>  *1'  db  . The  limitation 
at  u»xl  imposed  on  L2  (due  to  B^  ) is  much  more  serious  as  it  is 
similar  to  the  one  imposed  at  wx  (of  the  outer  loop)  in  a cascade 
design,  on  the  first  inner  loop,  due  to  the  bound  B^  on  the  outer 
loop.  This  is  seen  as  follows. 

From  (2.3-lc)  when  uncertainty  in  both  and  P2  is  considered, 

the  first  inner  loop  is 


PaGl 

Qi 


Pn0‘*'P30P2PlGl 


1 + P + 
nO  n 


+ Pn0* * ”P30G3  + Pn0*  * *P30P2G2 


Pn0*'*P30P2PlGl 


( 1 + PnOGn  + * • * + pn0  • 


■P30G3) 


1 + 


Pn0* " *P30P2G2 


1 + P -G  + 
nO  n 


+ Pn0’ * ,P30G3 


A L2P1G1 
“ G2(l+L2) 


(3.4-2) 


Consider  a cascade  design  of  2 sections  in  Figure  3-3.  The  outer  loop 
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L-  - 


pbpaga 


o 1 + P_G_  ' 
B B 


(3.4-3) 


Figure  3-3.  A cascade  system  of  2 sections. 


which  is  analogous  to  (2.4-2)  with 


G =G  p „p  p Pn0---P30P2 

B 2 ' 1 A ' B " l+Pn0Gn+...+P.n...P,rtG 


Hence,  the  bounds  on  the  second  inner  loop. 


n0*  30  3 

(3.4-4a-c) 


P „ P„G_ 

nO  20  2 


"20  1 + P nG  + ...  +P  A...P,„G, 

nO  n nO  30  3 


in  (3.4-2),  such  that  those  on  L 


are  satisfied,  are  precisely  of  the  same  type  as  those  on  * PB0GB  ' 


in  (3.4-3),  so  that  those  on  are  satisfied. 


In  the  case  of  the  third  inner  loop  ( L 4 in  the  notation  of 


Chapter  3) , the  bounds  on  must  all  be  considered.  For 


the  same  reason  as  in  the  above,  those  on  L3  dominate. 


3.5  Design  of  second  Inner  Loop 


Suppose  0m2-°  , so  the  universal  hf  boundary  of  L is  a 


vertical  line.  Only  the  constraint  on  L.  due  to  the  bound  on  L 

3 1 

need  be  considered.  The  uncertainties  in  p^  , p previously 


omitted,  are  now  included,  giving 


The  effective  plant  is 


which  is  comparable  with  P 


(here)  is  set  equal 


to  XX.  (there),  X 


(here)  to  X-  (there) 


design  of  the  second  inner  loop  here  is  similar  to  that  of  the  second 


inner  loop  in  Chapter  2 


However,  in  the  case  of  O 


constraints  on  L, 


these  dominate 


the  design  of  L30  , whose  desiqn  is  then  comparable  to  that  of  the 
first  inner  loop  in  the  parallel-cascade  design  of  chapter  2.  This  is 
seen  by  exhibiting  the  first  inner  loop  L-  in  terms  of  the  un- 


certainties in  P 


and  then 
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L1  (P3c*P3Pa,)G2 


(3.5-8) 


Compare  to  the  first  inner  loop  of  Chapter  2,  whose  effective 

plant  in  (2.3-2a)  is 


P + (P.P,n..  .P  _)P. 
c 1 20  nO  b 


r 

020  [ 1 * — 


<P1P20-“P„0) 


(3.5-9) 


Thus,  the  design  of  L^0  - P3oG3$4o  to  sat*sfy  tlle  bounds  on 
L2  , is  analogous  to  the  design  of  the  first  inner  loop  in  Chapter  2, 
to  satisfy  the  bounds  on  the  outer  loop. 


Higher  Inner  Loops 


For  the  case  of  i-th  inner  loop  design,  if  0.-0  for 

m3 

2 $ j £ i-1  , then  the  effective  plant  P^  can  be  written  as 


?i  _ (XiPiOPlb'l'Plc) 

1 + *i+1Lt+lf0  + XiLi,o(1+Li+l,0)  + ,,,+X2L20a+L30)  (1+L40)  * ' * (1+L1+1 ,0) 

p (xi  p p +p  ) (3*5_10) 

Lth  X*  , Xi  , J+*  3+1  f0 , j - 2 , 3, . . . i 

i+1,0  3 jO 


.i  P1R1  . PjOGj 

L.  -yr  , L.  - , j-2,3,...i 

1 <®i+2,0  j°  £)j+l,0 


(3.5-lla,b) 
( 3. 5-12a,b) 


In  the  case  the  0 . 0 , the  design  of  L ..  is  dominated  by  the 

mi  40 

need  to  satisfy  the  B.  , bound  on  L,  . Write 

n J J 


1 tP3c  * P3b(P4c  + P4bP50)]G3 
L3  " r (P.  + P ..  P__)G.  -i 
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4 <P3e  * P3bP«,S3 


(3.5-13) 


Compare  this  to  (3.5-8),  and  it  is  seen  that  the  design  of  the  third 
inner  loop  to  satisfy  the  bounds  on  the  second  inner  loop,  is  analogous 
to  the  design  of  the  first  inner  loop  in  Chapter  2,  to  satisfy  the 
bounds  on  the  outer  loop.  The  extension  to  higher  inner  loops  is 
obvious. 


3.6  Sensor  Noise  Effects  at  Plant  Input 


These  effects  are  found  here  as  the  nominal  plant  values.  The 
effect  of  was  given  in  (3.2-4),  repeated  here 


-X  L10/P10 


(3.6-1) 


The  effect  of  N,  is  given  by 


I*  2 
"2 


In  general. 


Ni‘£i 


20<1+L10) 


Li0/Pi0 


°2 L20/P20 

J0(1+L20)(1+L10>  " <1+L20)  (1+V 

(3.6-2) 


(1+Li0)...(l+L10) 


(3.6-3) 


3.7  Design  Example 

The  procedure  for  quickly  achieving  approximate  designs  is  fairly 
obvious,  in  view  of  the  above  noted  similarities  of  the  present  design 
problem  with  that  in  Chapter  2.  The  explanation  is  presented  by  means 
of  the  following  example.  The  detailed  design  will  be  presented  first, 
followed  by  Design  Perspective  in  Section  3.8. 


105- 


3.7.1.  Problem  Statement. 

The  structure  Is  shown  in  Figure  3-4a  with 

[ 1 , 40 ] 10  3 ..  [1.18] 100  „ _ [1,60]105 

P4c , ' P3b  ' P3c  2 

s 

[1,18)10  . [2,120] 10^ 

2b  s ' 2c  3 

s 

_ [1,60]  B , [4,800] 106 

Plb  s ' lc  s4 


The  numbers  were  deliberately  chosen  so  that  the  outer  loop  and 
first  inner  loop  designs  here  are  identical  to  those  of  the  design 
example  in  Chapter  2,  Section  2.2.1.  Time-domain  bounds  on  the  accept- 
able step  response  and  its  translation  into  bounds  on  |T(ju>)  | are  the 
same  as  Figure  2-3a,b  of  Chapter  2.  In  the  notation  of  (3.1-1) 


P4  " P4c  * [1,40]103/s  with  uncertainty  - 40 

P3  " P4P3b  + P3c  “ t2*7B0)105/32  , U3  - 390 

’ P40P3b  * P3c  • IS-78U°5^2  ' "I"  » 
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design  of  the  outer  loop  (L10)  is  the  seme,  because  the  design  of  the 
outer  loop  does  not  need  to  refer  to  the  constraints  of  inner  loops. 

3.7.2.  Outer  Loop  Design. 

As  discussed  in  Section  3.2,  let  handle  the  entire  un- 
certainty up  to  . But  for  a)  > o>xl  , handles  only  the  un- 

certainty in  P,  and  P_.  with  all  the  other  plants  at  their 
lc  lb 

nominal  values.  From  (3.7-2),  U°»130-the  outer  loop  uncertainty  of 

the  design  example  in  Chapter  2 given  from  (2.2-5b).  The  condition  in 

Section  2.2  (following  2.2-4a,b)  ^c/ac  ^ b^/e^  becomes  here 

b,  /a,  ■ 800/4  > b,.  /a,.  ■ 60/1  , so  the  outer  loop  of  the  design 

1C  lc  ID  ID 

example  in  2.2.1  can  be  used  here  as  L^q  . on  Nichols  chart  is 

shown  repeatedly  in  Figure  3- 5a. 


3.7.3.  Inner  Loop  Design. 

First  Inner  Loop. 

For  the  first  inner  loop  design,  L*  should  handle  the  uncertainty 
in  effective  plant  P^  of  (3.3-4),  with 

\ 

x1  - P.3-°P— . 2x19  * 130  . 39  . X . 40  of  (2.3-2b) . 

2 P30P2b0  * P2c0  2X1  * 2 1 

Again,  in  case  of  9 “0°  , the  designed  first  inner  loop  of 

2.3.1  can  be  used  as  the  first  inner  loop  L2Q  here.  However,  a small 
phase  margin  of  18°  was  actually  used,  as  shown  in  Figures  3-5b,7,10 
for  both  cases  A and  B,  with 


L20  ” 


0.168(1  + 2286) (1  + 18Q) (1  + 144Q) 


I1+%JT+  (^-)2]  f1+^-+  (ott>V(1+^>  (1+T^r) 


_L_)2’ 

'320' 


°-vfes  , 8 . \ 2 • 2 

2680  '2680' 


150' 


3200' 
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3.7.4.  Second  and  Higher  Inner  Loop  (Case  A) . 

2 

The  second  inner  loop  L--  for  handling  the  effective  plant  P^ 


30 

of  (3.5-6)  is  analogously  derived  as  in  Chapter  2 for  0m2 1 


to 


satisfy  the  constraint  on  the  outer  loop.  The  values  of  X2 


X"?  of  P^  in  (3.5-6)  are  now 


P P + P 
2 _ *40*3b  3c 

*3  p 

J *30 


lxia  + 60 


39 


( 3. 7-4a,b) 


X3P30P2b  * P2c 
P20 


39x2x18  + 120 


- 381 


The  bounds  on  L30  and  the  designed  L30  (also  allowing  18° 
phase  margin)  are  shown  in  figure  3-6a,  with 


L30  " 


0.068(1  + 1440Q)  U + 2250  + (2250^  1(1'*'7000) 


11  2300  '2300}  1 11  2250  '2250'  J 1 14000  14000 

(3.7-5) 


0.66  , 8 . 2 , 
2250  '2250' 


0.6s 


2,2 


(1  + s/23000) 

Similarly,  the  bounds  and  loop  transmission  of  the  3rd  inner  loop 
(i.e.  i-3  ) are  shown  in  Figure  3-6b  with 


0-068(1  +-—-)  [l  + (ttttf)2}  (1  + s 


t40  ll*fi5^*,n5oo>2"1* 


75200' 


11745  '11745' 


36540 


120000 


11745  '11745' 


73000  ‘73000 


2 2 
) J 


(3.7-6) 

is  designed  to  handle  the  uncertainty  of  the  effective  plant  P^ 


in  (3.5-10)  with  X*  - 40  , X^-  390  , X*  - 3540 


' Ljo 


obtained  from  (3.5-11,12)  for  i-3 


2(Ui  . 


The  Bode  plots  of  all  the  loop  transmissions  for  both  the  detailed 
design  and  the  perspective  design  (explained  later)  of  case  A are  shown 


y>r&> 


( *♦  yj 


I I I I iff?*/:  !\  i»$-|h|  ‘ 

: UT  7ne*i(-fa)J  : ! : !./  : / : I .J>  f.  S: 

: : ; • / : • : >:  B: 

’:c«* ; r'l'jfe37-; i7"'f 7 I 


23.08 


-46.15 


Figure  3-6a.  Bounds  and  loop  of  L. 


-23.08 


. i -210* 


-*46.15 


Figure  3-6b.  Bounds  and  loop  of  . 


in  Figure  3-7.  Figures  3-8a,b  show  the  sensor  noise  effects  of  the 


different  designs  in  logarithmic  and  arithmetic  scales  respectively 


vHl4J 


Desert,  lirspect'* 
Detailed  Detijn. 


Figure  3-7.  Bode  plots  of  loop  transmissions  for  perspective  and 
detailed  design. 


3.7.5.  Second  and  Higher  Inner  Loop  (Case  B) 


are  used.  It  was  shown 


in  Section  3.5  that  if  there  exists  the  constraint  of  maintaining  some 


value,  then  if  0 ^0,  the  design  of  the  i-th  inner  loop  to 


constraint,  is  analogous  to  the  design  of  the 


satisfy  this 


first  inner  loop  in  Chapter  2 to  satisfy  the  bounds  on  the  outer  loop 


(i.e.  i-1  ) of  case  A was  purposely  designed 


The  first  inner  loop  L, 


with  0 


so  it  could  be  used  here  (however,  in  case  A,  there  was 
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no  consideration  of  0^  in  designing  L30  ).  Actually,  for  i*  1,2,3  , 
- 40  ■ U°  * U°  = 39  , the  problem  becomes  that  of  Figure  3-4b  from 
Equations  (3.3-2),  (3.5-8),  (3.5-13)  for  each  i . The  bounds  obtained 
on  L^q  , L40  and  the  loop  transmissions  are  shown  in  Figure  3-9a,b 
respectively,  with 
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The  Bode  plots  for  both  the  detailed  and  perspective  designs  are  shown 
in  Figure  3-10.  Figures  3-lla,b  show  the  sensor  noise  effects  for  these 
designs  in  logarithmic  and  arithmetic  scales  respect  vely. 


3.8  Design  Perspective  Procedure 

As  in  Chapters  1,2  Design  Perspective  enables  the  designer  after 
picking  a specific  structure  and  specific  trade-off  values,  to  very 
rapidly  obtain  a good  approximation  to  the  actual  final  detailed  design 
for  that  specific  structure  and  trade-off  values. 

1)  The  first  step  is  to  design  the  single-loop  Lg0  to  handle  the 

entire  uncertainty  of  P3  (=100.5  db)  in  (3.2-5a).  Let  aQ  be  the 

trade-off  (overdesign)  margin  of  LQ  . Then  the  outer  loop  L^Q  of  a 
multiple-loop  design  is  obtained  by  shifting  the  UHF  characteristic  of 
Ls0  uPward  by  4g  , the  net  saving  due  to  the  multiple-loop  design.  In 

our  example,  aQ  = 0 , so  $g  = (l^-U^db  « 58  db  . Here,  the  L10  used 


5^310 


Figure  3-9b.  Bounds  and  loop  of  L 
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Detailed  Design- 


(at  which  Il2qI  has  its  P63*  value)  by  placing  a 
of  Figure  l-24a  over  Lin  of  Figure  3-7,10  and 


2)  Find  a) 


noting  the  arrc 


with  value 


of  Mj  ( “-10  db  ) obtained  from  Appendix  3 of  Chapter  2.  Point  Q2 
is  |L  ( ju>  ) | . Draw  a horizontal  line  in  Figure  3-7,10  at  A 


Figure  3- 11a.  Sensor  noise  effects  in  decibels. 
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Figure  3- lib.  Sensor  noise  effects  in  arithmetic  scale 
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magnitude,  with  A.  - Min(B.  _)  - a,  “ -35  db  for  the  example,  where 
2 ru  l 

Min(B.  „)  is  the  point  of  minimum  value  of  B.  _ . Here  A'  “ -74  db 
n2  n2  2 

for  | LJo I is  also  done,  because  LiQ  = 0 for  i i 3 . Lay  L^q(IF) 
of  Figure  l-24a  over  Figures  3-7,10,  so  that  Q coincides  with  Q . 

Find  the  intersection  point  C 2 (or  ) of  this  L^ (IF)  with  the 
A2  (or  Aj  ) horizontal  line  in  Figures  3-7,10. 

4)  Pick  the  L^Q(HF)  curve  in  Figure  l-24b,  according  to  the 

0^2  value  (<*18°  here)  being  used  for  L^q  . Lay  LiQ(HF)  on  Figures 
3-7,10,  such  that  C coincides  with  C2  (or  C2  ).  L2Q  consists  of 
LiO (HF)  of  Figure  l-24a  in  the  intermediate  u>  range  and  of  L^q(HF) 
of  Figure  l-24b  in  the  high  frequency  range.  Use  the  pattern  of  L.,g(HF) 
to  the  left  of  C to  obtain  a smooth  continuous  curve  for  L2Q  . 

5)  Steps  3,4  are  repeated  in  order  to  determine  L3Q  . Use  the 

arrow  on  L^q(HF)  of  Figure  l-24b  to  locate  u>m3  , according  to  the 

8m2  value  being  used  ( 0m2=o°  for  case  A,  18°  for  case  B) . Use 

(2. A. 3-6)  to  obtain  , giving  a new  point  Q3  in  Figure  3-7,10. 

Then  lay  L^q (IF)  of  Figure  l-24a  on  Figures  3-7,10,  so  that  Q3  = Q . 

A horizontal  line  of  value  A.  = Min(B.  ,)  - a,  is  drawn,  etc.  In  the 

3 n3  3 

example,  A3  at  -35  db  for  |l3Q|  40(3  A^  at  -55  db  for  |l*q|  • Min(Bh3) 

is  the  point  of  minimum  value  of  B^3  . The  entire  process  is  repeated 
until  all  the  loops  are  exhausted.  The  dashed  curves  in  Figures  3-7,10 
are  the  results  of  using  this  perspective  technique.  They  are  in 
excellent  agreement  with  the  actual  detailed  design  (solid  lines  in 
Figures  3-7,10),  obtained  by  using  the  computer  to  obtain  the  inner  loop 
bounds,  finding  a rational  L^Q(s)  to  satisfy  these  bounds,  etc. 
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6)  After  each  LiQ  is  obtained,  it  is  a good  idea  to  sketch  the 

effective  P values  to  use  for  the  sensor  noise  effect.  Thus,  after 

L _ is  obtained,  sketch  P,_  of  (3.2-5)  in  Figures  3-7,10.  If  there 
sO  10 

is  little  sensor  noise  amplification  ( | Lao^Pio ^ n0t  ^arqe  over  a lar9® 
u range) , there  may  be  no  point  in  using  more  feedback  loops.  It  is 
also  easy  to  see  the  saving  in  sensor  noise  effect,  by  using 

ILiq/PiqI  • sketch  |p2Q|  to  see  the  hf  N2  sensor  noise  effect 
(“|l20/P20|  ).  Similarly  |l30/P30I  gives  the  hf  N3  sensor  noise 
effect. 

I 

The  designer  has  to  decide  which  sensor  points  to  use  and  the 
corresponding  trade-off  values.  Design  perspective  anables  him  to 

try  out  various  designs  very  quickly  and  easily,  and  thus  arrive  at  a 
suitable  trade-off,  after  which  he  can  proceed  with  a careful  detailed 
design. 


J 
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CHAPTER  4 

THE  BASIC  PARALLEL-CASCADE  STRUCTURE 
4.1  Introduction 

This  chapter  extends  "quantitative  synthesis"  to  the  plant 
structure  (heavy  line)  of  Figure  4-la,  consisting  of  two  parallel 
branches,  each  with  one  internal  sensing  point  C^2  (i-1,2).  The 

system  input  r and  output  C can  be  sensed,  giving  4 degrees  of 
freedom,  in  the  form  of  4 independent  data  processing  of  these  measure- 
ments. It  is  assumed,  as  in  previous  chapters,  that  each  plant  section 
has  independently  uncertain  parameters,  with  known  ranges.  The 
tolerances  2. 2-1, 2 are  to  be  satisfied  V p^  . Emphasis  is  on 

the  minimization  of  net  sensor  noise  effects  at  the  plant  input  x . 
Again,  this  problem  obviously  does  not  lend  itself  to  a rigourous 
mathematical  theorem-proving  treatment.  So  again,  the  approach  taken 
is  to  find  the  principal  design  factors  and  trade-offs.  Several  design 
philosophies  are  presented,  and  except  for  one  of  these,  "Design 
Perspective"  (Section  4.5)  is  provided  by  relatively  simple  fast 
sketches,  based  on  only  a preliminary  single-loop  design  for  the 
problem.  This  enables  the  designer  to  decide  early  in  the  game,  which 
loops  to  use  and  how,  and  gives  him  even  then  a very  good  idea  of  the 
final  design.  This  chapter  paves  the  way  for  the  much  more  complex 
plant  structure  of  Chapter  5. 
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4.2  Design  of  the  Outer  Loop 

The  objectives  of  the  problem  are  the  same  as  Equations  (2.2-1), 

(2.2-2).  They  are  also  achievable  [H3)  with  a single  loop  Gj  “ Gs  » 

G,  - 0 , G_„  “ 0 in  Figure  4-1.  But  the  resulting  L *G  P may  then 
12  22  s 3 

require  very  large  bandwidth,  causing  great  amplification  of  sensor 
noise  , as  in  the  previous  chapters.  In  Figure  4-la,  let 

a cl s)  L.  (s) 

p & p P +p  p . t(S)  - ' --  * F Is)  — = 

P P12P11  22  21  ' ' R(s)  1 + L.  (s) 


PG, 


Lx(s) 


1 + P12G12  + P22G22 


P G. 

e 1 


(4.2-la-d) 


-X 


V 


"l  1 + ^12^12  * ^22*^22  + WL1 


In  (4.2-ld)  the  sensor  noise  effect  is  examined  at  the  plant  input 
X where  it  tends  to  be  large  (see  Figure  4-8) , causing  plant  saturation. 
As  shown  in  Chapters  1,2,  the  hf  range  is  the  major  trouble  source,  so 
the  major  effort  in  sensor  noise  effect  reduction  will  be  made  in  hf. 

It  was  also  shown  in  Section  1.4.4  that  the  saving  in  the  outer  loop  1^ 
(by  using  the  inner  loops)  is  very  minor  for  u>  < , so  the  same 

design  philosophy  of  Chapters  1 and  2 is  used  here:  Let  the  outer  loop 
from  C (of  Figure  4- la)  cope  with  the  uncertainty  of  the  entire  plant 
P up  to  , but  handle  only  the  uncertainty  in  P^  , Pj^  for 

^ wxl  * giv®s  40  much  more  economical  than  in  a single- 
loop design  in  which  L1  ■ Lg  must  cope  with  the  entire  P for  all  u>  . 

Outer  loop  design  in  the  hf  range,  is  now  a single-loop  problem  with  the 


equivalent  plant  P of  (4.2-lc),  denoted  by 

pi  A ? (P12°P“  V.  , L,  £ pX  , L 


* 1 +P120G12  + P220G22 


■ P*  *G 
10  eO  1 


(4.2-2a-c) 


Figure  4-la.  The  basic  parallel-cascaded  structure. 
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Figure  4-lb.  Specified  time-domain  Figure  4-lc.  "Equivalent" 
bounds  on  step  response.  frequency-domain  bounds. 


equivalently  a nominal  L1Q  so  that  (2. 2-1,2)  are  satisfied.  The  single 
loop  design  was  reviewed  in  Section  1.2,  and  is  not  repeated  here. 

4.3  Design  of  the  Inner  Loops 

4.3.1.  Existence  of  the  Inner  Loop. 

In  order  to  obtain  the  most  economical  outer-loop  and  thereby 
minimize  the  effect  of  sensor  noise  at  X , it  was  assumed  that  the 

local  feedback  G12  , G22  in  Figure  4-1  are  able  to  handle  the  un- 
certainties of  P^2  , P22  . It  is  next  proven  that  such  G12  ' G22 


so  in  hf, 


®11  +el2  “ e21  + e22 

, let 

1192  ' g22=b21g2  ' 

Gi2“ 

<klllt12  +k21k22)gl 

1 1 + <k12bn  +k22b21)92 


Compare  with 


A (bllk12  *b21k22)9l 
lx  1 + <k12bn  +k22b21)92 


*i2  ell+e12 

' G1  " gls 

(4.3-1) 

(4.3-2) 

(4.3-3) 


Arg  « Arg  Llx  , but  IHML1xI  V kll,k21  uncertainty. 

(4.3-4) 

For  u>  > the  problem  of  violation  of  specifications  is  of  at 

the  top  of  the  template  (recall  Figures  1-Ala, b) . Suppose  there  exists 
a G2  which  satisfies  the  system  specifications  V k12'k22  uncertainty 
at  fixed  k^-b^  * k2l“b21  * <nie  above  (4*3~4)  guarantees  that  any 
other  values  of  k^  , k21  give  points  lower  down  in  the  template,  in 


.1 

i s 

I 


I 


J 


I 


the  Nichole  chart  — further  from  the  locus.  Hence,  such  a g2  is 

completely  satisfactory.  So  now  we  need  only  consider  the  problem  at 
k^-b^  , *21-b21  . at  which  - Llx  of  (4.3-3). 


Let  S ■ "’jlkU*bJlk22>*2 


Li.-sr/u,s' 


(4.3-5) 


(4.3-6) 


and  the  problem  is  to  find  a nominal  LgQ  , 3 Lg/(1+Lg)  satisfies 


certain  bounds,  with 


• X, 


(b12bll  ~>‘b22b21) 
max  " (a12bu  ♦ a22b21> 


(4.3-7) 


This  is  precisely  the  same  problem  as  in  a 2-section  cascade  design  of 
Section  3.4,  using  the  structure  of  Figure  4-2,  with  Lg  ■ PgGg  * 

eg 

PB  “ ^/s  , G0  ■ ggS  . The  problem  there  is  identical  to  the 

present  problem.  Hence,  cascade  design  philosophy  with  its  Desgin 
Perspective  can  be  used  to  find  Lg  here,  from  which  g2  is  obtained, 
and  then  g^2  , g22  are  available  from  (4.3-1). 


Figure  4-2.  A 2-section  cascade  system. 


4.3.2.  Trade-off  Function  x(u)  between  . 

The  above  provides  a simple  design  procedure  but  is  based  on 


(4.3-1),  giving  G12//G22  " bll//b21  * 1*'is  is  not  desirable  for 

example  ~ b2l  wh*le  N22>>N12  in  Fi9ure  Return  to  (4.2-1) 


and  let 


a Pi2bil 

L2  ^ P12G12  + P22G22  ' °i2  “ P12bn  + P22b21 


(4.3-8) 


The  effect  of  sensor  noise  Ni2  , at  the  plant  input  X , is 

^X_  ^i2 Gi2  °i2L2 

i2  - Ni2  " l + P12G12tP22G22  + PG1  " <1+L2>  (1+I.j)  " Pi2 


for  i«l,2  . Hence,  to  minimize  this  effect  reduce  |g^2|  as  fast  as 
possible.  But  I*2  must  cope  with  the  uncertainty  in  P^2  , P22  , 
ignored  by  . So  there  is  a trade-off  between  G^2  , G22  . Introduce 


a trade-off  function  x(w)  by  defining 


b21g12  " bll922  “ (b11/b22)x(“) 


Equation  (4.3-1)  corresponds  to  x*0  . 


4.3.3.  Design  of  First  Inner  Loop. 


(4.3-10) 


Substitute  (4.3-10)  into  (4.2-lc),  replacing  P^G^  by  b^g^  , 
P21G1  by  b2igi  inasmuch  as  the  top  of  the  template  in  U)>U)X^  need 


only  be  considered,  giving 


_ b21P22egl 
1 1+P22eg22 


, p . b»‘klf u * V . ,4.3-U..b) 

22e  x^12bll + b21b22 


To  design  g22  , assume  x(o>)  can  handle  the  uncertainty  in  k12  , 
so  replace  k12  by  k12Q  - a12  and  define 


and  define 
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„o 

Li  S p Q 
22  22e’22 


b22(a!2bll  * b21k22)922 
Xa12bll  + b21b22 


(4.3-12) 


with  L°  - L°  at  k__  ■ a_. 
220  22  22  22 


Thus  L22  need  handle  only  the  un- 


certainty in  P-Joe  * *-n  which  only  k22  is  uncertain.  From  (4.3-llb) 
if  x*0  , then  L22  must  handle  the  uncertainty  in 

)C12bll  + k22b21 


which  gives  the  design  approach  of  Section  4.3.1. 


The  ratio  of  maximum  uncertainties,  x = 0 cf.  x/0  , is 
jc-0  ( bllb12  * b21b22  /(  al2bll  * b22b21  bllb12  * b21b22 


( P11P12  P21P22  ^ /(  1211  P22P21  \ = °11°12  T °21u22 
^blla12  +b21a22  V ^ al2bll + a22b21  ^ bllal2  + b21b22 


(4.3-13) 


In  view  of  (4.3-lla),  the  loop  L22  of  4.3-12,  is  designed  analogous  to 
the  inner  loop  design  of  a 2-section  cascade  system,  with  P of  Figure 

O 

4-3  analogous  to  P22e  of  4.3-12. 

4.3.4.  Design  of  Second  Inner  Loop. 

If  x(w)  /0  is  used  for  the  design  of  L22(x)  in  Section  4.3.3, 
then  x(w)  must  handle  the  uncertainty  in  k^2  ignored  by  L°2  . Re- 
write (4.3-11,12)  as  follows: 


bllk12 


)/(-W  ■ 


L22  " 922P22e 
(4. 3-14a,b) 


which  describes  precisely  the  structure  of  Figure  4-3.  Thus,  the  design 
procedure  for  the  loop 


k12bllX 
X b21b22 


(4.3-15) 


In  Figure  4-1 , at  the  nominal  values  of  the  P 


Substituting  4.3-10,11,14  and  manipulate  to  give  at  hf  (at  which  the 
sensor  noise  effects  dominate) , 


Similarly,  -X/N^  can  be  manipulated  (at  the  nominal  values  of  the 
P, . ) , at  hf  into 
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If  x - 0 


xO 


0 , so  (4.3-18)  gives  the  extra  effect  of  N 


12 


which  one  pays  for  easing  the  design  of  L^  by  means  of  x^O  , i.e. 
by  the  factor  given  in  (4.3-13). 


-X/N 


12 


'11 


-X/N 


(case  x *■  0 ) 


(4.3-19) 


22 


21 


but  I _X/N22  lx-0  > l"X/NooLwn  from  (4.3-9.  13,  17). 


22 1 x*0 


4.4  Design  Perspective 

As  emphasized  in  previous  chapters,  "Design  Perspective"  provides 
the  designer  an  easy  tool  for  deciding  which  sensor  points  to  use,  and 
the  trade-offs.  Since  it  is  so  useful,  it  is  again  presented  here  in 
detail,  although  the  procedures  are  very  similar  to  that  of  the  previous 
chapters. 


(1)  The  first  step  is  to  design  the  single-loop  LsQ  to  handle 


the  entire  uncertainty  of  P in  (4.2-la).  Let  be  the  overdesign 


margin  of  L 


10 


Then  the  outer  loop  L^q  of  a multiple-loop  design  is 


obtained  by  shifting  the  UHF  characteristic  of  Lg0  upward  by  db  , 
the  net  saving  due  to  the  multiple  loop  design. 


P P + P P ■» 
llx  12x  21x  22x  1 


♦ — 20  loq  < — 

s y I P P + P P 1 

l iivrnfi  r5i»rimJ 


a,  db 


llx  120  21x  220 

for  the  structure  of  Figure  4-la  ~ 36  db  in  the  example  of  4.5. 

(2)  Find  wm22  which  L 22Q  has  its  peak  value)  by  placing  a 
transparency  of  Loutar  of  Figure  l-24a  over  L._  of  Figure  4-5a  , and 


noting  the  arrow. 


(3)  Locate  Q22  (the  peak  position  of  L220  ^ at  wm22  ' w^th 

value  of  M22  (“  -14  db)  obtained  from  Figure  1-22,  because  L22Q 
design  is  analogous  to  the  inner  loop  design  of  cascade  system.  Point 
Q22  is  lL220^gm22^  * Draw  a horizontal  line  in  Figure  4-5a  at  A22 
magnitude,  with  A22  = Min(Bh22)  “ a22  ~42*-J-4*“24  db  for  different 

designs  in  the  example  of  Section  4.5),  where  Min(B^22)  is  the  point 
of  minimum  value  of  Bh22  • Lay  L^Q(IF)  l~24a  over  Figure 

4-5a,  so  that  Q coincides  with  Q22  . Find  the  intersection  point 
C22  of  this  L^q (IF)  with  A22  horizontal  line  in  Figure  4-5a. 

(4)  Pick  the  L^Q(IF)  curve  in  Figure  l-24b,  according  to  the 
®m22  value  being  used  for  L22Q  . Lay  LiQ(HF)  on  Figure  4-5a,  such 
that  C coincides  with  C22  . L22o  consists  of  LiQ(IF)  of  Figure 
l-24a  in  the  intermediate  u>  range,  and  of  LiQ(HF)  of  Figure  l-24b 
in  the  high  frequency  range.  Use  the  pattern  of  L^(HF)  to  the  left 
of  C to  obtain  a smooth  continuous  curve  for  L^Q  . 

(5)  Steps  (3), (4)  are  repeated  in  order  to  determine  LxQ  . Use 
the  arrow  on  L^Q (HF)  of  Figure  l-24b  to  locate  w . The  peak  value 
Mx  of  LxQ  is  obtained  from  Appendix  3 of  Chapter  2,  because  LxQ 
basically  is  the  first  inner  loop  design  of  structure  in  Chapter  2. 

Determine  Qx  by  and  M . Then  lay  L^q (IF)  of  Figure  l-24a 

\ 

over  Figure  4-5a,  so  that  Qx  = Q . A horizontal  line  of  value 
Ax  - Min(Bhx)  - <*x  (t*  -42  db  in  the  example)  is  drawn,  etc.  Figure  4-5a 

are  the  results  of  the  design  example  in  Section  4.5  by  using  this  fast 
perspective  technique.  They  are  in  excellent  agreement  with  the  actual 
detail  design  (in  Figure  4-5b) , obtained  by  using  the  computer  to  obtain 
the  inner  loop  bounds,  finding  a rational  LiQ(s)  bo  satisfy  these 
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bounds  , etc . 

(6)  After  each  LiQ  (i.e.  L22Q  , Lx0  ) is  obtained,  it  is  a good 
idea  to  sketch  the  effective  P values  to  use  for  the  sensor  noise 
effects.  Thus  after  Lg0  is  obtained,  sketch 

|P0|  = lPiioP120  +P210P22C)I  in  Fi9ure  4-5a*  If  there  is  little  sensor 
noise  amplification  ( |Ls0/pol  not  large  over  a large  w range) , there 
may  be  no  point  in  using  more  feedback  loops.  After  has  been 

obtained,  it  is  easy  to  see  the  saving  in  sensor  noise  effect,  by 

comparing  the  area  between  | l-gg/Pg  I and  | L10^P0  I ' sltetch 

lP22ol  = lP120(Pllx/P21x)  +P22ol  (=  1,1/s2  in  the  examPle>  to  see  the 

hf  N22  sensor  noise  effect  (<*  I L22</  ^P120  ^Pllx^P21x^  +P220^I  ^ • 

2 

Similarly  sketch  |p12Q|  in  the  examPle)  to  9et  lLxO/,pi2ol  ' 

then  see  the  hf  N^2  sensor  noise  effect  from  (4.3-15b). 

This  fast  design  perspective  enables  the  designer  to  try  out 
various  designs  very  quickly  and  easily,  and  thus  arrive  at  a suitable 
trade-off,  after  which  he  can  proceed  with  a careful  detailed  design. 


4.5  Design  Example 

The  structure  is  shown  in  Figure  4-la,  with 


P12  = [l,90]/s2  , Pn  = [1, 10]/s 
P22  = [0.1,2]/s2  , P21  = [1, 10]/s 


(4.6-1) 


The  time-domain  bounds  on  the  acceptable 


translation  into  bounds  on  |T(jw) | are  shown 

L 


disturbance  specification  is 


1+L 


£2.3  db. 


step  response  and  its 
in  Fioures  4-lb,c.  The 


4.5.1.  Outer  Loop  Design. 

For  ui  t ui^  , we  let  the  outer  loop  handle  the  entire  plant 

of  P in  (4.2-la),  with  uncertainty  ^ " 836.36  . But 

for  w i w , , L.  handles  only  the  uncertainty  U°  in  P*  of  (4.2-2a) 

XI  1 * e 

with  U°  =*  * o'Txi  = i0  * The  savin9  *c  thereby  obtained  by  a 

multiple  loop  design,  is 

4>g  - 20  log [83Yq36  1 - c»1  db  *=  (38.45  -a1)db  <*  -36  db  . 

The  design  of  is  exactly  the  problem  of  single-loop  design.  The 

resulting  is  shown  in  Figure  4-4  on  the  Nichols  chart,  and  in 

Figure  4-5b  its  Bode  plot. 
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Figure  4-4.  Outer  loop  L, _ on  Nichols  chart 
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Figure  4-5a.  Design  Perspective. 


Figure  4- 5b.  Detailed  Design. 
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4.5.2.  Inner  Loop  Designs. 

So  far,  the  only  two  design  techniques  given  are  those  of  x*0 
of  4.3.1  and  x 1 0 of  4. 3. 2-5.  If  |n22|  » |Nn|  , and  since 
bjj“b21“10  , then  Equations  4.3-13,18,19  indicate  that  x^O  should 
be  used.  However,  if  ~ N22  ' t*len  x = 0 would  be  the  reasonable 

approach.  In  both  cases,  the  design  L220  ana^°9°us  to  that  of  the 
first  inner  loop  of  a cascade  system,  with  effective  inner  loop  plant 
of  (4.3-12)  for  w > , 


P22e(x*0) 


P22xlP120Pllx  * P21xP22] 
xP120Pllx  + P21xP22x 


If  x » 0 is  used. 


^226 (x-0) 


_P22xtP12PHx~>'P21xP22] 


P P 

*21x*22x 


The  uncertainty  of  P 
is 


22e  (xj*0) 
90x10  + 10x2 


22e (x=0) 


10  + 1 


is 

83.63 


1X10  + 10x2 
1X10  + 10x0.1 


2.7  , but  that  of 


The  loop  transmissions  and  the  bounds  on  L220(x^0)  for  trade“off  value 
a22  = 5 db  , 15  db  , and  that  on  L220(x=0)  are  shown  Figures  4-6a,b,c 
respectively.  This  completes  the  design  for  the  case  x=0  because 
g12/,g22  ■ bn/b21  • F°r  x / 0 , the  design  of  LxQ  is  analogous  to 
the  first  inner  loop  design  of  Chapter  2.  The  loop  transmission  and 
bounds  of  the  corresponding  LxQ  (for  a22  ■ 5 db  , 15  db  ) are  shown 
in  Figures  4-7a,b.  (In  this  example,  the  bounds  on  0m22  wer®  not 
maintained  in  the  course  of  finding  the  bounds  on  L .)  Bode  plots 
are  shown  in  Figure  4-5b,  including  perspective  design  results  in 
Figure  4-5a  for  comparison.  The  rational  functions  of  , L22f)  , 

I*x0  are  given  in  Appendix  1 of  this  chapter. 
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Sensor  Noise  Effects. 

These  are  shown  in  Figure  4-8.  If  x«0  is  used,  |TN12|  » 
|-X/Ni2|  » |-X/N22|  « | tN22 I which  are  significantly  iess  than 
lTN12lxyo  ' as  e*P«cted.  However,  It^I^  « |tn22|^0  over  most  of 
the  frequency  range.  One  decides  which  is  best  by  evaluating 
lTNij  Xij  ( ju) | , etc.  But  it  is  obvious  that  there  are  N22  , N^2 
functions  for  which  x-0  is  best,  and  others  for  which  x ? 0 is  best. 
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P11  ' P21  ^°r  W > Wxl  * The  inner  ^00PS  handle  the  uncertainty  of 

P12  ' P22  u > “xl  ’ was  a*so  shown  Section  4.3.1  that  the 

inner  loops  may  be  designed  at  fixed  (largest)  values  with  = 

b^g^  , P2iG2l  = p21g21  ' Th^-S  retained  in  the  present  section. 

The  new  approach  of  this  section  is  explained  by  means  of  Figure  4-9. 

At  b11  , b21  , the  outer  loop  is  denoted  by  Llx  , and  sketched  in 
Figure  4-9.  This  is  at  fixed  nominal  values  (a12  , a22)  of  , 

k22  . Due  to  the  uncertainties  k^  € [a^b^]  , k22  £ [a22,b22]  , 

L^x  will  vary  but  g, 2 , g22  must  be  designed  so  that  the  resulting 
Llx  does  not  penetrate  into  y^  of  Figure  4-9.  We  simply  partition 
this  reserve  (between  Llx  and  ) » between  L^2  and  L22  . 

U1V1W U4V4W4  show  four  different  partitions.  The  first  inner  loop 

is  designed  such  that  due  to  its  plant  uncertainty  at  fixed  second  inner 
loop,  L^x  does  not  penetrate  beyond  . Then  the  second  inner 

loop  is  designed  so  that  due  to  its  uncertainties  in  both  inner  loops 
L^x  does  not  penetrate  into  y^  (and  if  there  is  also  a constraint  y 
on  the  first  inner  loop,  so  that  this  constraint  is  also  not  violated) . 
One  may  choose  either  (12)  or  (22)  as  the  first  inner  loop.  The 
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L220/P220 


22 


(1+L220) (1+L120) (1+L10> 


<4.6-2c) 


From  (4.6-2b,c),  to  reduce  effect 


Li20  ' 1-1,2  * 


Ni2 


, minimize  the  bandwidth  of 


4.6.2.  First  Inner  Loop  Design. 

To  minimize  the  bandwidth  of  L^2Q  ' P22*P220  used*  i,e* 
L120  ^andles  the  uncertainty  in  P^2  only.  Let 
P, 


12 

'l2  “ p ' 

A *120 


(4.6-3) 


p p 4.  p p 

.1  £ *12  ilx  * *220*2  lx  „ & „ 


(4.6-4) 


1 ^ X12P120Pllx  * P220P21x 
6 (1 + P220G22) (1 + X12L120) 


(4.6-5) 


Equation  (4.6-5)  is  analogous  to  (2.3-2),  with  Pj^o  here  ~ pa0  of 
(2.3-2),  Pllx  here  ~ pb  there,  P22oP21x  here  ~ pc  there,  and  X^2 
here  ~ there.  The  design  of  L^20  therefore  identical  to  that 
of  L10  in  Section  2.3.1,  except  that  the  constraint  on  L*  here  is 
the  curve  * instead  of  there. 


4.6.3.  Second  Inner  Loop  Design. 

The  next  step  is  to  design  G22  to  handle  the  uncertain  plant 
P22  ' wh:*-ch  was  ignored  in  the  design  of  L220  * Let 


22 


P22/P220 


(4.6-6) 


..  .2  A P21xP22  *P12Pllx 

*"  L1  " 1 +P..G^  +P,,G,.  G1 


£ -2 


22  22  12  12 
P"  * G.  . 


X22P21xP220  * X12P120Pllx 
1 + X22L220  + X12 (1+L220) L120  1 


(4.6-7) 
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L220  *s  des^9ne<^  so  that  does  not  penetrate  into  over  the 

range  of  X12  , X22  . 

4.6.4.  Design  Example. 

The  same  structure  and  specifications  in  Section  4.5  are  used. 
The  extra  margin  of  the  outer  loop  below  is  increased  by 

2 db,  because  of  the  present  approach.  The  resulting  Llx  {at 
k^  » b^  , *21  = ^21  ^ s^own  Fi9ure  4-9.  L^2  is  taken  as  the 

first  inner  loop.  Figures  4-10a-d  show  the  bounds  on  Lj^o  and  t^ie 
designed  L^2Q  for  the  4 different  cases  in  Figure  4-9.  The 

corresponding  bounds  on  L220  and  result^-n9  L220  are  s^own 
Figures  4-lla-d  (in  which  the  bounds  on  6,^2  were  not  considered) . 
Since  the  bounds  on  L22o  for  cases  1*3  are  very  close  to  those  of 
cases  2,4  respectively,  the  same  L22Q  used  for  cases  1,2,  and 
another  one  for  cases  3,4.  Bode  plots  of  L^x  , L^2Q  , 1*220  ^or  a*^ 
these  r»<?es  are  shown  in  Figure  4-12.  The  rational  function  values 
are  9-.A..1  in  Appendix  1 of  this  chapter.  The  sensor  noise  effects 


are  shown  in  Figure  4-13. 
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4.7  Partitioning  with  Reverse  Order 

In  Section  4.6,  the  inner  loops  via  G^2  , G22  were  defined  in 
Equation  (4.6-1),  such  that  L12Q  was  first  designed,  and  then  L22Q  . 
In  this  section  the  order  is  reversed,  so  the  loops  are  now  defined  by 


120 


P120G12 


L220  P220G22/ (1  + L120} 


(4. 7-la, b) 


The  noise  effects  at  plant  input  X at  nominal  plant  values,  due  to 
noise  sources  , N^2  , N22  are  now 


X 

N, 


L10/P0 


L220/P220 


1 +L. 


10 


N22  (1 +L220)  11  + ^.O* 


L120/P120 


N 


12 


(1+L120)(1+L220)  (1+110) 


(4.7-2a-c) 


^10  here  is  the  same  as  in  Section  4.5.1,  which  was  also  used  in 
Section  4.6.  The  design  procedures  for  L220  ' L120  are  exactly  the 
same  as  in  Sections  4. 6. 2, 3,  except  for  the  change  in  subscripts 
(12)  (22)  , (11)  (21)  . For  example,  Equations  (4. 6-4, 7)  become 


.1  A P22P21x  + P120Pllx  „ A 
L1  * 1 + P,  _G, . + P "G’“  G1  Pe°l 


120  12 


22  22 


P P + P P 

2 A llx*12  22  21x 
1 “ 1 + P,„G,„  + P.„G_  1 


12  12 


22  22 


X12PllxP120  + X22P220P21x 
1 + A12L120  + X22(1  +L120)L220 


G = P2  • G 
lei 


(4.7-3) 


respectively.  The  same  numerical  example  in  4.6.4  is  designed  in  this 
manner  for  comparison  with  4.6.4.  Figures  4-14a-d  show  the  bounds  on 
L220  ^es^9ned  first)  , with  the  resulting  loop  L22o  ^or  4 different 
UiViWi  cases  of  Pi9ure  4-9 . The  corresponding  bounds  on  L12Q  and 


Bounds  and  loop  L, 


Figure  4-15c.  Bounds  and  loop  L 
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result  ing  L120  are  shown  in  Figures  4-15a-c  (in  which  the  bounds  on 
Qn^  were  not  considered) . The  rational  function  expressions  are 
given  in  Appendix  1 of  this  chapter.  Bode-plots  are  shown  in  Figure 
4-16,  and  sensor  noise  effects  in  Figure  4-17. 

4.8  Discussion 

Comparison  of  the  noise  effects  in  Figures  4-8,13  and  17  for  the 
different  approaches  shows  that  the  best  cases  in  each  approach  are 
very  close.  The  small  differences  may  be  ascribed  to  the  departures 
from  optimal  design,  as  all  L _ Q were  shaped  by  hand.  In  the  first 
approach  with  x * 0 , the  fixed  relation  between  G^2  , G22  in 
Equation  4.3-10  does  not  permit  trade-off  between  | TN22 ^ ' ^12^  ’ 
Trade-off  is  possible  if  x^O  is  used,  by  means  of  a22  . Or,  by 
changing  the  order  and  designing  the  (12)  loop  first,  trade-off  is 
achievable  by  means  of  a^2  . The  above  approaches  have  the  benefit 
of  Design  Perspective.  The  final  Partitioning  approach  presents  another 
means  of  trade-off  between  the  two  loops.  It  appears  to  be  superior  to 
the  x / 0 , a22  or  a^2  trade-off  method  in  sensor  noise  effects  in 
the  example  above.  It  is  also  more  flexible,  because  the  partitioning 
(Figure  4-9)  may  be  shaped  so  as  to  change  the  trade-off  vs.  u>  . Thus, 
different  parts  of  u^viWi  are  critical  in  different  u)  ranges. 
Obviously  assignement  of  some  value  to  a22  or  <*12  is  a cruder  trade- 
off technique.  Of  course,  Design  Perspective  has  not  been  developed 
for  the  Partitioning  approach. 

It  must  be  admitted,  however,  that  a more  intensive  study  is 
needed  of  the  various  techniques  in  order  to  determine  the  best  approach 
for  any  specific  design  problem.  For  example,  in  the  numerical  design 
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problem  of  Section  4.5,  the  much  larger  level  of  P^2  vs.  P22 
undoubtedly  an  important  factor.  Nevertheless,  the  techniques 
presented,  especially  Design  Perspective,  enable  the  designer  to 
quickly  try  various  designs. 


is 
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The  prefilter  F and  the  outer  loop  for  all  the  different 

approaches  are: 

.2447  x 107(S2  + 3.226s  + 3.715) (1  + s/8) 

n * 5 2 2 s 

s[s  + 2.5682s  + 3.1705]  [s  + 36.73s  + 1645] * (1  + 

„ . 1861  x 105 (s4  + 29 . 17s3  + 217.38s2  + 764s  + 947.33) 

s7  + 778.16s6  + 81823.65s5  + .10123  x 10?s4  + .6277  x 107s2  + <p 

(0  = .16128  x 108s2  + .2467  x 10®s  + .1756  x 108 

The  transfer  functions  L_„  , L _ for  Sections  4.3  to  4.5  are: 

220  xO 


L220(a22“5)  =' 


0.645(1+ W) 

..  0.6s  s . 2.  ..  0.3s  . s >2,  , 2s  ^ s 

ll+~+(25)  1 [1+H[00^l00)  lll  + 400^400)  1 
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Case  2. 


0.9(1  ^>(1^) 


120  [i  0 • 3s.  + (-£.)  2 ] (i  + _Jt_)  [1  + 0 8s  + (-%—)  2 ] 2 

1 25  V25J  J V 420*  1 900  k900'  J 


L220  “ L220  °f  Case  1‘ 


Case  3. 


1.44(1  +y^Q)  (1  + 5§o) 


120 


[1  + — + (— )2]  (l  +— ) [l  + 0 ' 7s- + ( — 2 — )2]2 
1 22  '22'  M 340M  1250  V1250  J 


0.07 


220  [1+£J*  + (±)2][1+Oi3s  ji.)2j 

1 25  '25*  J 1 90  v90;  J 


Case  4. 


1*56(1  +I§o><l  + gfo* 


120 


[1  + — + (— )2]  (i  +-2— ) [i  + 0 • 7.s  + (—2 — )2]2 
1 22  V22'  M 340M  1350  V1350'  J 


L220  " L220  °f  CaS®  3< 

The  transfer  functions  l22q  ' L120  for  Section  4,7  are: 


Case  1. 


0.06 


220  g__)  2]  [i  + 0-3a  4.  (_g_)  2] 

1 32.5  V32. 5*  100  ^100;  J 


^ 0-645(1  ^43o)  d+np*  x (1  + 2fo} 

120  [i  (_®_) 2]  n +2i§®  + (_S_) 2j2 

1 27  '27*  1 1 800  '800*  J V 420^ 


Case  2. 


0.07 


220  [ 1 + Q--.8iL  + (_ g— ) 2 j [ 1 + 0 ».3.s. + (-2—)  2 ] 

1 32.5  V32. 5'  J 1 100  'lOO'  J 


°.9(1+-2_)(l+-2_) 


l20  [i +£i®i+  (-L.)2]  (i +_®_)  [i +£t§®.  + (_®_)2]2 
1 24  l24;  M 400M  855  ^855*  J 


Cases  3,4.  L 


0.07 


220  _s_2  0^3s  _s_2 

1 32.5  V32. 5'  J 1 110  V110'  J 


L120  " L120 


of  Case  2 
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CHAPTER  5 


THE  GENERAL  PARALLEL  CASCADED  PLANT  STRUCTURE 


5.1  Introduction 

Previous  chapters  presented  quantitative  synthesis  techniques  for 
the  structures  of  Figures  2-1,  3-1  and  4-la,  with  detailed  examples. 

The  final  objective  is  to  develop  a quantitative  synthesis 

technique  for  the  more  complex  and  general  structure  of  Figure  1-14, 

which  is  redrawn  here  as  Figure  5-la.  There  are  m parallel  branches. 

Branch  one  has  n^  cascaded  plant  sections,  etc.,  and  branch  m has 

m 

n cascaded  sections,  so  that  there  is  a total  of  1+  V (n.-l)  possible 
m . l 

1*1 

independent  feedback  loops.  Each  plant  P„  has  its  given  parameter 
uncertainties  with  the  known  ranges.  The  objectives,  as  before,  are  to 
satisfy  specifications  of  the  form  (2. 2-1, 2),  for  each  P„  „ . 

The  feedback  compensation  functions  are  to  be  chosen  to  satisfy  those 
objectives  despite  the  uncertainty,  but  to  do  so  with  minimum  net  effect 
of  the  sensor  noise  sources  N , N^ ^ at  the  plant  input  3 . It  is 
obvious  that  no  analytical  solution  is  available  for  this  highly  complex 
problem.  A transparent  design  theory  which  gives  the  designer  a great 
deal  of  insight  into  the  trade-offs  has  been  developed.  Furthermore,  as 
before,  this  insight  is  available  without  a detailed  design.  This  design 
theory  will  be  recognized  as  based  upon  the  work  of  previous  chapters, 
especially  Chapters  2,4. 

5.2  Design  of  Outer  Loop 

The  design  technique  is  developed  by  means  of  a specific  example 
in  Figure  5-la,  m*4  , n^*n2*2  , n3=n4=3  , giving  Figure  5-lb.  Let 
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then 


m p p p 

+ P,J 

43  42  41 

33 

S P4  3<*4  3 * 

P33G33 

" P43P42G42 

+ P 
*33 

P,  G, 

A C A p 

11 

3 2 11 


a Li 

- F -r-r- 


L1  - 1 ^ G1  • Pel'Gl 


VP1  . L1 

“ — in  hf. 


»i  i + h3  + h2  + i+h  p3 


(5.2-la-c) 


(5.2-2a-c) 


Hence,  to  reduce  the  effect  of  sensor  noise  source  on  plant  input 

X , minimize  the  bandwidth  of  . From  Equation  (5.2-2b),  should 
cope  with  the  uncertainty  of  Pel  . To  minimize  the  bandwidth  of  , 

minimize  the  template  of  Pgl  . It  was  shown  in  previous  chapters  that 

Gij 

the  major  noise  effect  is  at  hf  (high  frequency)  with  pij  * ki/5  • 

It  was  also  shown  that  the  saving  by  using  inner  loops,  is  minor  for 
ui  < o>x3  , the  frequency  at  which  L1Q  turns  the  corner  of  its  "universal 
high-frequency  boundary"  . So  the  design  philosophy  for  the  outer 

loop  of  previous  chapters  applies  here  again:  Design  the  outer  loop  L° 
with  the  uncertainty  of  entire  plant  P^  for  u>  < iox3  . For  u i wx3  , 
let  L°  handle  only  the  uncertainty  U°  in  P°^  , i.e.  of 

P° 

P°.  - , with  L°  - P°  • G 

el  jQo  1 el  1 

P1  “ P430P420P41  + P330P320P31  + P220P21  + P120P11  (5-2'3a-c) 

L10  • PelO  • G1  ’ p«10  ■ P10'©!  ' ■Q?  ■ WH30*H20  ’ 


This  gives  a much  more  economical  outer  loop  L1Q  than  that  (Lgg)  of 


-159- 


single-loop  design,  which  must  handle  the  entire  uncertainty  for 

all  u>  . Thus,  in  the  design  example  below  Lg0  must  handle  un- 
certainty of  54.5  db,  whereas  need  handle  uncertainty  U°  of 

5.2-3c  only,  which  is  20  db. 


Numerical  Example 
A)  Problem  statement 

(a)  System  structure:  In  Figure  5-la,  let  m-4  , , n^-n^-3  , 

giving  Figure  5- lb. 


(b)  Plant: 


k12 

n 

^11 

12 

“ S(1  +S/A) 

' P11 

S ' 

22 

D 

k21(l  +S/z) 

22 

S 

' P21 

2 * 
S(S  + BS  + C) 

ij 

.iu 

s 

* i * 

3,4  , j - 1,2,3 

(c)  Plant  Uncertainty:  ^n'^l  ^ (4 , 40] 

k12'k22  6 l25'75°l 

k31,k41  € 14,201 

k32,k42  6 15,401 
k33,k43  € [5,75] 

A 6 [1,2]  , z € [1,2]  , B € [0,1]  , 

C € [0.04,1]  . 


B)  Specifications 

(a)  Bounds  on  unit  step  response  C(t)  : The  same  as  shown  in 
Figure  4-lb. 

(b)  Bounds  on  | T ( ju»)  | : Shown  in  Figure  4-lc  were  derived  from 
time  domain  bounds  of  Figure  4-lb. 

L 

1+L 


(c)  Disturbance  Response: 


* 2.3  db 


r 
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The  resulting  L1Q  and  LgQ  in  the  Nichols  chart  are  shown  in 
Figure  5-2a,  their  Bode  plots  are  shown  in  Figure  5-2b.  The  rational 
functions  of  Ls0  , L^0  with  same  prefilter  F are: 

(1  +S/8) 


F1  - 


(1  + S/2)  (1  + S/3) 


2.2(1  + S/0.5)  (1  + S/100)  (1  + S/30) (1  + S/4200) 

Js0  ^2  2~T 

S (1  + S/20)  (1  + S/1000) (1  + S/45)  (1 + 0.5s/2800  + (S/2800)  V 


2.2(1 + S/0.5)  (1  + S/100) (1 + S/30) 


10  2 o o 

S (1  + S/20) (1  + 0.7S/115 + (S/115)  ) (1 + S/45)  (1 4 S/1000) 


The  noise  effects 


X 

X 

N. 

9 

N, 

1 

s 

1 

at  plant  input  X at 


II, III 


nominal  plant  values  due  to  noise  source  for  both  single-loop  and 

multiple-loop  designs  are  shown  in  Figure  5.5. 


5.3  Design  of  Inner  Loops 

5.3.1.  Design  of  the  First  Inner  Loops. 

The  next  step  is  obviously  to  allow  for  the  uncertainty  in  Pi2 
(i * 1,2, 3,4)  to  be  handled  by  the  loops  via  G^2  • For  this  purpose,  let 

P2  = P430P42P41  + P330P32P31  + P22P21  + P12P11 


H2  ~ P430P42G42  + P330P32G32  + P22G22  + P12G12 


H30  = P430G43  + P330G33 


(5.3-la-c) 


pS«i 


P2G1 


1 ltK30  + H2  + 


,o  4 4 

j’1+b3o'4); 


h 


20 


20  1 + H 


•i2 


N 


30 


P L°  /P° 
llxzO'  2x0 


i2 


1+H30  + H20+P10G1  (1  + L°0)  (1  + L°Q) 


(5. 3-2a-c) 


(5.3-3) 


-~U-  — - 
->rf- 




figure  5- 2a.  LgQ  , and  bounds. 
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(with  P°2Q  - P°2  of  Equation  5.3-5  at  Pi2  * Pi2Q  for  i € [1,4]  ). 


At  this  point,  we  can  use  any  of  the  techniques  developed  in 
Chapter  4.  We  choose,  because  of  its  simplicity,  the  method  of  Section 
4.3.1,  i.e.  let 


Gi2  " PilxG2  f°r  ^ “ 1,2, 3, 4 , 


(5.3-4) 


6 ' 1 A 

where  at  hf,  Pilj£G2  - (bi]/s  1 )G2  * bii^2  * Then 


. o A (P430P42P41x  + P330P32P31x  + P22P21x  + P12Pllx)G2 


1 + H 


30 


A „o 


“’SxVfiSo  “ Pe2G2  ’ 


(5.3-5) 


In  Chapter  4,  Section  4.3.1,  it  was  shown  how  it  suffices  to 
consider  the  first  group  of  inner  loops  at  the  maximum  values  (k^  - b^  , 

21  * b21>  °f  P11  ' P21 


k^,  *b_.)  of  P, , , P01  sections  of  Figure  4-la.  Obviously,  the 


situation  is  the  same  here,  so  we  need  to  examine  for  uncertainties 


in  Pi2  only  at  P^ » Pilx  (i.e  at  kil*bil)  for  i"l,2,3,4. 


Now  at  P^  * Pilx  becomes 


Lix- 


P2xG1 


P2xG1 


Pe2Gl  (L2/G2> 


(1+H30)  (1+L°)  <jQ°0(l  + L°)  1+L°  1+L°  1 


G1 
(5.3-6) 


Equation  (5.3-6)  is  analogous  to  (3.4-3)  of  the  2-section  cascaded 


system.  The  design  of  L2Q  is  therefore  identical  to  that  of  the  inner 


loop  of  the  cascaded  system.  The  bounds  on  the  nominal  inner  loop  L2Q  , 


such  that  those  on  the  outer  loop  are  satisfied,  are  precisely  of 

the  same  type  as  those  on  Lgg * PB0GB  *n  * so  that  those  on  LQ 

o * 

are  satisfied.  In  our  present  example,  the  bounds  on  L2Q  (L2Q)  and 


L20  ^L20^  are  sbown  Pi-9ures  5-3a,b  respectively,  with 


2 

J 
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* 

L20- 


0.59(1  + s/700) 


[1  +0.4s/90  + (s/90)2]  [1  + 0.7s/2000  + (s/2000)2]2 


20 


0.59(1  + s/500) 


[l  + 0.4s/90+  (s/90)  2]  [1  + 0.7S/1250  + (s/1250)2]2 


where  L^q  is  used  if  no  inner  loops  via  G33  , G43  are  used,  and  L^q 

* 

is  used  if  inner  loops  via  G33  , G^3  are  used.  The  Bode  plots  of  L2C)* 


20 


are  also  shown  in  Figure  5-2b.  The  noise  effects  at  plant  input 


X at  nominal  plant  values  for  using  L2Q  or  L2Q  are  shown  in  Figure 


5-5,  denoted  as 


X 

2 

X 

N12 

9 

II 

N12 

respectively. 


Ill 


5.3.2.  Design  of  Second  Group  of  Inner  Loops. 

The  function  of  second  group  of  inner  loops  via  G, _,G__, . . . ,G  _ 

±z  ZZ  n\z 

in  Figure  5-la  is  to  handle  the  uncertainties  in  P, _,P__, . . . ,P  _ . 

lz  zz  m2 

Therefore,  in  Equation  5.2-3b  replace  the  pi30  by  pi3  » i m . 

Again,  as  in  Section  4.3.1,  we  need  consider  the  uncertainties  of  P^3 
only  at  the  maximum  values  of  P^  , j < 3 , i.e.  at  P^x  (corresponding 
to  j - b^ j ).  Following  again,  because  of  its  greater  simplicity, 
the  design  philosophy  of  Chapter  4,  Section  4.3.1  (note  that  we  can 
switch  philosophies  at  any  stage) , let 


Gi3  - Pi2xPilxG3  for  i“3'4  of  pi9ure  5-lb. 


(5.3-7) 


Then  for  the  structure  in  Figure  5-lb,  the  "new"  outer  loop 
(originally  L*  in  5.3-2a)  is  now 


(P43P42xP41x  + P33P32xP31x  + P22xP21x  + P12xPllx) G1 


,2  A 

i ” i + <paA  "♦P„G„)  + (PJ1„P. 


43  43  33  33' 


43‘42xG42  + P33P32xG32  * P22xG22  + P12xG12} 


from  5.3-7 
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2 (P43P42xP41x*P33P32xP31x'*'P22xP21x+P12xP11x)G1 

^ = 1+ (P43P42xP41x+P33P32xP31x) G3+ (P43P42xP41x+P33P32xP31x+P22  Alx*P12xPU^ °2 


Pe2  • G2 


L2  G1 


(‘*P«2'G2>  ^ 


(5.3-8) 


T1  A 1 „ 

L2  " Pe2  * G2 


(5.3-9a,b) 


?1  A P43P42xP41x  * P33P32xP31x  + P22xP21x  + P12xPllx 
®2  1 + (P43P42xP41x  + P33P32xP31x)G3 


Compare  (5.3-8)  with  (5.3-6).  If  LX  (with  effective  plant  PX_  ) 

< el 

satisfies  the  bounds  on  L°  , then  L2  also  satisfies  the  bounds  of 
L^x  . The  problem  is  to  choose  G3  to  accomplish  this.  Let 

L3  * lP43P42xP41x  + P33P32xP31x) G3  “ P3x°3  = Pe3  * G3 


P = (p  p p +PPP) 
e30  1 430  42x  41x  330  32x  31x' 


(5.3-10a,b) 


We  must  choose  G3  (or  L3Q  = Pe3QG3  ) such  that  {L^}  c (L°>  , 
V P43>P33  . Let  X3  - P°3/P°3Q  , and  L3  = Pe3G3  , then 


'2  ' Pe2V  [ 


Pe3  + (P22xP21x  * P12xPllx) 
1 + Pe3G3 


X3Pe30  + (P22xP21x  * P12xPllx) 
1 + X3L30 


(5.3-11) 


X Gi3  Pi2xPilxL30/P3x0  (5,3“12) 

■'  1 — = . ■ — — — ■ ■ — — — as  - ■ - — - 

Si3  1 + H3q  + H20  + P10G1  + **30^  ^ + ^20^  + **10^ 


With  P«30*P.3  o£5.3-10a,t  P<3-P430  , p33  - p330  • 


Equation  (5.3-11)  is  precisely  the  same  as  Equation  (2.3-2) 
of  the  first  inner  loop  design  of  triangular  structure  in  Chapter  2, 


which  is  redrawn  here  as  Figure  5-4a.  The  design  of  L 


is  therefore  identical  to  that  of  the  first  inner  loop  design  of  the 


triangular  structure  in  Chapter  2.  The  bounds  on  the  nominal  L 


such  that  those  on  L! 


and  its  bounds  are  shown  in  Figure 


5-4b,  with 


[1  + 0.8S/1000  + (s/1000)  ] [1  + 0.2s/3000  + (s/3000)  ] 


is  also  shown  in  Figure  5-2b.  The  noise  effect 


at  plant  input  X at  nominal  plant  value  due  to  N 


is  shown  in 


5.4  Design  Completion 


The  design  for  the  structure  of  Figure  5-lb  is  completed  as 


follows.  The  previous  steps  gave  the  loop  transmissions  L 


From  (5.3-10),  G,  is  available  from  L,./P 


From  (5.3-7) 


G,  are  found.  With  G 


known,  H 


is  known  from  (5.3-lc).  G is  then  obtained  from  (5.3-5) 


from  known  P 


Substitute  G into  equation  (5.3-4)  to 


for  i*  1,2, 3, 4 . G,  is  obtained  from  (5.2-3c) 


Time-Domain  Simulation 


The  actual  step  responses  of  the  numerical  example  for  various 


combinations  of  parameter  values  are  shown  in  Figures  5-6a-d.  The 


disturbance  responses  — actually  response  to  steps  at  the  output  side 


of  the  prefilter,  are  shown  in  Figure  5-7a-c 


Ill  CO*V  IhUOUS  SVSfF.H  »OOKLl*ft  Pfcl»0*»4*  Ilf  V|«3  fUFCUllUN  UylPUl  III 


-170- 


•»*** 

009 

III 


9 - t 
•a  a 

2 F«- 

3 9 « • • 
l »cm; 


9a«\**»>99n.'9a>» 
*-  C9-*A£--*A 

- - t«rf 


* 


9009 
» • • • 

* M 4 1 

«rv««  ■>9  — /**-*'*'/'**'^.."*  • 

• !«>•  *«.0-9y>>94©*'>lt*/»9^>»9 

9>/'?  — c «• 

so  ••  • c***4c«nj>>»?» 

a 


OA*>m 
0«t" 
fc  V * 9 
??■*> 


9 •**»*• 
O-VJ* 

« r c * 
X * >9 
>J?7 


9 9*1 

a-* 

9 00 
909 


II 


*»xiv 
coo 
I I • 


9 »•» 

• 99 
9 «• 


-»^9«9r4%9‘0>.'*©90*4^^00-'.''9-*9^N^94tf>-V*99 


V 

o 


VI 

n 


90  90 
0009 
«***»>• 
• • • « 


to 

ei 


5 • 

IC 


zxrz 

3333 

& a aa 

Mil 

Ch«« 


* a 


5% 

a 9 


5°. 

a a 


a 

a 

u»* 


• 

1 

ft. 

♦ a3: 
• a 3 3 
a 33 

33 

♦ 

a a as  3 33 300 
333 

i 

o° 

a 

O 

o 

o 

a 

3 

O 

33 

' 

c 

a 

♦ 

3 

3 

a 

♦ 

O 

a 

♦ 

O 

M 

3 
• 
a 
c ♦ 

a 

3 

a 

39 

a 

3* 

3-> 

A -I 
B=l 
C=l 


Lb 

LA 

■C 

'O 

t 

* 

X 

0 

ft. 

41 

4o 

4- 

4? 

*» 

*5 

i5 

75> 

ft. 

4 

4o 

4 

AD 

ft 

>5 

o5 

P> 

75o 

ft. 

4 

JO 

4 

JO 

& 

5 

5 

4o 

AD 

%>} 

75 

75 

75 

75 

ft. 

«3D 

4 

<*i  O 

ft* 

? 

S' 

&D 

40 

?5 

15 

15 

75 

:<  «.,.*p».*;:r%\i»(ri.‘!  k»*fj»''c»,>f":»:\;i.,C3ceo 
SJ'C*  1N0*)J*3M1».'3»>MmJJCJ3OO 
• • t<cf  fjf  .*?7f  r • • • • • 


OOOC OO  6 70000 00900000 OC  090000000959000909 

c :o::c:ccooocoucooooooo;oc*>or,oo':  rooococc 

W 000000009CC9090&0990000C0030999C9COC  5900 

i ccocccccoo-^raco***?  c-*v". ••  9*-«ec*-N*;cr«*-cac 
-o-N-i#r  «nc  v>  ••••_••••  • 


Figure  5-6a.  Simulation  results  of  step  response. 
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Figure  5-6c.  Simulation  results  of  step  response. 
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Figure  5-6d.  Simulation  results  of  step  response 


s Ifriti*  t 1 1»4  f*  S»Cflf»N-j  A d»»A*»C*«tS  SYS! 


Figure  5- 7a.  Simulation  results  of  disturbance  response 
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Figure  5- 7b.  Simulation  results  of  disturbance  response 
(same  parameter  set  as  Figure  5-6c) . 


Sf«M4.#rm*i  of  sicri'tN-j  or  4 m»a»> c^s  srSf 


176 


• n^NTi^-o?  - * - r n # c v r c t c 4 rr  4-  r « c c * 

j>/\f  ,.s>  >?o  j»)‘>:3  j3::o  jj: 

r — e:?a.‘.\'C30cc:Tm?j‘C90POJcjo:?.*f?JOs. ooec  .*.  ■> 

30  

I * — • • • • • • • • • 


4 


oe  900 
1 1 • • 1 

-».»*fc.fcocM-;5C  ‘V'vjck  -«a  ff-»  — 

s — *♦  — #c-o«-4^^»x^r»oM^^4o«V4/»^*o  — ®c  00000  oc  oooocooooo  : o 90  5 

re  • •••  •— ctre?e??c:  • • • • • • • » 


Figure  5-7c.  Simulation  results  of  disturbance  r 
(same  parameter  set  as  Figure  5-6d) 
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5.5  Design  Perspective 

It  is  very  important  that  one  can  get  good  perspective  of  the 
design,  without  performing  the  detailed  design.  Since  the  design 
philosophy  of  each  loop,  as  shown  in  previous  sections,  is  the  same  as 
either  cascade  (in  Chapter  1)  problem  or  the  problem  of  Chapter  2,  the 
design  perspectives  in  those  two  chapters  are  applicable.  Figure  5-8 
shows  those  loops  obtained  from  design  perspective,  they  are  in  excellent 
agreement  with  those  obtained  from  detailed  design  shown  in  Figure  5-2b. 


5.6  Summary  of  Design  Procedure  for  General  Structure,  Based  on 
the  Design  Philosophy  of  Chapter  4,  Section  4.3.1. 

(1)  Design  the  outer  loop  as  in  the  previous  chapters. 

(2)  Choose  G12  - pllxG2  * Gi2  “ PilxG2  ' etc*  Defi-ne  L2  in 
a manner  similar  to  (5.3-5).  The  design  of  L°Q  is  the  same  as  that 
of  the  first  inner  loop  in  the  cascade  problem,  for  which  the  solution 

is  known.  Similarly,  choose  Gi3  - Pi2xPilxG3  411(1  define  L°  as 
5.3-10  (in  the  example,  the  process  terminated  there,  so  was  denoted  as 
L3  ) . The  design  of  L°0  ,is  equivalent  to  that  of  the  second  inner 
loop  in  the  cascade  design  problem,  for  which  the  solution  is  known. 

This  process  continues  up  to  n , the  minimum  value  of  n.  , 
i-l,2,...,m  in  Figure  5-la,  i.e.  until  at  least  one  of  the  m branches 
has  been  exhausted. 


(3)  At  this  point  the  cascade  design  philosophy  is  no  longer 
valid,  but  by  continuing  in  the  same  manner  the  design  philosophy  of 
Figure  2-1  (or  Figure  5-4a)  is  henceforth  valid  step  by  step,  until  the 
last  step  at  the  maximum  value  of  n^  of  the  n^  , i*l,2,...,m  . 
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Suppose  for  example  in  Figure  5- lb,  there  existed  an  additional  section 
P34  , permitting  feedback  to  X via  G34  . The  in  (5.3-9)  will 

then  be 


.1  [P43P42xP41x  + P340P33P32xP31x  + P22xP21x  + P12xPllx3  . 

■•o  ■ " , „ r r : r r — - r . „ 1 G„ 


r 

2 

/ 1 4.  p ft  \ 

r ( (P43P42xP41x  + P340P33P32xP31x)G3  ' 

L 

(1  P340G34} 

l 1 * P340G34  J 

(5.6-1) 


PP  P +P  PP  P p 

r o A *4 3*4 2x  4 lx  340*33  32x*31x  „ A 3x  „ A o 


3 1 * P340G34 

*°3  - ^30  ' then 


G-  G, 

3 3 

*30 


e3  3 

(5.6-2) 


.1  Pe3  * (P22xP21x  + P12xPllx)/^: 


1 + Pe3G3 


V3Pe30  + (P22xP21x  + P12xPllx)/c^30 


1 + X3L30 


(5.6-3) 


P P iP  /vP 

X ^ i2x  ilx^c/^xO 

“i3  (1  . L°0) (1 4 L°0) (1  . L°0) 


(5.6-4) 


The  design  of  L30  is  exactly  the  same  as  the  first  inner  loop  design 
of  Chapter  2,  as  discussed  in  Section  5.3.2,  except  the  definition  L° 
of  (5.6-2)  is  different  from  L3  of  (5.3-11). 


For  the  next  step  of  allowing  the  uncertainty  in  the  additional 


plant  P34  , let 

G34  = 

P33xP32xP31xG4 

(5.6-5) 

P3x 

P43P42xP41x  + P34P33P32xP31x 

Pe3  * P3V«0lo 

Sx" 

P34P33xP32xP31x  ' Pe4  “ P4x/^40 

* * P°.G. 

4 e4  4 

X3  Pe3/Pe30 


P°  /p°  = P /P 

e4'  e40  34'  340 


(5.6-6a-g) 
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then 


L2  « P3x  * (P22xP21x  * P12xPllx)  Q 

2 1 * Y^O  * 2 

_ X3Pe30  * (P22xP21x  * P12xPllx) 
i - a4l40  ♦ 


(5.6-7a,b) 


P P P L /P 
i3x  i2x  ilx  40/  4x0 


N.  . 

i4 


(1+L40)(1  + L30)(1+L20)  {1  + L°10) 


(5.6-8) 


Equation  (5.6-7)  is  analogous  to  (2.3-5).  So  the  design  of  L^q  is 
exactly  the  same  as  the  second  inner  loop  design  of  Chapter  2. 

5.7  Summary  of  Design  Procedure  for  General  Structure,  Based  on 
the  Design  Philosophy  of  Chapter  4,  Sections  4. 3. 2-4. 

The  various  techniques  of  Chapter  4 can  be  used  until  one  of  m 
branches  is  exhausted.  But  thereafter,  only  the  triangular  technique 
of  Chapter  2 can  be  used.  Design  perspective  allows  the  designer  to 
quickly  decide  which  philosophy  to  use  in  the  first  part.  For  example, 
if  N42»  N12'N22'N32  ' one  wants  to  try  the  Philosophy  of  Sections 
4. 3.2-4  at  the  very  first  inner  loop  stage.  Let 


b41912  - 

bllg42  ' 

(bll/b42)y12 

b41g22  “ 

b2lg42  = 

(b21/b42)y22 

(5.7-la-c) 

b41g32  * 

b31g42  = 

(b31/b42)y32 

Equation 

(4.3-10) . 

Then  L^  in  (5.3-2a) 

can  be  written  as 

l}  

P2G1 

P2G1 

(5.7-2) 

l30Tn2  *^20+H2 


.1 


J 


i 


I 


Suppose  n12  “ N22  “ N32  ' let  yi2  “ y “ y for  i_1'2'3  » ^ 
substitute  (5.7-1)  into  the  denominator  of  (5.7-2), 

<®20  + H2  “<^20  + a43k42g42  + a33  bT7  [ b31942  + bTT  Y ] 


20  2 «*20  43  42^42  33  b. 


f 

I 


1 I 

1 

i1 

1 

C 


k22 

b2i 

*12 

r.  . bn  „ I 

+ b41 

- 2lg42+b42  Y . 

[bng42+b^yJ 

{b41b42^20  + y[a33k32b31 +k22b21 + k12bll] 


41  42 


+ b42g42  ta43k42b41  + a33k32b31  + k22b21  + k12bll] } 


“ b b {b41b42^20 + yK2 +b42g42la43k42b41 +K2}  } 
41  42 

Replace  K2  by  its  nominal  value  K2Q  , since  y will  handle  the  un- 
certainty in  K2  , so 


_ 1 . , b41P42egl 

Vy) 3 

1 + P42eg42 


(5.7-3) 


where 


r1  ■ ,5-7-4a’bi 

b41b42&20  + y‘K20 


with  K20-K2  at  k.2-ai2  , i = 1,2,3  . 

In  view  of  (5.7-3),  the  loop  L°2  of  (5.7-4b)  is  designed 
analogously  to  the  inner  loop  design  of  a 2-section  cascade  system  with 
P0  of  Figure  4-3  analogous  to  P°2e  (5.7-4a) . 

The  function  y(oi)  is  next  designed  to  handle  the  uncertainty 
in  K2  ignored  by  L°2  , such  that 
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and  define  L 


b41b42^20 


(5.7-7) 


Thus,  the  design  procedure  for  the  loop  Ly  in  (5.7-7),  such  that 
(5.7-5)  is  satisfied,  is  precisely  the  same  as  that  of  the  first  inner 

loop  in  the  triangular  structure  in  Chapter  2.  The  design  perspective 
there  also  applies. 

The  procedure  is  the  same  if  n42  « ^ , instead 

of  n42  >>  n12  “ n22  “ N32  ' excePfc  n°«  take  y12 * 0 in  (5.7-la),  and 
K2  * a33k32b31  + k22b21  * 


Rewrite 


,1  A b42ta43k42b41  + K2> 


b b 
41°42’ 


:®20  * V 


a43k42 


^20  _ b41(02O 


4142*-*20 


(5.7-5) 


(5.7-6) 
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